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ABSTRACT 


Local mass transfer due to the impingement of a 
confined laminar two-dimensional air jet on a flat surface 
has been studied. The influence of the jet Reynolds number 
and the jet-to-plate spacing on the local mass transfer were 
investigated. The range of Reynolds number was 100 to 400 
with the slot width, b, taken as the characteristic length. 
The jet-to-plate spacings were 2b, 4b and i2b. Experimental 
mass transfer studies were made using a swollen polymer 
method coupled with laser holography interferometric 
techniques. 

The local Sherwood number along the impingement piate 
was found to exhibit a local minimum and a local maximum in 
the region away from the jet centre. The locations of the 
extrema points were a function of the jet Reynolds number 
and jet-to-plate spacing. 

A two-dimensional numerical study was also made 
Simulating the experimental set-up. The momentum and 
transport equations were numerically solved using hybrid 
differencing schemes (upstream-weighted and upstream 
schemes). The numerical study confirmed the presence of the 
local extrema in the Sherwood number. Contours of the flow 
stream function indicated that the presence of the Sherwood 
number extrema is due to flow recirculation in the region 


between the confinement plate and the impingement surface. 
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Excellent agreement was obtained between the measured local 
Sherwood numbers and those computed using the numerical 


solution of the transport equations. 
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1. INTRODUCTION 


Jet impingement flows are frequently used fer their 
excellent heat and mass transfer characteristics, where 
localized and controlled surface transfer is desirable. The 
drying of textiles, veneer, paper or film material, the 
annealing of metal and Oidstac Sheets, the tempering of 
glass, and the cooling of gas turbine blades and miniature 
electronic components are some of its more important 
practical applications. 

Heat and mass transfer characteristics of the various 
types of impinging jets have been studied rather 
extensively. Under technically realistic conditions, these 
Studies have been mainly concerned with relatively high 
velocities so that the flow, developing from the exit of the 
nozzle from which the jet issues is turbulent. The case of 
heat and mass transfer in impinging flow with a jet that is 
laminar at the nozzle exit has recieved lesser attention. 
This is especially true for the case of two-dimensional jet. 
Within the knowledge of the author, there has been only one 
publication dealing with the laminar two-dimensional 
impinging jet in the presence of the confinement plates. Van 
Heiningen et. al. (96) predicted numerically the flow field 
and impingement heat transfer due to a laminar 
two-dimensional jet with an upper confinement plate, 


including the effect of uniform suction at the impingement 
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plate. No experimental results of either heat or mass 
transfer due to a confined laminar impingina two-dimensional 
jet are available. 

Reasonably accurate analogy exists between heat and 
mass transfer provided that the mass transfer rate is low 
end Che iormal, surface velocity is nearly zero. In such a 
case, a single transport process, either heat or mass 
transfer, can be studied in isolation. Usually determination 
of local heat transfer coefficients is not very reliable 
Owing to the relatively large errors involved in the heat 
transfer sensors. Therefore, it iS more convenient to resort 
to mass transfer studies rather than heat transfer 
experimentation. According to the analogy for stagnation 
flows, it is possible to use the usual (Pr/Sc)°-* factor to 
convert the fairly precise mass transfer measurements 
presented by the Sherwood number, into heat transfer 
measurements presented by the Nusselt number. 

The present work is concerned with a confined 
two-dimensional jet having an initially laminar fully 
developed profile at the nozzle exit. Impingement mass 
transfer due to this jet has been studied both 
experimentally and numerically in order to emphasize not 
only the engineering applications but also the macroscopic 
mature of the transport phenomenon. A set of empirical 
equations which are obtained from the mass transfer 
experiments can be used immediately for engineering designs. 


On the other hand, numerical predictions verified by the 
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experimental results provide a complete understanding of the 


flow and mass transfer characteristics. 


* wan 
i : oe) 
; a 
; > & ’ 
Py oie 


os is eats ab arbres "S, afta 
aa | cee¥etaatoabg sry 


fi 7 
ae ‘= Se a » het ee 
> ss pes 3 a Z mE io a 


1.1 OBJECTIVES OF THE WORK 


In 


this work it was desired to: 


examine experimentally the effects of the jet Reynolds 
number and jet-to-plate spacing on the local mass 
transfer due to the impingement of a confined laminar 
two-dimensional jet, | 

study the effect of the presence of confinement plates 
by comparing the experimental results from this work to 
those obtained by other investigators using unconfined 
two-dimensional jets, 

develop a numerical model which would predict the flow 
field and the mass transfer due to a confined laminar 
two-dimensional jet, and finally 

investigate numerically the effect of nozzle exit 
velocity profile, (i.e. flat and parabolic velocity 
profiles), on the flow and mass transfer characteristics 


of a confined laminar impinging two-dimensional jet. 
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1.2 CHAPTER CONTENTS 


A brief review of related literature in impinging jet 
flow, heat and mass transfer is given in Chapter 2. The 
choice of experimental techniques is discussed in Chapter 3 
with experimental set-up and procedure given in Chapter 4. 
Calibration and validity of the experimental set-up are 
given in Chapter 5. The theoretical development which leads 
to the finite difference equations is presented in Chapters 
6 and 7. The adopted numerical procedure is given in Chapter 
8 and the validity of the numerical solutions are tested in 
Chapter 9. Both experimental and numerical results are 
discussed in Chapter 10 where the effects of jet Reynolds 
number, jet-to-plate spacing, nozzle exit velocity profile 
and presence of confinement plates are studied. 

Finally, conclusions are drawn in Chapter 11 and 
recommendations for further study are outlined in Chapter 
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2. LITERATURE REVIEW 


Detailed investigations of the flow field, heat and 
mass transfer due to single as well as multiple impinging 
jets have been given by numerous authors. Because of the 
large number of possible governing parameters, it is hardly 
Surprising that a eapucgte ges between the results and 
correlations of different investigators. A complete list of 
literature, together with suitable editorial comment, are of 
themselves a major undertaking. Therefore, no attempt will 
be made to review all literature in this area in detail. 
Only the more relevant references dealing with a single 
laminar impinging jet similar to the one studied in this 
work will be reviewed in detail. 

In general, there are two types of impinging jets: 
te® Liquid@jete (or iree jet): liguid to gas, i.e. water to 

air referred to as "water jet". 
2. submerged jet: liquid to liquid or gas to gas, i.e. air 
to air referred to as "air jet". 
The two types of jets differ substantially from each other. 
For a liquid jet, entrainment is negligible and the jet 
forms a free surface at ambient pressure imposed by the 
surrounding gas. This is the reason why it is also called a 
"free jet". For a submerged jet, entrainment is important 
and there is a substantial amount of mixing between the jet 


and the surrounding fluid. 
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There are only few studies of impinging jets that deal 
with liquid jets (11,58,63,98). On the other hand, impinging 
submerged jets are studied in more detail. In particular, 
turbulent submerged jets are studied in detail due to their 
wide industrial applications. Although literature survey of 
turbulent impinging submerged jets is not the main objective 
in this chapter, a brief review will also be given here. 

For a single turbulent impinging submerged jet, 
numerous investigations have been made for both axisymmetric 
ela el Oy i 20 ied ne 324,31 33730, o4,06,67 ,68,70,86587,95) 
and two-dimensional (6,9,10,18,19,42,101,102) cases. 
Different variables such as jet flow rate, size of the 
nozzle and jet-to-plate spacing have been considered. For 
more detail on the studies of single turbulent impinging 
submerged jet, the reader is referred to the literature 
reviews given by Cartwright and Russell (10), Gauntner et. 
al. (21) and Martin (48). 

For multiple turbulent impinging submerged jets, 
numerous investigations have been made for both axisymmetric 
(20,33,41,59) and two-dimensional (19,42) cases. An 
additional variable to those mentioned above for a single 
submerged jet is the spacings between adjacent jets. For 
more detail on the studies of multiple turbulent impinging 
submerged jets, the reader is referred to the literature 
review given by Martin (48). 

Another interesting topic in the study of turbulent 


impinging submerged jets is the impingement heat transfer 
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with crossflow. A few investigators have studied the 
impingement heat transfer with crossflow for Singlets, 69), 
and multiple (17,40,56,57) axisymmetric submerged jets. 
Their studies were concerned with the effect of jet-to-plate 
Spacing, jet flow rate and controlled crossflow flow rate on 
the impingement heat transfer. 

Turbulent submerged jets have been studied rather 
extensively as mentioned above. However, laminar submerged 
jets have recieved lesser attention. Most studies of laminar 


submerged jets deal with a single jet. 


2.1 FLOW CHARACTERISTICS OF A SINGLE LAMINAR IMPINGING 


SUBMERGED JET 


The flow pattern produced by a single submerged jet 
(laminar or turbulent, confined or unconfined) impinging 
normally on a flat plate can generally be subdivided into 
three characteristic regions: the free jet region, the 
Stagnation flow region and the wall jet region (15,21,48). 
The flow field of an unconfined laminar impinging submerged 
jet is shown schematically in Figure 2.1. The flow field of 
a confined laminar impinging submerged jet is very similar 
to the one shown in Figure 2.1 and will be discussed later 


in this Section. 
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FIGURE 2.1 : FLOW FIELD OF AN UNCONFINED LAMINAR 
IMPINGING SUBMERGED JET 
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2.1.1 UNCONFINED SUBMERGED JET 
Most of the studies of a laminar impinging submerged 
jet deal with an unconfined jet. In this case, no 


confinement plates of any kind are present. 


2.1.1.1 FREE JET REGION 

The free jet region is defined as a regicn of the 
submerged jet not influenced by ere impingement plate. A 
submerged jet has been examined at some length (27.10) temnn 
general, near the nozzle exit, the jet is decelerated by 
tangential shear stress. At the same time, the surrounding 
Still fluid is accelerated producing a "mixing region". The 
width of the mixing region increases continuously, and at 
some distance downstream it is wide enough to have 
penetrated to the centerline of the jet. Up to this point 
the centerline velocity is practically unaffected by mixing 
and is substantially equal to the nozzle exit velocity. 
Beyond the end of the so called "potential core” the 
centerline velocity decays as the jet shares its momentum 
with more entrained Copan 

Extensive reviews of literature concerning laminar 
submerged jets were presented by Schlichting (78). He 
studied the decay of centerline velocity and the spread of 
the jet for both axisymmetric and two-dimensional cases by 
applying the boundary layer theory. For the sake of 
simplicity, Schlichting assumed the jet exit is infinitely 


small and the velocity at the exit is infinity in order to 
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retain a finite volumetric flow rate as well as a finite 
momentum. For a laminar two-dimensional submerged jet, the 
Gecayrot) center lines velocity ,. vj x=s,;0and the spread»of the 


jet, ¢, are given by 


(v/ Vi) ixse = fot iy/b aed? (25 4) 
to/ bye cen. (y/b)2 ~? es :- 02:, 2)} 


where b is the slot width, and d, and e, are constants. 
Their values depend on the flux of momentum at the jet 
nozzle exit. For a laminar axisymmetric Submerged jet, the 
decay of centerline velocity and the spread of the jet are 


given by 


(W/ FVipse = da (y/a)~' (2e2)) 
Y/d =e, (y/d) (2.4) 


where r is the radial coordinate, d is the nozzle diameter, 
and d, and e, are constants. Their values are dependent on 
the flux of momentum at the jet nozzle exit. 

The length of potential core has been investigated by 
several researchers for jets of finite size. The core length 
in general is defined as the distance from the nozzle exit 
to a position where the centerline velocity is 98% of its 
initial value, (v/ Vj) |x=0=0.98 (22 enor thescententine 
pressure head is 95% of its initial value, (v/ Vi)? [x=0=0.95 


(48). Core lengths of about 4b for turbulent two-dimensional 
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submerged jets (48,101), or 4d for turbulent axisymmetric 
submerged jets (15,22,48) are to be expected. For laminar 
submerged jets, core lengths are expected to be longer and 
depend upon the jet Reynolds number. Since the behavior of 
turbulence increases the rate of mixing between the jet and 
Surrounding fluid, the core length decreases sharply. 
According to Hrycak et. al.'s (32)-experiment for an 
axisymmetric air jet, the core length is. about 15d at Re gy 


=500 and increases with jet Reynolds number to a maximum 


value of about 20d at Re ,=1000. 


2.1.1.2 STAGNATION FLOW REGION 

As the submerged jet approaches the impingement plate, 
the axial velocity component, v, is decelerated and 
transformed into an accelerated streamwise velocity, u. At 
the stagnation point, the velocity 1s zero and the pressure 
attains a local maximum. Stagnation flow of this type is a 
typical boundary layer flow with the influence of viscosity 
being restricted to a thin layer near the impingement plate. 

The axial extent of the stagnation flow region can be 
measured as the distance from the impingement plate where 
the axial centerline velocity drops to 98% of that in the 
undisturbed submerged jet at the same distance from the jet 
nozzle exit (22). In other words, it is where the submerged 
jet flow deflection begins. For a turbulent axisymmetric 
submerged jet, the beginning of stagnation flow region is 


about 1.2d away from the impingement plate (18,22,32). 
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Within the knowledge of the author, no such measurement has 
been reported for a laminar submerged jet. 

Schlichting (78) studied analytically the stagnation 
regions of laminar axisymmetric and two-dimensional flows 
against an infinite plate. The velocity components of the 
inviscid flow within the stagnation flow region are linearly 
proportional to the distance from the stagnation point. For 
two-dimensional stagnation flow, these velocity components 


are given by 


Vy me Ve (ieee) / 1b (eo) 


J 
Won a, (x/b) (25.6) 
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where a, 1S a constant value. For axisymmetric stagnation 


flow, the velocity components are given by 
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where a, is constant value. The viscous boundary layer 
thickness, o,, is defined as the distance from the 
impingement plate where the streamwise velocity, u, reaches 
99% of the value of the inviscid flow in either Equations 
2.6 or 2.8 depending on the geometry of the stagnation flow. 


The boundary thickness is given by Schlichting (78) as 


Ay SIT ENS 0 Oe (2.9) 
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for two-dimensional stagnation flow, and 


Tevdne = 19957 Via Y Reeyoe* (2590) 


for axisymmetric stagnation flow. 


2.1.1.3 WALL JET REGION 

The wall jet is the boundary layer flow formed by 
deflection of the submerged jet through ee Stagnation flow 
region. Due to the exchange of momentum with the surrounding 
fluid, the fluid which is accelerated in the stagnation flow 
region must eventually decelerate in the wall jet region. 
Therefore, the streamwise velocity, u, initially increases 
linearly in the stagnation region mentioned above must reach 
a maximum value at a certain distance downstream and finally 
decreases in the wall jet region. For a laminar axisymmetric 
jet, this maximum value is located at about 1d away from the 


Stagnation point for Re,=1. As the jet Reynolds number 


a= 
increases, this location moves further away from the 
Stagnation point (45). 

The wall jet has the characteristic of zero velocity at 
both the impingement plate and the outer jet edge, and hence 
it exhibits a maximum velocity. The flow pattern in the wall 
jet region can be divided into two parts according to the 
location of the maximum velocity: an inner layer which has 


features common to the ordinary boundary layer and an outer 


layer which has features common to a submerged jet. 
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Important parameters in the analysis of flow characteristics 
of wall jets are the growth of the wall jet boundary layer 
and the decay of the maximum velocity. The boundary layer 
thickness of the wall jet, o, generally used is defined as 
the thickness in wall jet outer layer where the streamwise 
velocity is 50% of the maximum velocity. 

The first complete theoretical analysis was performed 
by Glauert (25) who eeuaied both laminar and turbulent, 
two-dimensional and axisymmetric wall jets. Since then 
numerous investigations of turbulent (7,10,14,25,48,60,61, 
70,83,84) and laminar (25,71,79,82) wall jets have been 
reported. For a laminar two-dimensional wall jet, the decay 


of maximum velocity, u and the growth of the wall jet 


max’ 


boundary layer, o, are given by (25,82) 


deena (iy) cofer(z/>)o a6 (peta) 
a/b = g, (x/b)°-75 (2m) 


where f£, and g, are constants. For a laminar axisymmetric 
wall jet, the decay of the maximum velocity and the growth 


of the wall jet boundary layer are given by (25,71,79) 


i as Viut—et oc hE/0) ake” C2 eal Sa 


C/ Dee Geer (On ant (2.512) 


where £, and g. are constants. 


Uo? Bite hy ape -E ery hs fe Te Fe a \ oe 
Lh ; . ; A ‘ ul Aut s 7 » 
) 7 i oy oe ks f ) , f m 
ae oe ee 
‘ ' 


eeu, eeaiiovi 0h hae a ng ‘ae coset 
aed repaiad St totSe Lay umenbnath a: da. 


Bla 2 Iga es bsauyuitasaane . sai sty thaw ods: te 


oeintaaias ah aqadmy FNAL: sao rrva sisw. mh .2 
; phhsa ter musinam ads 19 aoe ag 
Bago lyse. .2eé.; ag hane fsais aa gay ae malgnce, 2 a hea 7 
Sasa te eh Lee Jit athe, Goyhie ue oh 
Aad ane |. ata The9 4) cy ing e, Ers Lurcus 
0 he, GR Bes SEO Sig heals Se arb isorreavnk, 
at ave 24st bias USS Ver. T Gane shh irens Sas i. 8, ee 
(soS5 of3 ., 190 ul eee sengs Pfhasit-awt satin: % +63 «+d ~ 
ioe Dlaew Ade 26 daigie ans a6. . i etlootey muse tied 
or aa | ‘aa yg ee 
($8, es) . Navi ots > , e¢et v7) 
pe 4 | ; . \ of 
i i : ® F y ) , 
. { :® 
AN 6 RM Ae IN ly 
is ..£) TaN)! 1 = CVS 
wanes «t 
oi FIR ER ERD tad tae A 702 +. Snes en09 - st .6 bie 4% 
aninig adit Bae gi ie1S4 womixan ef) t2 yuceh an? ed tls 
oo 
HET PE al 4s genie e746 toysl gtedeved sot finw + 
“4 ai 
= ’ ] ; re 
- ss Oi ; »a€e > 
_ - aut | (6\\2) 3 = # \ nuite” ‘UN a ‘ 
7 - — . = 
{ H #+. (Bye) we de 


16 


2.1.2 CONFINED SUBMERGED JET 

The confined two-dimensional submerged jet used in this 
study 1S a submerged jet with a confinement plate parallel 
te and at a distance, h, from the impingement plate. The 
length of this confinement plate is the same as that of the 
impingement plate. In addition, two spacers are used to fit 
between the impingement and the confinement plates so as to 
form a rectangular channel to which the jet flow would be 
confined. No such confined submerged jet has been studied. 
Actually, only a few investigators (34,45,77,96) who studied 
a laminar impinging submerged jet numerically, considered 
the effect of the presence of a confinement plate. In all 
these studies, the submerged jet was only partially confined 
by a confinement plate parallel to the impingement plate. 
Such a jet, with no spacers between the impingement and the 
confinement plates, is referred to as "semi-confined" 
Submerged jet by the investigators. 

For a semi-confined submerged jet, the flow field is 
very Similar to that shown in Figure 2.1 for an unconfined 
Submerged jet except a recirculation region is induced 
between the impingement plate and the confinement plate 
(45,77,96).. Thisscecirculation region can affect the spread 
of the submerged jet in the free jet region. The inertial 
effects of the upper part of this recirculation derives the 
greater part of its momentum from the jet in the wall jet 
region (96). An interesting result was observed by Van 


Heiningen et. al. (96) in their study of a semi-confined 
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laminar two-dimensional air jet is that the jet contracts 
slightly below the nozzle exit for an initially parabolic 
jet velocity profile. 

The "effects of the velocity profile at the nozzleexit 
on the spread of a submerged jet and the decay of centerline 
velocity in the free jet region were noted by Van Heiningen 
et. al. (96) and Li (45) in their numerical studies of 
semi-confined laminar two-dimensional and axisymmetric 
submerged jets, respectively. A submerged jet issuing with a 
Flat velocity profile spreads and decays faster than that 
with a parabolic velocity profile. For the case of flat 
velocity profile, the momentum just inside the free 
streamline is an order of magnitude greater than that of the 
Surrounding fluid just outside the free streamline. 
Therefore interaction occurs immediately between the high 
momentum fluid at the outer edge of the jet and the still 
fivid Surrounding rt; cauSing a higher rate “of ‘spread of the 
jet. This spreading effect is less important for the 
Submerged jet with a parabolic velocity profile. As 
mentioned above, the jet actually contracts slightly below 


the nozzle exit in this case. 


2.2 HEAT AND MASS TRANSFER DUE TO A SINGLE LAMINAR IMPINGING 


SUMBERGED JET 


Heat and mass transfer characteristics of laminar 
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impinging axisymmetric submerged jets have been studied 
rather extensively. The theoretical and experimental results 
are well correlated. However, heat and mass transfer 
chacteristics of laminar impinging two-dimensional submerged 
jets have recieved lesser attention. In this section, 
laterature related to confined and unconfined submerged jets 
will be reviewed separately with emphasis on the 


two-dimensional case. 


2.2.1 UNCONFINED SUBMERGED JET 

Most of the studies of heat and mass transfer due to 
unconfined laminar impinging submerged jets cover both the 
Stagnation flow and the wall jet regions. But there are also 
a few studies which cover the wall jet region only (79,82). 
In these studies, correlations for heat and mass transfer in 
terms of simple governing parameters of the flow were 
obtained. In general, at a given jet Reynolds number, the 
heat and mass transfer characteristics in the stagnation 
flow region are very much affected by the jet-to-plate 
Spacing, but those in the wall jet region are not so much 


affected. 


2.2.1.1 UNCONFINED AXISYMMETRIC SUBMERGED JET 

Scholtz and Trass (79) studied mass transfer due to a 
laminar impinging axisymmetric submerged jet both 
theoretically and experimentally. Their theoretical 


expression, which has been experimentally verified for mass 
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transfer in the wall jet region is given by 
Cesc, Beet oli? fe2- 31: 


Bere Cc, = cout —-g! (0 )) (25-16) 

and Can OUpe BemtCSertt/3)/(0 (Sc) I (1/3) Coes) 
(-g'(0)) is defined as a dimensionless gradient of 
concentration at the wall. c, 1S a constant and its value is 
dependent on the exterior flux of momentum flux defined as 
bysGlauert, (25). For a flat velocity -profile.at nozzle exit, 
C2=0.426. On the other hand, for a parabolic velocity 
profile, c,=0.458. For comparison with the theoretical 
expression of Equation 2.15, experimental data were obtained 
by measuring the shrinkage of coatings of acetanilide and 
benzoic acid in the range of 1000sRe ,<3000 and 970<Sc<4400. 
Agreement with theory was good. 

Later, Scholtz and Trass (80,81) have given a 
theoretical solution in the stagnation flow region of a 
laminar axisymmetric submerged jet with a parabolic velocity 
profile, at nozzle exit.,The, solution. for inviscid flow in 
the body of impinging jet was first obtained (80) and then 
was used as the boundary condition to solve the viscous 
boundary layer flow near the impingement plate (81). For 
impinging jet with a parabolic velocity profile at nozzle 
exit, the Sherwood number in the stagnation flow region was 


given by 
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f Rey’ = 1.6484 Sc°®-356! 


Sete OOOO Aa) 7 SC Crit es lien Oe tod 


for h/d=0.5, r/d<0.5 and 1<Sc<10. By setting r=0, the 


Stagnation point Sherwood number is given by 


sha = 1.6484 Reg Soe vires 4 C2599 

In a Similar manner, the mass transfer from the 
impingement plate in the stagnation flow region using an 
impinging jet with a flat velocity profile at the nozzle 
exit was calculated by Scholtz and Trass (81). The results 
indicated that, at a given jet Reynolds number, the mass 
transfer at the stagnation point is less than half that 
observed for a jet with a parabolic velocity profile at 
nozzle exit. However, no experimental data were obtained to 
confirm this theoretical finding. 

Mass transfer experimental data were obtained by 
Scholtz and Trass (81) from measurement of the sublimation 
rates of a naphthalene coating exposed to air jet witha 
parabolic velocity profile at nozzle exit (Sc=2.45). For the 
ranges of 5970SRe 41970 and 0 5sh/ds6) ‘good’ ‘agreement was 
reported between the experiment results and the theoretical 
expression of Equation 2.18. In addition, a smooth 
transition from the stagnation flow region to the wall jet 
region was observed in the region of 0.5<r/ds2.25. 


Experimental data in this transition region were correlated 
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in terms of jet Reynolds number as 
Sacre O Re Ce Cr /d yet s C2220) 


for 970SRe,<51970 and Sc=2.45. 

A typical curve showing the variation of local Sherwood 
number on the impingement plate for an unconfined laminar 
axisymmetric submerged a can be plotted from the 
correlations for different regions in Equations 2.15, 2.18 
and 2.20. Since the dependence of local Sherwood number on 
Reynolds number differ in the stagnation flow, transition 
and wall jet regions, the relationship between the local 
Sherwood number in these regions on a plot of Shy versus b/d 


will depend on Re A typical variation of the local 


d° 
Sherwood number on the impingement plate for an unconfined 
laminar axisymmetric submerged jet with a parabolic velocity 
profile at nozzle exit is shown in Figure 2.2 for Re,=1740 
and Sc=2.45. 

Kapur and Macleod (36,37,38) have applied the 
techniques of holographic interferometry to the 
profilometric measurement of mass transfer rates at 
solid-fluid surface exposed to an unconfined laminar 
impinging axisymmetric air jet. Experimental results for the 
wall jet region in the range of 255<Re,<1870 were in 
excellent agreement with the theoretical predictions of 


Scholtz and Trass (79). Mass transfer coefficients near the 


Stagnation point were also reported for Re y= 1340 and 
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FIGURE 2.2 : VARIATION OF LOCAL SHERWOOD NUMBER ON 
IMPINGEMENT PLATE FOR AN UNCONFINED 
LAMINAR AXISYMMETRIC SUBMERGED JET 
WITH INITIAL PARABOLIC PROFILE 
(Reg=1740, Sc=2.45) 
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h/d=0.5d, id and 1.5d. Within the narrow range of the 
jet-to-plate spacings used, the mass transfer coefficients 
in the stagnation flow region were found by Kapur and 
Macleod to be independent of jet-to-plate spacing. 

Masliyah and Nguyen carried out a number of 
investigations on mass transfer due to laminar axisymmetric 
M50;,51), “square ©(50;5'1) » w~ectanquilar °(50,52) and 
two-dimensional (53) air dt by uSing a laser holography 
interferometric technique. In their study of axisymmetric 
jet, experimental results for the wall jet region of Re 


d 


=1145 and 1420 were found to be dependent on Re,° °° and 
Ged io 3) «Thetexponent s ‘of Re, and (r/d) were in good 
agreement with the theoretical findings by Scholtz and Trass 


(79). 


2.2.1.2 UNCONFINED TWO-DIMENSIONAL SUBMERGED JET 

Schwarz and Caswell (82) studied the heat transfer 
Characteristics of a laminar two-dimensional wall jet by 
solving the flow equations analytically for a wall jet as 
developed by Glauert (25). Their eee a ien for mass 


transfer in the wall jet region is given by 
Shas Came, Ga oe pen a C2e2.\s) 
Wheres ¢, = c, (-g",(0))) (2522) 


c, iS a constant whose value is dependent on the exterior 


flux of momentum flux as defined by Glauert (25). For a flat 
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velocity profile at nozzle exit, c,=0.446. On the other 
hand, for a parabolic velocity profile, c,=0.465. 

Gardon and Akfirat (19) studied experimentally the heat 
transfer characteristics of unconfined impinging 
two-dimensional air jets for both laminar and turbulent 
cases. Although the nozzle exit velocity profile of the jet 
was not specified, due to the short nozzles used by Gardon 
and Akfirat, the nozzle paae velocity profile was most 
likely flat. In their studies of a single unconfined laminar 
two-dimensional jet (Re, =450, 650 and 950), their main 
interest was the effects of jet Reynolds number and 
jet-to-plate spacing on the stagnation point Nusselt number. 
For jet-to-plate spacing less than the length of the 
potential core of the submerged jet (0.5<h/b<5), Nusselt 
numbers at the Ee ney point depend on the jet Reynolds 
number only. They are approximately proportional to Re,°-°*?. 
At large jet-to-plate spacings, stagnation point Nusselt 
numbers’ varied» with (h/b)7°:* for Re, =450. However, when 
Reynolds number was increased to 650 and 950, this monotonic 
decreasing behavior of stagnation point Nusselt numbers with 
increasing jet-to-plate spacing does not exist. Instead, a 
maximum of stagnation point Nusselt numbers occurs when 
jet-to-plate spacing is in the range of 10<h/b<15. The 
behavior was explained by Gardon and Akfirat (18) and it is 
caused by the penetration of mixing induced turbulence to 
the centerline of initially laminar jets. The interactions 


between the centerline turbulence and the approaching 
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centerline velocity in the stagnation flow region is the 
reason why this non-monotonic behavior of the stagnation 
point Nusselt numbers was observed. It is because the 
centerline turbulence increases with increasing jet-to-plate 
Spacing and the approaching centerline velocity decreases 
with increasing jet-to-plate spacing due to the decay of 
centerline velocity. Gardon and Akfirat (19) concluded that 
for small jet-to-plate eee ice it is likely that the 
potential core velocity will be independent of the 
jet-to-plate spacing, and the stagnation point Nusselt 
numbers should show a corresponding trend. Once the core has 
been engulfed by the mixing region, the decay of the 
centerline velocity would tend to decrease the heat transfer 
rate at the stagnation point, but centerline turbulence 
would have an Bene eens ce ett At lower Reynolds numbers the 
former should be more important, whereas at the higher 
Reynolds numbers the influence from centerline turbulence 
Should become more important instead. 

Miyazaki and Silberman (58) have analysed theoretically 
the flow friction and heat transfer characteristics on a 
heated or cooled flat plate with unconfined laminar 
impinging two-dimensional jet with initial flat profile. A 
potential flow solution was obtained to provide the 
distribution of the main-stream velocity. The boundary layer 
and energy equations were then solved numerically by a 
finite difference method. The correlation obtained by 


Miyazaki and Silberman in terms of mass transfer variables 
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at the stagnation point is given by 
Shy SIE OUO MRE. Sow. 8 ts U2 323) 


for h/b21.5 and 0.7<Sc<10. Variations of local Nusselt 
numbers on the impingement plate were also evaluated, though 
no correlations were given. For h/b21.5, Nusselt number 
exhibits a maximum at ojcege erie point and remains 
constant within the stagnation flow region. Away from the 
Stagnation flow region, Nusselt number decreases 
monotonically and finaily tends to vary proportionally with 
x~°-*5 as for the flow over a flat plate. 

Sparrow and Wong (91) studied the impingement mass 
transfer due to an unconfined laminar two-dimensional 
Submerged jet with initial parabolic velocity profile. By 
using the naphthalene sublimation technique (Sc=2.5), 
variation of the local mass transfer coefficient on the 
impingement plate was determined for Re,=150, 300, 450, 650 
epar7o0, and h=2b, 5b... t0be 15b and 20b. In general, the 
local mass transfer coefficient was found to decrease 
monotonically with increasing distance from the stagnation 
point, but correlations were not given. Sparrow and Wong 
also studied the effect of jet-to-plate spacing on the 
Stagnation point Sherwood number. The non-monotonic behavior 
of stagnation point transfer coefficients with increasing 
jet-to-plate spacing mentioned by Gardon and Akfirat (19) 


for higher Reynolds numbers was again observed. Stagnation 
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point Sherwood numbers measured by Sparrow and Wong (91) 
were converted into Nusselt numbers by using the analogy 
between heat and mass transfer, Nu, = (Pr /Scu <4 Shp, in 
order to compare with the stagnation point Nusselt numbers 
measured by Gardon and Akfirat (19). The results of Sparrow 
and Wong are about 30% higher than those of Gardon and 
Akfirat. The main cause of this different is the nature of 
the initial velocity profiles of these two studies. Similar 
observation is made by Scholtz and Trass (81) for an 
unconfined impinging axisymmetric submerged jet that the 
Mass transfer at the stagnation point for a jet with initial 
flat velocity profile is less than that for a jet with 
initial parabolic velocity profile. 

Sparrow and Lee (90) analysed fluid flow and heat and 
mass transfer characteristics associated with the 
impingement of an unconfined laminar two-dimensional 
Submerged jet with an initial parabolic velocity profile. 
The velocity field within the impinging jet was solved 
within the framework of an inviscid flow model and the 
results were used as input for the analysis of the boundary 
layer heat or mass transfer on the impingement plate. The 
transfer coefficients were found insensitive to the 
jet-to-plate spacing within the range investigated, which 
was 0.375Sh/bs1.5. The dependence of stagnation point 


Sherwood number on Reynolds number was given as 


Shp HON AD oe (2.24) 
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for Sc=2.5. Fair agreement between Equation 2.24 and the 
experimental data of Sparrow and Wong (91) for h/b=2 was 
obtained. 

In the experimental study of an unconfined laminar 
impinging two-dimensional air jet with initial parabolic 
velocity profile by Masliyah and Nguyen (53), a regression 


equation for the local Sherwood number was given as 


SU mea Ooo gke Fs ley (2225) 
Pom eoC=2.05, 90<Re, <300 and 1£$x/bs30. All the experimental 
data reported were taken with jet-to-plate spacing equal to 
4b. They also pointed out that the experimental results with 
a jet-to-plate Spacing equal to 8b gave similar results. An 
overall general agreement between Equation 2.25 and the 
experimental results from Sparrow and Wong (91) was 
obtained. In addition, the exponent, -0.73, of (x/b) in 
Equation 2.25 is in fair agreement with the theoretical 
analysis of Schwarz and Caswell (82) for a two-dimensional 
wall jet, which is -0.75. However, their exponents of 
Reynolds number do not agree. 

It would be appropriate to conclude that there are no 
unified correlations of the variation of local Sherwood 
number for unconfined two-dimensional submerged jet. In 
general a typical distribution of local Sherwood number on 
impingement plate for the two-dimensional case is similar to 


that for the axisymmetric case. A plot is shown in Figure 
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2.3 for Re,=400 and Sc=2.85 by using the correlations of 
Equations 2.21 2.24 and 2.25 and assuming that the local 
Sherwood number in the stagnation flow region remains 


constant. 


2.2.2 CONFINED SUBMERGED JET 

Heat and mass transfer due to a confined laminar 
impinging submerged jet have been studied numerically. All 
Studies cover both the stagnation flow and the wall jet 


regions. 


2.2.2.1 CONFINED AXISYMMETRIC SUBMERGED JET 

Saad et. al. (77) have simulated the flow and heat 
transfer characteristics of a semi-confined laminar 
impinging axisymmetric air jet. The vorticity-stream 
function formulation of the Navier-Stokes and the energy 
equations were solved numerically. Effects of uniform 
Suction and nozzle exit velocity profile on the flow and 
heat transfer characteristics were studied. Saad et. al. 
also observed the finding of Scholtz and Trass (81) 
concerning the effect of the initial velocity profile even 
with the presence of a confinement plate. By comparing their 
results with the experimental data obtained by Scholtz and 
Trass (81) who used an unconfined axisymmetric jet, Saad 
et. al. showed that the presence of a confinement plate has 
only a minor influence on the stagnation point heat transfer 


for the range of 2<h/ds<4 and 450<Re,<2500. 
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30 


ro 
7 =i ; 
i ke er ae 
: - ‘jee : Tee 
, site - ; ae # - ip sone Re is 
7 Pag " 
HAs 7 Ne i» . on: way ho oe 
Me te ut Pad 
i 7 * 
: ad Snr eae, gaa * 
’ 
t a a 1 a 


ge ay S at ey |e DA 
. ? 


ee | ee 


* 
bs 
72 


. 
eet sina 
| 
. 


! 
* 
™ 
x 
5 
4 


~~ . © em 
\ 
, 


Cd 


oo, 


3 


Other numerical studies of semi-confined laminar 
axisymmetric jets included the studies by Huang et. al. (34) 
and Li (45). Huang et. al. (34) have numerically solved the 
Navier-Stokes equations in their primitive form with the 
energy equation for a semi-confined laminar axisymmetric 
Swirling jet. Li (45) has studied numerically the 
Simultaneous heat and mass transfer under a semi-confined 


laminar impinging axisymmetric jet. 


2.2.2.2 CONFINED TWO-DIMENSIONAL SUBMERGED JET 

Only one study of heat transfer due to a laminar 
impinging two-dimensional submerged jet by Van Heiningen et. 
al. (96), considered the effect of the presence of a 
confinement plate. By using the same numerical technique as 
Saad et. al.'s (77), Van Heiningen studied the effects of 
uniform suction and nozzle exit velocity profile on the flow 
and heat transfer characteristics. They found that for a 
jet-to-plate spacing equal to 4b and 100<Re, <950, the 
stagnation point Nusselt number for an initial parabolic 
velocity profile was between 1.5 and 2 times the value for 
an initial flat velocity profile. For a flat velocity 
peotitve, Simider behaviors of Skin friction. factors and 
Nusselt numbers along the impingement plate mentioned by 
Miyazaki and Silberman (58) were also observed, although the 
results obtained by Miyazaki and Silberman were consistently 
higher than those obtained by Van Heiningen et. als. nus 


difference as pointed out by Van Heiningen et. al. is due to 
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the assumption of potential flow for the flow outside the 
boundary layer used by Miyazaki and Silberman. For the case 
of a parabolic velocity profile, Sherwood numbers measured 
by Sparrow and Wong (91) were converted into Nusselt numbers 
by using the analogy between heat and mass transfer as 
mentioned above in order to compare with the numerical 
results from Van Heiningen et. al.. Good agreement was 
obtained both in the stagnation flow and in the wall jet 
regions indicating that the presence of a confinement plate 
has only a minor effect on the heat and mass transfer 
Characteristics. It is worthwhile to mention that the effect 
of jet-to-plate spacing was not studied by Van Heiningen et. 
al. as only one jet-to-plate spacing of 4b was used in their 


study. 
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3. BASIC EXPERIMENTAL TECHNIQUES 


Measurement techniques used to study problems in heat 
and mass transfer are numerous and often quite complex. Mass 
transfer techniques have usually been used to avoid the 
errors inherent in heat transfer measurements; for example, 
conduction errors ina ferecd convection heat transfer 
System even with near-ideal design often cannot be reduced 
to a sufficiently low level due to the finite conductivity 
of even the best of insulators. Mass transfer techniques, if 
properly utilized, do not have an analogous conduction error 
and for this reason, they are very valuable when precise 
local measurements are required. 

A novel profilometric technique for determining local 
mass transfer coefficients at solid surface was first 
described by Macleod and Todd (47). Instead of a volatile or 
soluble surface coating (i.e., naphthalene coating) 
customarily employed in profilometric work, the coating used 
is composed of non-volatile polymer capable of absorbing 
volatile or soluble swelling agents. Rates of transfer of 
swelling agent to or from the surface by the experimental 
fluid stream can then be evaluated from measurements of the 
swelling or shrinking rate of the coating. In most practical 
cases, the polymer coating swelling or shrinkage is of order 
of 10°-*‘m and consequently mechanical means for these 


measurements are difficult and unreliable. Holographic 
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interferometry in the conjunction with the swollen polymer 
technique of Macleod and Todd (47) proposed by Kapur and 
Macleod (36,37,38) and Masliyah and Nguyen (50,51,52,53) was 
found to be a powerful and convenient means for measuring 
such a small change in polymer coating thickness. The unique 
advantage of this optical method is that it does not 
influence the process examined. In. addition, instead of 
Dornteby jpoint WS acueen hee information about a whole field 
of interest can be obtained by the evaluation of 


photographs. 


3.1 SWOLLEN POLYMER TECHNIQUE 


The essence of the swollen polymer technique is to coat 
the flat surface under study with a thin layer of an 
non-volatile silicon polymer and to swell this coating ina 
bath of ethylsalicylate. Transfer of ethylsalicylate from 
the polymer coating when the coating is subjected to an air 
jet results in local changes in the degree of swelling. 
These changes cause the polymer coating to shrink (47). It 
has been shown that the local shrinkage of the polymer 
coating over a known time within the "constant rate period” 
provides an accurate measure of the local mass transfer 
coefficients over the surface (47,51). In other words, in 
order for the swollen polymer technique to be useful in mass 


transfer studies, it is essential that a "constant rate 
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period” is maintained while measurements of mass transfer 

are taken. 

Macleod and Todd (47) have pointed out in their study 
of this technique that the following important conditions 
hold within the constant rate period: 

1. Swollen polymers having virtually zero heats of mixing, 
exhibited no total volume change on swelling, i.e., the 
swollen polymer volume is the sum of the volumes of the 
Gry polymer and that of the pure swelling agent. Thus, 
the shrinkage of the polymer coating is proportional to 
the amount of swelling agent transferred from the 
coating. 

2. Although the vapor pressure of the swelling agent over 
the surface depends both on temperature and on the 
composition of a Caen and must change as 
volatilization proceeds at constant temperature, it is 
reasonable to assume that a change in the partial vapor 
pressure of up to 5% at the polymer coating during an 
experiment has a negligible effect on the determination 
of mass transfer coefficients compared with other 
sources of error. The»so-called."constant. rate period” 
is defined as the time period when the vapor pressure 
remains within 5% of its initial value during 
volatilization. Within this period, the partial vapor 
pressure of the swelling agent at the surface can be 
assumed to be effectively the same as the vapor pressure 


of the pure swelling agent. Such a constant rate period 
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was found experimentally and theoretically to be as high 
as 2500 s. 

The overall resistance to mass transfer is predominantly 
in the gas phase. 

The effect of diffusion of the swelling agent parallel 
to the surface from regions of low mass transfer to 
regions of high mass transfer can be neglected. The 
behavior of the coating at any given point is therefore 
governed by the local mass transfer in the direction 


normal to the polymer coating only. 


Because of all these conditions mentioned above the swollen 


technique is a powerful technique in the profilometric 


Secermination of local mass transfer coefficients. 


The practical advantages of the swollen polymer 


technique are summarized as follows: 


1. 


The coating itself can be used repeatedly; for each new 
experiment it can be re-activated by re-immersion in the 
swelling agent bath. 

The Schmidt number of the system can be varied; this can 
be accomplished by a changing the swelling agent alone, 
without changing the experimental fluid or the polymer 
coating. 

Because the mass transfer rate is independent of the 
coating thickness, no uniformity of thickness other than 
that implied by the need for hydraulic smoothness at the 
polymer surface of the coating is required. 


Unlike the surface of a subliming solid, the surface of 
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the shrinking polymer maintains a nearly constant 
optical quality. It is therefore possible to record its 
shrinkage by fairly standard techniques of holographic 
interferometry. 
It 1S propriate to mention that the swollen polymer 
technique is not suitable for determining local mass 
transfer coefficients for a curved surface and also when the 


local mass transfer rate is high. 


3.2 LASER HOLOGRAPHY INTERFEROMETRIC TECHNIQUE 


Many optical methods have been used in heat and mass 
transfer studies; reviews of this methods have been 
presented by Hauf and Grigull (29) and Goldstein (26). Up to 
a decade ago, the most common methods used were Mach-Zehnder 
and Michelson interferometry. In 1948, Gabor invented a new 
optical recording technique which he called "image formation 
by reconstructed wavefronts" at that time. Not until two 
years later, the word "holography" was apparently first 
applied by Roger (43) to describe this new photographic 
procedure. In contrast to ordinary photography, by which 
only the amplitude of the reflected light from the object is 
recorded, holography allows the recording and reconstruction 
not only ene amplitude but also the phase distribution of 
wavefronts. As holography demands a highly coherent light 


source, there is a technological break-through after the 
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invention of the laser about twenty years ago; since then, 
the holographic techniques have played an important role in 
the studies of heat and mass transfer (54). 

The general theory of holography is so comprehensive 
that for a detailed description the reader is referred to 
Bue ebiterature -(13743744,73). In this:workyonly. the 
principles necessary for understanding the holographic 
measurement techniques Sie mentioned. 

To record a hologram of an object, a beam from a 
coherent monochromatic light source (usually a laser) is 
Split into two separate beams. One beam, so called the 
"object beam or wave", illuminates the object, and the 
reflected scattered light falls directly onto a photographic 
plate. The other beam, so called "reference beam or wave", 
travels by a separate path bypassing the object, and falls 
on the same area of the photographic plate. If both these 
two waves are mutually coherent, they will interact and form 
a Stable interference pattern. This interference pattern can 
therefore map an exposure pattern across the surface of the 
photosensitive emulsion on the photographic plate. When the 
photographic plate is developed and fixed, the plate with a 
pattern of dark fringes is called "hologram". The fringe 
pattern which in general consists of about 1000 lines per 
millimeter, contains all the visual information about the 
object. The amplitude of the light from the object is 
recorded in the form of different contrast of the fringes 


since the intensity of light incident on the emulsion will 
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be the square of the sum of the amplitudes from the 
reference and the object waves. The phase of the light from 
the object is recorded in the spatial variation of the 
fringe pattern since when the pathlenagths between the 
reference and the object waves differ by one-half 
wavelength, interference will be destructive and no light 
energy will be available to expose the emulsion. 

When the hologram Pan eed by the reference wave, 
the fringe pattern acts like a diffraction grating. A number 
of light waves are generated by the interaction of the light 
with the grating. The "zero order" wave propagates in the 
Same direction as the incident light. Beyond the two first 
order waves, second, third and higher order waves also 
occur. One of the first order waves travels in the same 
direction as the original object wave and it is responsible 
for the virtual image of the original object. The other 
first order wave is responsible for the real image of the 
object which usually appears unsharp and highly distorted. 
Since the scattered light from one point of the object is 
recorded in each part of the hologram, this point can be 
Seen even through only a small fraction of the plate. It can 
therefore be observed from various angles, limited only by 
the size of the hologram. Thus, in the reconstruction stage, 
a truly three-dimensional picture of the image is obtained. 

Making use of these recording properties, several 
exposures can be made on the same photographic plate. 


Illuminating this multiple exposed hologram with the 


eh 


3 


leet sag moe, aebug i laps fet Se ete ht eS 
mex? aieed, ag? $8. paadaredt +Bavey ate 


ad to nepeises fs oLinge ane ap besiesies ee 22 


at? asewdod en {He shadagyi cree aps 19038 aa? 


j gd SBS 1B eayaw: Palanied: onic 


3g02l: on ond. sv bday sdRsi ond hin cnuretnan 


i 
} - 
> 
ri, 
s 
7 

ia? 


Ae 
FI 


HOLS hua a69. Seeks ie, 6 del sere. od iim 


- 


SiS1-su9 Ja: So aee iether ii Pee LOL OS edt 


(i, — 
smn) A sae SHAD Pees ee: ac0e caeieg ame 


pili sai 26 necaestetal ens ve wala teiee ets 22a 


12927 owe of 2 SnoysA AGB +) SRR . aie os: noi soe 


o8he SVs", 2eht0.1sigin eas ph. gone » RON RES 
4] - 


amee Seice it  ehewsit ceyvsr sabre dazed. sae 30- ee 
nod i 
os {ceaaggyeea: et ..71 Se over sSeede, Ser toes ere 
mic a? ,tgera¢ [serieew was 26 epens) Levioary 


of] 2 wiswt | tse% $a3' 393 eletdeuGaes ve eva, 19f%048 
Serictase ¢higéd Gre Hsétenu) 6teeaas idl tmawies. dstibaea 


és 393060 60 20 20Le¢ ano meet? jfoi : ote tteae od7 


u A. 


$G 163. HOG Sine yRaueolod one to z2an dsad at, 

isa 45° ,e3 eta seg to notes? lisva 2 ylno dgeentg nerf 
tt ving Setinit 2Sipas evolisv nes? Seviseds ed e107 

, Seta: ni sau azaodes atts Ai ,29dt .neipolod. edz. ao 0 et ? 
‘Beatsida a2, Septl ae jc stussiq Lanotensaib-seeds eh a ie 
tasever 4afregeta Qrabiiesy. se644.%¢ sav a 5 i. 
-9tela oidganpevatic auaz ang nD SREM ot meDs 7 

- be site apaled: arnt wena aap mre 


: 
> 


40 


reference wave after processing, the object waves due to 
different exposures will all be reconstructed and if they 
differ only slightly from one another they can interfere 
with each other to form fringes. This is the basic idea of 
holographic interferometry. 

Normally holographic interferometry is made using a 
double exposure. The interference fringe pattern from the 
two reconstructed object waves after two different exposures 
can be photographed and therefore can be analysed at any 


time. 


3.2.1 DOUBLE EXPOSURE HOLOGRAPHIC INTERFEROMETRY 

In double exposure holographic interferometry, a 
hologram of the object, which in this work is the polymer 
coating, is made before subjecting the polymer coating to 
the air jet. After the polymer coating is subjected to the 
air jet, a second exposure is made onto the same 
photographic plate. After processing the photographic plate, 
the double exposed hologram is repositioned and illuminated 
by the reference wave. Now both object waves, the first one 
is due to the original undisturbed polymer coating and the 
Second is due to the shrunk coating, are reconstructed. The 
interference from these two object waves forms fringes 
indicating the change in thickness of the polymer coating 
between the two exposures. These interference fringes, so 
called "frozen fringes", are contours of equal coating 


shrinkage or equal mass transfer and can be photographed 
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with a 35mm camera. 

The main difference of this technique from classical 
interferometry, such as Mach-Zehnder and Michelson 
interferometry, is that the object wave is compared to 
itself. Since both waves pass through the same optical 
set-up, any imperfections of mirrors and lenses are 
eliminated. Another advantage of this technique is that it 
can provide a complete record of the local mass transfer 
coefficients at all points on the surface. No other 
technique of heat or mass transfer measurement hitherto 
reported has this capability. Therefore, mass transfer data 
can be provided both in precision and completeness. 

However, there are also some drawbacks. The 
photographic process of this technique, particularly at the 
initial stages, requires a high degree of specialized skill. 
Its success is greatly dependent on the correct choice of 
several optical and photographic processing variables, such 
as ratio of beam intensities, length of exposure and 
conditions of hologram development. More importantly, in 
this technique rt is difficult “to “identify ‘the “order “ot “any 
Particular fringe. The "frozen tringe” ‘pattern’on ‘the *double 
exposed hologram only provides information about differences 
between the shrinkage of the coating at different points. 
The actual value of the shrinkage can only be found if the 
order of fringe at a given point is known. This shortcoming 


can be overcome by real time holographic interferometry. 


o sit eee oe ae 


(aa teemts: nanarandgt Ae phe? toe eels — 
| dia ara Gite ~ebiipae seamed ion ie 


in a ‘ 


a2 Saga ar Sew Thane wt Aber s ib oe 
fost igo ome ‘aes aes 2 8 esveu ated 
28 Saens it vee Sow te) "aoe! P2873 wind 


nis et Suet oe GEE ‘ep Haelteg | 938 708K me 
isiengt? ceet Sets — ens” te aivte dy. Ingod (6 
a) wy ™ "i e aaah! i 
tH4v0 OW eos iyor SS Ae siting ile 26 o28 
le a 4 5 " ~ . < as ah —— a 
stamieiet Yneeetsakin: deSengrd BRT. " teed “77 
45e5 swataetsts |etee stetart . git Pitangs aids ba 
sons te lye Soe Soir i ‘ayied' Bot is 
wT Ledoedwash aise! bale sia‘ oieit? 
t ff r on "\ z 7 tiene % : se to I Bas . 
wis te yh al vo sea . Supoaye Ss Pas %6 mas aot ping 


‘Bi 


bitde peed ‘his baie a: ho tenga. eri ‘ mata toni 
io eolodds Sasiiteo Say aa yritrcaigeds’ piven ae 22, oes 
dogs ealge- "av pales Bile pony aa nek: 7 

as Sue bans % Wieast .adisiive fal mad’ we i 

nt , vi oneindeiae: Say. TseQsiavss> cancoled io | 

Yin tb subse ats BPs 43 Halstexig 22 4: on 

eldund siz ne: ee ses" ety mores 


a26 i tvO70 elas congbtelt 


oe f : : ro - 
Srsincan sis! ecuiaalal 
7 >" = " 


ve 


£2 


3.2.2 REAL TIME HOLOGRAPHIC INTERFEROMETRY 

Real time holographic interferometry can be done by 
using a precise and adjustable hologram mount. After the 
first exposure, the hologram is processed and replaced in 
its Original position. The "reference image", which is 
actually the image of the object before the mass transfer 
experiment, is then reconstructed continuously by 
illuminating the hologram with the reference wave. Then the 
mass transfer experiment is started, the reconstructed light 
can be superimposed onto the changing object wave. The 
changes of the interference pattern can then be continuously 
observed and the order of any particular fringe at any given 
time can be easily determined by counting fringes at a given 
location as they appear in real time. 

Another apparent advantage of this technique is that it 
is possible to distinguish between a monotonic anda "U" 
type coating shrinkage between which the double exposure 
technique fails to distinguish. In double exposure 
technique, the appearance of "frozen fringes" is simply due 
to change in depth regardless of whether such a change is a 
Shrinkage or a rise. But in real time technique, the ability 
to determine the order of each fringe gives a complete 
information on how the coating shrinkage proceeds. Therefore 
a "U" type coating shrinkage can be easily observed and 
recorded. 

The optical set-up used in this work which is suitable 


for both double exposure and real time holographic 
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interferometry will be discussed in Chapter 4. 
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4. EXPERIMENTAL SET-UP AND PROCEDURE 


The experimental set-up in this work is shown in Figure 
4.1 which includes the optical set-up for recording a 
hologram and also the experimental apparatus for the mass 
transfer experiment due to an impinging submerged jet. With 
the hypothenuse surface a a right-angle glass prism as the 
mass transferring surface, where the polymer coating is 
applied, it becomes possible to have the entire optical 
set-up to be located behind the mass transferring surface. 
Since both the prism and the polymer coating are 
transparent, the mass transferring surface can be observed 
from the optical system all the time during the progress of 
the experiment. Therefore, such an arrangement is suitable 
not only for double exposure holographic interferometry but 


also for real time holographic interferometry. 


4.1 OPTICAL SET-UP 


A concrete table with a concrete slab of size 2.74 x 1.37 x 
O-13amis used as the working surface for the optical set-up. 
The function of heavy slab is to prevent floor vibration 
from being transmitted to the optical components. 

The optical components include: 


1. a Spectra Physics model 125A 50mW He-Ne laser which 1s 
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FIGURE 4.1 : ARRANGEMENT FOR THE EXPERIMENTAL SET-UP 
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used as the source of the monochromatic coherent light 
oftwavetengthfequal.to 632.8 x 10-*m, 

a Prontor-Press shutter which is used to set the 
exposure time for making a hologram, 

an Elomag beam splitter-attenuator, model VBA-200, which 
is a variable reflectivity aluminum mirror for use in 
splitting the laser beam into the reference and the 
object beams and in adjusting the ratio of the beam 
intensities, 

a Spectra Physics model 332 spatial filter and expanding 
lens assembly (expanding lens L3, aperture A3) mounted 
on a Spectra Physics model 306A precision optical mount 
with base and vertical post which is used to eliminate 
Spatial noise and produce a smooth intensity profile 
across the eeoadoee object beam, 

two Spectra Physics model 386-11 utility mirror mounts 
(includes model 564 mirrors) which are used to reflect 
the reference beam in order that the reference beam can 
travel bypassing the object, 

a Spectra Physics model 332 spatial filter and expanding 
lens assembly (expanding lens L4, aperture A4) mounted 
on a Spectra Physics model 306A precision optical mount 
with» base and vertical post. which is used to eliminate 
spatial noise and produce a smooth intensity profile 
across the expanded reference beam, 

a piece of 0.3 x 0.2m frosted glass mounted on a heavy 


metal base, of which its reconstructed image from the 
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hologram is used as a white background for the 
interference fringes, 

8. a photographic plate holder which is designed to hold a 
half size Agfa-Gevaert 10E75 AH (0.1016 x 0.0635m) 
photographic plate for double exposure holographic 
technique, or an Elomag immersion-type X-Y 
micropositionable photographic plate holder, model 
MPH-45W, including PC-45 plate carriage which is used to 
hold and to reposition a full size Agfa-Gevaert 10E75 
(OS1016 x" O-i27m)) phetographic® plate’ for real time 
holographic technique. Photographs of these holographic 
plate holders are shown in Plate 4.1. 

The positions of the optical components and the polymer 
coating surface are determined according to the Holo-Diagram 
technique (1,2,3,4). The polymer coating is placed along the 
ellipse whose foci are the beam splitter and the center of 
the photoghaphic plate and away from the line joining the 
two foci. In this work, the distance between the beam 
splitter and the center of the photographic plate is fixed 
at approximately 1.2m. The length of the object surface, 
which is the hypothenuse of a right-angle glass prism is 
0.125m. Because of the geometry of the optical set-up in 
this work, the angle of illumination and observation of the 
object surface is restricted to 45° by use of the 
right-angle prism. Abramson (1) has mentioned that in 
holographic interferometry, if the object is illuminated and 


looked upon in a direction that makes an angle B, to the 
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(b) REAL TIME HOLOGRAPHY 


PHOTOGRAPHS OF HOLOGRAPHIC PLATE HOLDERS 
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object surface, the distance between two dark fringes will 
correspond to a difference in thickness of half the 
wavelength of the light used multipled by a constant k,. 
Abramson defined the constant k, as the reciprocal of cos 
Boek in, this work,.thesvalue of k, can be easily evaluated as 
meovenoretherreciprocalsof.cos) 45°. According) to the 
Holo-Diagram, the object in this work must be placed along 
the ellipse passing through the point at which the k, value 
is 1.414 and the distances to both the two foci are equal. 
Due to the finite length of the object surface, the light 
that is reflected from its different parts must travel 
different distances. In other words, the k, values along the 
object surface at different parts are always different. 
Consequently, the pattern of fringes moves relatively to the 
object when viewed from different directions (2). In this 
work, the length of the object is only 0.125m compared to 
the distance which is 1.2m between the two foci, therefore 
the movement of the fringe pattern is not serious. In 
addition, this shortcoming can be overcome by looking at and 
photographing the fringe pattern at a fixed direction all 
the time. This can be accomplished by mounting an aluminum 
plate with a 0.013m diameter hole at the center on the 
photographic plate holder (see Plate 4.1). The fringe 
pattern is then viewed through the circular hole. 

The positions of the two mirrors are adjusted so that 
the path lengths of the object and the reference beams, 


measured from the beam splitter to the center of the 
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photographic plate have the same value. Consequently, the 
maximum coherent length of the laser is utilized so as to 
give the maximum object field of depth. 

An overall view of the optical set-up is shown in the 
phOotographs—in Plate 4.2. 

To mee oae a hologram of the polymer coating from this 
optical set-up, the laser beam with diameter about 0.002m is 
split into a reference and a object beams. The ratio of beam 
intensitives can be altered at any time by rotating the 
mirror of the beam splitter-attenuator to different surface 
reflectivities. The reference beam is directed by two 
mirrors bypassing the polymer coating and is expanded by the 
Spatial filter and expanding lens assembly. The expanded 
reference beam falls on the photographic plate. The object 
beam 1S expanded by the spatial filter and expanding lens 
assembly so that it is wide enough to cover the whole 
object. The object beam travels through the prism and the 
polymer coating and it experiences a total internal 
reflection at the surface of the coating. Finally, the 
object beam leaves the prism and falls on the photographic 
plate. The interference pattern generated by over-lapping 
object and reference beams is recorded to form the hologram. 

In both the double exposure holographic technique and 
the real time holographic technique, on reconstruction of 
the hologram, the object images recorded before and after 
the change in coating thickness interfere with each other to 


form fringes. The absolute magnitude of the change of 
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PLATE 4.2 


OVERALL VIEW OF THE OPTICAL SET-UP 
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coating thickness is given by 
r' =n BA | (4.1) 


where r' is the displacement of a given object point from a 
fixed point, » is the wavelength of light, and n is the 
Eringe order of the point under consideration. B is, in 
General, a function of arial position of the object point 
relative to the beam splitter and the photographic plate. In 
this work, all the positions of the optical components 
remained the same, therefore B is considered a constant. By 
Simple geometrical manipulations, B can be easily evaluated 
from the optical set-up by calculating the difference in the 
path lengths of a iight ray before and after a change of 
coating thickness. Figure 4.2 shows the paths of a light ray 
before and after a mass transfer experiment for a given 
change of coating thickness, r'. The change in the light 
path length between a-a and b-b is the same between c-c and 
d-d, and is given by 
Zh ie 5 Sine 6, 0SiNn=67)) f COS ps (a2) 
where n andn_ are refractive index of the swollen polymer 
and the glass prism respectively (n,=1.428 and nj=1.5). The 
angle £, as shown in Figure 4.2 is equal to 45°. The angle 
8B; = sin-' (np sin £, / ng) is calculated as 47°58" .0AS the 


appearance of a fringe is due to a change of light path 


f a 

6 me i gaieg toaisy mig ye sdasiabieli sal 
effi 2t a nab seid io dtpaaipvey ite ae 

| gmt @ pneige? ogiendh, ‘yehee anisq ado pope 

mag 3 | i sit 7 pe aI ba 3) or to. a 

sewrameicneteg “ie bine | aint ee oad 563 oe 

esoeneqre> Laataqe os Ve RAossiaeq per a 


vi .dAmgeicdioo «& pesebi ued 2) ff. etetezets men MOS 


‘| 
3 
‘ 


1 oy ; 


»S364 


re 


me 
4 


ssisyisvs yliges et ase @ , act yelogenas Leots sets 
ots ni suneweltib edd. palyaleealas yd gu-see two 
& a yen ) & te32e Gee acoted at ‘pit so te 
ot aipdl @ to kWdag wity ude See. amyl asoeil 
perce & 308 trem? regis rvingess ater 's se2dn 4 


a 
i‘ 


fei) a2 2, dpeets sft). *S .eeendahds ane vod 


5 
| 
th 


o-o peawter amee ete ef d-d Sta.e-s nese red 


ot ve mage 


os 7 &,a08 \-4,8 aia... ota aft ne, 


\ s 
jeezlcg aollowe df So asbei svisoesiey ofe " 


qT ./ 2. tage Bos $88. myn wevtsooqea maton Ws 


aisns ant "22 29 teuge at 5,0 — ® 
-. ) oa an twats * 2s bataluoles htt 5 | Gahan 


e9s3 a 


_ aaa’ rapt te, aeneda a 98. 


My, 
at 


a 
fy. 
ae 7 


S36) 


Polymer Coating 


B2 
x. 
- Glass Prism 
By 
a is b 
\ ne Ue 
se 7 
' 7 
illuminating no viewing 
beam beam 


———> before mass transfer experiment 


— -~ after mass tranfer experiment 


FIGURE 4.2 : LIGHT PATH LENGTH BEFORE AND AFTER 
MASS TRANSFER EXPERIMENT 
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lengths of n\ /2, then the change of coating thickness r' 


for nth fringe is given by 
Pee De COSs) Bao 74 (n 2 = Np Sone isin) Bs) C453) 
Comparing Equations 4.1 and 4.3, the constant B is given by 
Baa sCOSUNE 4 / (44n. - np Mn, Sine ps (4.4) 


Hence, by knowing the incident angle £, and the refractive 
angle B., the value of B can be calculated by using Equation 
4.4. If the value of B is known and the fringe order n is 
determined, the absolute magnitude of the polymer coating 


Shrinkage at a given point can be evaluated. 


4.2 MASS TRANSFER EXPERIMENTAL SET-UP 


The mass transfer experimental set-up is shown in the 

photographs in Plate 4.3 and it includes the following: an 

impingement plate, a confinement plate from which a slot jet 

issues, a top spacer and also a bottom spacer. 

imeeTne mass transferring surface iS the hypothenuse surface 
of a right-angle glass prism (0.089 x 0.125m) which is 
fitted into a window cutting of an aluminum plate. The 
aluminum plate, 0.149 x 0.625m, has a thickness of 


0.012m. The surface of the prism together with the 
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THE OVERALL VIEW OF THE MASS TRANSFER 
EXPERIMENTAL SET-UP 
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aluminum plate form one flat surface. This surface is 
then coated with a thin layer of Silicone rubber to form 
the impingement plate. The front and top views of this 
impingement plate are shown in Figure 4.3. 

An aluminum plate having a thickness of 0.012m and the 
Same size as the impingement plate is attached normally 
to the nozzle of an aluminum slot tube. The aluminum 
tube with dimensions 0.0015 *% Dstis2Z % 084m ase long 
enough to provide a fully developed parabolic velocity 
profile at the nozzle exit for the highest Reynolds 
number used in this work. The large ratio of plate span 
(0.132m) to nozzle width (0.0015m) minimizes end effects 
and leads to a situation of nearly a two-dimensional 
jet. The front and top views of the confined plate and 
the aluminum tube are shown in Figure 4.4. 

Three sets of spacers are used so as to obtain three 
different jet-to-plate spacings. A set of spacers 
includes two aluminum plates of the equal width. One 
aluminum plate is mounted on top of the impingement 
plate and the confinement plate, so called top spacer. 
Another aluminum plate is mounted on the bottom, so 
called bottom spacer. The gaps between the impingement 
plate and the confinement plate for these three sets of 
Spacers are 0.003m, 0.006m, and 0.018m, respectively. 
These correspond to jet-to-plate spacings of 2b, 4b and 
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(a) FRONT VIEW 


taLaem 


(b) TOP VIEW 


FIGURE 4.3 : FRONT AND TOP VIEWS OF IMPINGEMENT PLATE 
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(a) FRONT VIEW 
Slot Tube 


40CmMm 


Confinement Plate 


(b) TOP VIEW 


FIGURE 4.4 ; FRONT AND TOP VIEWS OF CONFINEMENT PLATE 
WITH SLOT TUBE 
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The mass transfer experimental set-up is mounted on a 
heavy metal base as shown in Plate 4.3. The spacers are used 
to fit between the impingement and the confinement plates so 
as to form a rectangular channel in which the jet flow would 
be confined. With such an arrangement of the aluminum slot 
tube and the impingement plate, the air impinges normally on 
the impingement surface at a line which is 0.043m from one 
edge of the glass prism aad 0.082m from the other. Since 
only the mass transfer on the hypothenuse surface of the 
prism can be studied by holographic interferometry, the 
location of the slot nozzle allows the effective downstream 
distance for which mass transfer can be measured to be 55b. 

Before each experiment, the impingement plate is easily 
repositioned with the help of the top and bottom spacers. 
The impingement plate and the top spacer can be removed 
after each experiment while the confinement plate and the 


bottom spacer remained on the heavy metal base. 


4.3 EXPERIMENTAL PROCEDURE 


The experimental procedure for double exposure holographic 
interferometry is discussed in Sections 4.3.1 to 4.3.4, 
while the experimental procedure for real time holographic 


interferometry is discussed in Section 4.3.5. 
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4.3.1 PREPARATION OF POLYMER COATING 


1. 


The polymer coating is formed as follows: 

Liquid silicone primer, SS-4120, from General Electric 
is applied by bushing uniformly on the cleaned mass 
transferring surface. It is important to allow the 
primer to dry in air for at least one hour before the 
polymer is applied. 
Ten parts of liquid polymer rubber, RTV-615A silicone 
tubber, and one part of catalyst. RTV-615B (both from 
General Electric) are mixed throughly in a styrofcam cup 
tomaboutal0sil-kg of ym?xture. 

The polymer mixture is poured evenly on the mass 
transferring surface with the mass transferring surface 
being horizontal. 

In order to obtain a uniform coating thickness, the 
freshly coated surface is allowed to semi-harden at room 
temperature for about ten hours. 
The semi-harden surface is then cured in the oven at 
65°C for three hours. 


For complete swelling of a fresh polymer coating with 


ethylsalicylate it took at least ten hours, while 


re-swelling a partially exhausted coating after an 


experimental run required about three hours of immersion in 


the swelling agent bath. 


The refractive index of silicone rubber swollen to 


various degrees was found by MasSliyah and Nguyen (51). It is 


evident that the refractive index is a weak function of the 
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degree of swelling when the swollen polymer is near its 


equilibrium value. 


4.3.2 MAKING OF DOUBLE EXPOSURE HOLOGRAM 

Each full sized Agfa-Gevaert 10E75 AH photographic 
peare <0; 1016" 0.727m)" is cut) into twovhalvesi(Oet016 x 
0.0635m) in the dark room by using a glass cutter and an 
aluminum guiding plate. To avoid damage to the emulsion side 
of the photographic plate, an incision is made only on the 
back side with minimum stress. 

The 50 mW He-Ne laser is turned on at least half an 
hour before an experimental run. By rotating the mirror of 
the beam splitter-attenuator to the reading of 110, the 
intensity of the object beam is selected as about three 
times the intensity of the reference beam. The shutter is 
then set at 1/30 s. The selections of the exposure time and 
the intensity ratio between the two beams were found by 
trial and error, and remained the same for all the 
experimental runs. 

The impingement plate is removed from the swelling bath 
and is dried carefully by using clean tissue paper. The 
impingement plate is then placed on top of the bottom spacer 
and is mounted with the confinement plate and the top spacer 
to form a rectangular channel. The laboratory temperature 
and pressure are measured, and the laboratory lights are 
then switched off. The photographic plate is placed in the 


plate holder. 
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Since handling the mounting of the rectangular channel 
and the photographic plate can create a temperature 
gradient, the equipment is left for several minutes in order 
to equilibrate with the laboratory. The shutter is then 
activated with an exposure time of 1/30 s and the first 
exposure 1S made on the photographic plate. The mass 
transfer experiment is then started. Air is supplied from a 
compressed air cylinder and its volumetric. flow rate is 
measured by a rotameter. Two rotameters are used in this 
work, Fischer & Porter Co. Rotameter and Brook Rotameter. 
Each rotameter covers a different range of volumetric air 
flow rate. The calibrations of these two rotameters are 
given in Appendix B. After a given period of time, the air 
flow is stopped. The shutter is reactivated with the same 
exposure time to obtain a second exposure on the 


photographic plate. 


4.3.3 PROCESSING OF HOLOGRAM 
The exposed photographic plate is removed from the 
plate holder for processing. The processing sequence 1S as 
follows: 
The photographic plate is 
1. developed with Kodak D-19 developer solution for 2 to 3 
minutes, 
2. stopped with Kodak indicator stop bath for half a 
minute, 


3. fixed with Kodak rapid fixed solution for 2 minutes, 
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washed with tap water, 
and finally dried by spraying methanol on the 


photographic plate. 


4.3.4 PHOTOGRAPHING OF RECONSTRUCTED IMAGE 


'To record the virtual image of a double exposed 


hologram, the following procedure is used. 


i. 


Br 


4. 


The double exposed hologram, after processing, is 
repositioned in the photographic plate holder. 

The intensity of the reference beam is selected as about 
one half of that of the object beam by rotating the 
mirror of the beam splitter-attenuator to the reading of 
60. In other words, the intensity of the reference beam 
is higher than that when the double exposed hologram was 
made. 

The object beam is then blocked off. 

A Nikon F2A Photomic camera with Nikon Motor Drive MD-3 
and Micro-Nikkor 105mm £/4 Lens is placed behind the 
hologram and the view guiding aluminum plate. The lens 
aperture is normally set at f£/4 and exposure time varied 
from 8atot20eseconds. “Focusing, the camera,on the 
reconstructed virtual image, the photograph of this 
image can then be taken. 

Kodak. Tri=X: pan film. (TX 402)..with.ASA.speedeof 400, a 
fast black and white film, is used. 


The, processingsof, the.Tri-X pan. film is «proceededsin 


room temperature as follows: 
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The film is 

1. developed with Kodak D-76 developer solution for 7 
minutes, 

2. “stopped” with Kodak indicator stop bath-for*half°a 
minute, 

3. fixed with Kodak fixer solution for 4 minutes, 

4. washed with running tap water, 


5. and finally dried by hanging in a dust free room. 


4.3.5 PROCEDURE FOR REAL TIME HOLOGRAPHIC INTERFEROMETRY 

In this work, real time holographic interferometry is 
used for qualitative study only. This technique is used for 
determining the order of a particular fringe of the "frozen 
fringe” pattern obtained from double exposed hologram. 

For real time holographic interferometry, an 
immersion-type X-Y micropositionable photographic plate 
holder, model MPH-45W, is used instead of the simple plate 
holder used in double exposure holographic interferometry. A 
full size Agfa-Gevaert 10BE75 (0.1016 x 0.127m) photographic 
plate is placed in the PC-45 plate carriage and is allowed 
to normalize in the plate holder water cell for 
approximately 15 minutes. This operation pre-sensitizes the 
emulsion and also allows strain release, swelling, or other 
emulsion excursions to occur prior to exposure. The mass 
transfer set-up is prepared as previously discussed. The 
shutter is activated with an exposure time of 1/60 s. This 


exposure time is half of that of the double exposure 


- 


Py tei iy age 2 eg ag pos SNS 

wine ee ao Peete. asian” ieee’ “an 

DS asta, Wee p Gees se hens die oi 

mig? sos isu oat Qapae aa ¥ea oc ‘ital 

TAPAMORE ARYA T  DIRAATOONION: ag? then ads Stueaooi 

27. ya tometsiss30 i. haga teas a emad Ther’ oe adds 

jaya? supsadest eid? phn eu ts Be is3 ttsig"y 14 

rgaere” at? Yo sem Si zalcoreiea aiid ohio ~ efi /pabnkae 
ith Sele EBBOGks* > Ties ees 'S vei? aide eT 

fe ¥isanoeetestad “widguayolcs amie tess" 

atthig- sinteieadeda vidanotas igeuerrs im rit ‘aque 

minis ofgite bie. is ‘pee epi ‘Bead at (wae reson 


“ 2 4 
A Widow erage} sPaigareot ot subactas aiduok a 5seu es 
re 


shigerpotcdg: (an BE ude haliial s0) @UROH. sreaved“eiiga 
bevolis #7 tim aget tino ajaltq 86-374 427 ni ssa 
| io2 itso” té6os0 seBiod eietg “ade at asl fs 

a2 oouh sieves -ony Wdisaregs we2T .@asenia dF tis: 
redec yO ,eefiiawe yeasotex cieste awelta celia! 
220m adT SWWEEGRO OF wig TWase 7 enoisniose'n 
‘ott ,heveunei&b yisvetve:q se Beréeni1g af qu-s6e | . | 
ceiak- se bt tea odivasld setae etal ne dsie basal sam: - 194: ~— 


* 


65 


holographic interferometry because the photographic emulsion 
becomes more sensitive while soaking in water. But the 
intensity ratio between the reference and the object beams 
is maintained the same as that of the double exposure 
holographic interferometry. The plate carriage with the 
photographic plate is then removed from the water cell for 
processing. The processing sequence is as previous Section 
4.3.3 except in this case the developing time is shorter. 
The hologram after processing is placed back into the 
water cell to within one fringe alignment by using the 
precision X-Y controls on the plate holder. With both the 
reference and object beams striking the hologram, the 
fringes which occur between the holographic virtual image 
and the actual subject image can be eliminated by adjusting 
the plate position with the precision X-Y controls. When no 
interference fringes on the object are observed, this 
indicates that the photographic plate is repositioned in its 
Original position. The mass transfer experiment is then 
proceeded with. The shrinkage of the polymer coating can be 
viewed behind the hologram as real time interference fringes 
across the image. The order of any particular fringe can be 
determined by counting fringes at a given location on the 


impingement plate as they appear in real time. 


ode dus port “ik Lana doh tec 
zines ©! Toatae ede. baie oatpaatey dit waiienad * oli 
eee ee od ala benket 

at miv agaasiag aeeit eat sc rdoare oan pe 
12 ay aha Set pug seecomag: satsets: ast lnaaetig 2 
no lg nb, sum evesS gS onshore vase WT a 
a4e7G2@ 21 one cineca tas Sb sy id ae ne ae 
pis atx; ‘pee fera-g + oleae whe By a) MBIPOS 
ois sniaw goosebntp ele grew tsz: efe pitt of i 

ot4 fed C2ie..ceh lol-saaes y afk ae ahovsaes 39 
wad. = motgoran ats pris ieee seater bis. 9 | 
sogn? Japiary aeeieooicd ants wWao0 “tea a 
mridentia ye ha narrate! a oat. peet apse’ “ osatdue: nya 
on wed? ,ehbwertoe Prd nasectwetey. on ee mop taiony 
23603 .bevisedse Ge jas de sg a sngntes : 
sci ci Henobel doghs HE weary wstqayheseds aif? etd | “9 ag 
end »i 2edni tags teraans 2 daa ee” silaiabaehciaas Pr 
od mea gnitega Tey Ley ets to sQperaiuia salt. -daiw i 
rapt seg Saou eA suis fest Seeergoled 242.4 


“ te npeits®. raids, wre to 396% oF. - ssid sa 


‘rr? 


etd 1 notmecot ewig « ¢@ sepnizt pnsamuns de oo le: 
aaa: cee se no eae 


Ss ' ne = a > 


7 a - 


5S. CALIBRATION AND VALIDITY OF THE EXPERIMENTAL SET-UP 


Mass transfer characteristics due to an unconfined 
laminar impinging two-dimensional air jet is studied by 
measuring the local Sherwood number along the impingement 


plate. Local Sherwood number in this work is defined by 


Sh, kb /D (bie 


or Sh 


tf} 


pRSet RDS) 0 (5 $2) 


where k and k' are the local mass transfer coefficients, and 
Pers the diffusion coefficient. The Jocal mass. transfer 


coefficients k and k' are defined by 


si) oI O2/ asl at Mw / P (5ri35) 


and N k' (p, - P,) p®? Mw / P (5.4) 
where N is mass flux, p°® is the molar density of the gas 
mixture, Mw is the molecular weight of swelling agent, and P 
is the total pressure. The partial vapor pressures of the 
Swelling agent, Pr P and p, are located at the coating 
Surface, the jet nozzle exit and in the bulk flow, 
respectively. 

For an igontdnsa ainjeth the “partial vapor pressures 


of the swelling agent at the jet nozzle exit and in the bulk 


flow are equal to zero (p =p,=9). It leads to k=k' from 
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Equations 5.3 and 5.4. However, for an confined air jet, P, 
is no longer equal to zero and it is not easy to measure in 
the experimental set-up of this work. Therefore, only the 
Sherwood numbers evaluated from Equation 5.1 are 
investigated in the experimental study for confined air jet. 
On the other hand, p, can be evaluated numerically. Both the 
Sherwood numbers evaluated from Equations 5.1 and 5.2 are 
investigated in the numerical study. 

For both confined and unconfined air jet, Dae and the 
partial vapor pressure of the swelling agent at the coating 
surface, Poy is that of its vapor pressure, P°, during the 


"constant rate period" and Equation 5.3 simplifies to 
N = k B® °° Mw/ P (S85) 

The mass flux, N, can also be given as 
Nee 4c toe aaah (5...6.) 


where r' is the change of coating thickness given by 
Equation 4.1, Ps is the density of the swollen polymer and T 
is the duration of the mass transfer experiment. The product 
rip. is the mass of swelling agent transfered per unit area, 
since the swollen polymer volume is the sum of the volume of 
the dry polymer and that of the pure swelling agent (47). 


Substituting Equation 4.1 into Equation 5.6, one obtains 
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Ny = yD Ag Pou, T (5.7) 


Define A = B} pats then Equation 5.7 becomes 


Naas 4p) YT) | (5.8) 


where A is a constant for a given mass transfer experimental 
and optical set-up. Constant A is also referred to as the 
calibration constant of the experimental set-up. Combining 
Baueations 5S. th. 5.5 andy 5.8,) one. obtains 


Sh (AG.Piy/ 20° PoeMWve tb «f Divin: / Tp 


Ge (na /.T) (5.9) 


where G is dependent only on the physical properties of 
working fluids. 

Therefore, to evaluate the local Sherwood number it is 
necessary to know the physical properties under the 
operating conditions, the local fringe order and the 


duration of the experimental run. 


5.1 CALIBRATION PROCEDURE 


In order to evaluate the mass flux of the swelling 


agent, it is necessary to know the value of the constant A 
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of Equation 5.8. The calibration constant A can be obtained 
by studying the mass transfer due to an unconfined laminar 
impinging axisymmetric air jet with initial parabolic 
velocity profile using the same optical set-up as shown in 
Figure 4.1. 

This calibration is made using the theoretical 
expression Obtained py scholtz-and’ Trass 479) for’ an 
unconfined laminar impinging axisymmetric submerged jet in 


the wall jet region given by Equation 2.15 as 
SDipeprte WG yi SISA inet ia Gab Asli iis aloe CSiat 0) 


meereercy ="0.,e50 ie(scet 175 y7 tl Sc) 1 (17 o2) Shes ba 

FOr initiad parabolic velocity profile. Equation’ 5.10 has 
been shown by Scholtz and Trass (79) to be fairly accurate 
in’the wall*jet region.” Equating’ Equations 5.9 and 5.10 by 
using the nozzle diameter, d, instead of the slot width, b, 
asithe characteristic length in’ Equation’ 5.9 for"the case of 


axisymmetric and rearranging, yields 


neem(eu Dp pear? Mweyedoa DT Re, e-PRtryd)oe**) 


(5712) 


If a mass transfer experiment is performed using an 
unconfined impinging axisymmetric air jet with initial 
parabolic velocity profile, then a plot of n versus 
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Slope (c, D 0° P* Mw / d A P). The constant A can then be 
evaluated from the knowledge of the physical properties of 
the system. 

The experimental set-up and procedure for an unconfined 
axisymmetric jet are similar to those for a confined 
two-dimensional jet mentioned in Chapter 4. Here, the 
confinement plate and slot tube referred to in Chapter 4 is 
replaced by a circular Piss with a diameter equal to 0.003m. 
This tube is mounted on a heavy stand. The tube is set 
perpendicular to the mass transferring surface with a 
jet-to-plate spacing of 1.5d. Double exposure holographic 
interferometry is used throughtout this calibration. 
Experimental runs with durations equal to 90s, 180s and 360s 
and for Re y=1210 and 1470 are made. The average operating 
temperature and pressure of the experimental runs were 
20.7°C and 93.87kPa, respectively. 

Themusrozenatuinge, pattenn forerunyno., Ciil4ei1. 5B for 
Re,=1470 and T=90s iSushnown un Plate 5). la.) The wsirozen 
fringe” pattern for run no. CJ14-6A for Re ,=1470 and T=360s 
is shown in Plate 5.1b. Determination of fringe order of 
these "frozen fringe" patterns is rather simple. In Plate 
S.la, the outer bright region is considered to be the zeroth 
order fringe where no mass transfer as yet has occured. 
Since in this case the zeroth order fringe is known, the 
determination of the order of other fringes is 


Straightforward. By counting the fringe order from the outer 


zeroth order fringe towards the stagnation point, it is 
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PLATE 5.1 


(b) T = 360s 


CONTOURS OF EQUAL MASS TRANSFER RATE 
AN UNCONFINED AXISYMMETRIC AIR JET 
(Reg=1470) 
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possible to label all the fringes. It is thus this short 
duration experiment that gives the base on which the long 
duration experiment, such as the one shown in Plate 5.1b, 
are interpreted. In each experimental run, the variation of 
the local fringe order, n, with dimensionless radial 
distance, r/d, is determined. Experimental results of local 
fringe order, n, and dimensionless radial distance, r/d, are 
given in Appendix C. 

A confirmation of the "constant rate period" can be 
obtained by combining all the experimental runs for a given 
Reynolds) numbermin the form of n/T (51). A plot of n/T 


versus r/d for Re,=1470 is shown in Figure 5.1. Here, the 


d 
Gata for different duration collapse onto one curve 
indicating that the driving force, in this case, the 
swelling agent vapor pressure at the coating surface, P°, is 
the same for both the shortest and for the longest mass 
transfer experiment. 

ASpLOt OfenQuation= 521251 s.'given an “Frqures 5.2.) The 
slope of the least squares linear regression is 634.08s™' 
which gives Avas 1-6 149105 “Mko/m>using-physical properties 
in Appendix A evaluated at average operating conditions of 
20.7 Geand, 93..67kPa.. 

It is possible to proceed in a different manner to 
obtain a calibration constant A. Using the appropriate data, 
Equation 4.4 gives B=0.2615. Since the calibration constant, 


A, is defined as Bo., by using the values of swollen 


coating density, epg, and the wavelength of the di ghty yn 
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FIGURE 5.1 : VARIATION OF n/T WITH DIMENSIONLESS 
RADIAL DISTANCE FOR Re ,= 1470 
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FIGURE 5.2 : CALIBRATION PLOT USING AN UNCONFINED 
AXISYMMETRIC AIR JET 
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Appendix A, A becomes 1.67x10-* kg/m?. 

The agreement between the two methods for evaluating 
the constant A is within 3.3%. This agreement confirms the 
success of the method based on optical analysis mentioned in 
Section 4.1. This also indicates that the coating surface is 
reasonably smooth throughtout the region of interest during 
the experimental period. The value of A as obtained from an 
unconfined axisymmetric ae jet is used for the evaluation 
of local Sherwood number throughtout the experimental study 


of this work 


5.2 TEST OF CALIBRATION VALIDITY 


The validity of the value of calibration constant A 
evaluated in Section 5.1 is tested by studying the mass 
transfer due to an unconfined laminar impinging 
two-dimensional air jet with an initial parabolic velocity 
profile. The same optical set-up as shown in Figure 4.1 is 
used. Local Sherwood numbers in the wall jet region 
evaluated by using Equation 5.9 with A=1.614x10°* kg/m?’ are 
compared with the regression equation given by Masliyah and 


Nguyen (53) 


= 0.55 = 6,73 ya 
Shy. 0.55 Re, (x/b) (5 ) 


for an unconfined impinging two-dimensional air jet 


eer ne ea | — 
| maa seater ee abit Re oe es arenes: | 
arts rt cn Tests abt Se nidyi vier ave 
pacgidaen abeyiare iapbsagene beaded bodsom-. oa 
vue pataeee ard. sent YWesietbat ‘pale aint .t 
val woe: ade tunwhga ets dscome- 


par Jae 


hen pa? “t+ 2a of a - aq hee wit ieeaag wen 


oisgu leva ada p03 et: ag % “ie os stenwys ine fe J 
ae | a 

bie tosdeivarea wet yh cefaasee AT 2 pects Boow she - 48 
“a on 
ay - 


; = 


yTstdat HoT TART IAD 101 124 


i An) 
’ 
A a 


" | r ; 4 


rszeineo ndi¢es@iles Jo site nas 3° qth baw 


A 7hke a 


es paivybese ya bestest adi 8 satdont ts 


& 


exipaige! seckmnd bent idesau a6 03 sub 3 


~~ 


vyiiseltav siiagwmegG Ssivint ag Atiwv wt Sea tenoler 


ai t..¢ siwete we avawe ea, cay 3er fag 3qo 21g ott tha 


noha pe bicawe giz ti eiedaun Ssoedede 7 
eo Cael 


wa "aren br aera ea Hetw ez nolturet patae yd Sosy 


bas devilogt vd nay eottcups coizaezps: sd2 4atw 


ma 1 it 
ere Cd\ x) 54.9 af 227.0 * dz 
aT ; e re ¢ c “7 pao 


teneianenid- ows tiny nh 
7 _ | b 40 
: 7 Se 


76 


(Sc=2.85). 

; The experimental set-up and procedure for an unconfined 
two-dimensional jet are similar to those for an unconfined 
axisymmetric jet mentioned in Section 5.1, except a 
plexiglass slot tube with dimensions 0.0015 x 0.075 x 0.35m 
is used. This tube is mounted on a heavy stand and set 
perpendicular to the mass transferring surface with a 
jet-to-plate spacing of 4b. Double exposure holographic 
interferometry is used. Experiemental runs with durations 
equal to 120s, 240s and 480s and for Re, =94 and 204 are 
made. The average operating temperature and pressure of the 
experimental runs were 20.7°C and 93.98kPa, respectively. 
Under these operating conditions, Sc=2.74. 

The @frozen arrange satterny Lon mun ono maSTia=Z2h “for 
Re, =94 and T=120s is shown in Plate 5.2a. The "frozen 
fringe” pattern forsrun no. Sv i-4r tor Re, =94 and T=240s is 
shown in Plate 5.2b. In both cases, the zeroth order fringes 
are no longer on the mass transferring surface. A simple 
trial and error method introduced by Masliyah and Nguyen 

(51) when the zeroth order fringe cannot be identified is 
used to determine the bernde order in. this ‘case. In their 
study of unconfined axisymmetric air jet, Masliyah and 
Nguyen pointed out that, according to Equation 5.12, a plot 
of log(n) versus log(r/d) for experiments with various 
durations of a given Reynolds number gave curves which are 
parallel to each other. Similar behavior of fringe order 


with dimensionless streamwise distance for the 
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two-dimensional case is expected. By choosing the correct 
tocal fringe order, the curves from a plot of log(n) versus 
log(x/b) for experiments with various durations of a given 
Reynolds number should therefore parallel to each other. 
Such a plot for an unconfined two-dimensional air jet is 
shown in Figure 5.3 for Re, =94. Using such trial and error 
method, the outer darken Peyton's in botbePlates So .canaud 
OeeD are determined £6 be the first order fringes. By 
counting the fringe order from the outer first order fringe 
towards the stagnation point, it is possible to label all 
the fringes. AS soon as the local fringe order is known, 
local Sherwood number can be easily determined by using 
Equation 5.9. In each experimental run, the variation of 


local Sherwood number, Sh with dimensionless streamwise 


b’ 
distance, x/b, are determined. Experimental results of local 


Sherwood number, Sh and dimensionless streamwise distance, 


b? 
x/b, are given in Appendix D. 

It is worthwhile to mention that the regression 
equation given by Masliyah and Nguyen (53), Equation 5.13, 
Mswonmly walid tor Sca=2..85, which vs slightly different from 
Phe sSc=2474 of this work. Considering the effect of Schmidt 
number on local Sherwood number by using the correlation 


given by Scholtz and Trass, Equation 5.11, Equation 5.13 


becomes 


Sh. = 0.54 Re,’ °° (7 Oiedoee C5 314) 
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FIGURE 5.3 : VARIATION OF FRINGE ORDER WITH 
DIMENSIONLESS STREAMWISE DISTANCE 
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for Sc=2.74. 

A plot to test the validity of the experimental local 
Sherwood numbers for an unconfined two-dimensional air jet 
is shown in Figure 5.4. For a perfect fit with Equation 
5.14, all data points should lie on a straight line having a 
Slope of unity. From Figure 5.4, although it is observed 
that the data points are consistently above the perfect fit 
line, there are still within the scattered range (+9%) of 
the experimental data obtained by Masliyah and Nguyen (53). 
Therefore, the validity of the calibration constant A is 


reaffirmed. 
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FIGURE 5.4 : COMPARISON OF MASS TRANSFER DATA DUE TO 


AN UNCONFINED IMPINGING TWO-DIMENSIONAL 
AIR JET WITH EQUATION 5.14 
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6. NUMERICAL MODEL-MATHEMATICAL FORMULATION 


A two-dimensional numerical model is used to simulate 
the experimental set-up. The impinging jet system considered 
in this work is shown in Figure 6.1. The air jet issues from 
a two-dimensional slot tube of width b with an average 
velocity of Vi. The confinement plate is located parallel to 
and at a distance h from the impingement plate. The 
impingement is normal to the impingement plate. For the 
description of the flow and the concentration fields, a 
two-dimensional rectangular coordinate system is used with 
the origin at the center of the jet nozzle exit. The 
x-coordinate is parallel to the impingement plate and the 
y-coordinate is normal to it. The outflow region is chosen 
at a location sufficiently far away from the stagnation flow 
region to ensure that the velocity and the concentration 
profiles at this location are developing as those for 
parallel plates. 

The governing equations for this two-dimensional 
impinging jet system are presented in Section 6.1. The 
boundary conditions for this system are presented in Section 


6.2. 
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FIGURE 6.1 


: COORDINATE SYSTEM AND BOUNDARIES OF THE 
IMPINGING JET SYSTEM 
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6.1 GOVERNING EQUATIONS 
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The pressure gradient is so small that the change of 


density at each point in the system can be neglected. 


Therefore, air can be treated as an incompressible Newtonian 


PLUG, 


even though it 1S a compressible fluid itself. The 


two-dimensional momentum and transport equations can be 


reduced to the corresponding vorticity-stream function form 


with the assumption of steady state, 


incompressible viscous 


Newtonian fluid flow with constant physical properties. The 


pertinent equations are (69): 


9(u w) 


ox 


a(v w) 


oY 


where w is the vorticity and is defined by: 


dW 
) (6.1) 

oy? 
(67923) 

ac 
peek (6.3) 

oy’ 
(6.04) 


» is stream function and is defined by 
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— = U and aaae Gaels ev (G75) 


c is the molar concentration of the swelling agent, and R is 
the molar rate of production of the swelling agent per unit 
volume. For a system without chemical reaction, R=0. 


When dimensionless variables are introduced as follows: 
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a(U C) a(v Cc) 1 d7C yeh & 
ae ee eee ee | ) (6.9) 
ox OY Re, nol ole nT Pen TD Ge 
au au 
where — = U and — = -V (6.10) 
oY ox 


The objective of the numerical study is to solve 


Poauet 2ons 6.7, 6.8-and 6.9 for, 2, Pand.c. 


6.2 BOUNDARY CONDITIONS 


Because the governing equations are elliptic in nature, 
boundary conditions must be specified at all the boundaries. 
The boundaries are classified in five regions: the nozzle 
exit, the confinement plate, the impingement plate, the axis 


of symmetry and the outflow region. 


S.2. 1 NOZZLE EXIT 
Fonvam initaa) iparabolic’ velocity profile atthe nozzle 


exit, the velocity components are 
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From these velocity components, one can show that 


b= -1.5 % + 2 * (6.12) 


ta 


/ he 
| mw? neko. 
2%, oe n sy * ‘eal 


sipjen az - faiite aap end i sawp guinsevag afd, Bhi Sr 
.asiveboued adv iis 7a bai st gage ed faum aac 
sigue agz tehotean avit ab Gertreneds s16 waive 
abne ‘ets .s2atg Sdenepaigar = ree snemani dade: 
een siiahinaat ais baw: 7) 


v PINS BR: 


gixson odd zs silteve gaioeior Gidadetaq feiiink ae 


ae atacdogwss sttootey 


er eee oe ee | 
le ms rea ie 


- 7 ¥ ' 
7 : . oe i / 
ie jie os 


— ahaaaqnoa eoaaley sna + me 
ay ; nf ion oa we ae une 


7 - 


87 


N= 12 xX (6.13) 


where the stream function, $, at the axis of symmetry (X=0) 
is taken as zero. 
Fo@ene initials tlat velocity profile at nozzle exit. 


the velocity components are 


U = 0 (6.14) 


W = -X (65,15) 


Q = 0 (6.16) 


The boundary condition for the concentration at the 


nozzle exit is 


C = 0 (6547) 


for an air jet. 


6.2.2 CONFINEMENT PLATE 

Because the confinement plate is impermeable, the value 
of the stream function does not change along the plate. This 
value can be determined by substituting X=0.5 in Equations 


6.12 and 6.15 for the case of parabolic and flat velocity 
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profile, respectively. For both cases, the boundary 


condition of the stream function at the confinement plate is 


Y= -0.5 (6.18) 


The boundary condition for the vorticity at the 
confinement plate is evaluated by using no-slip boundary 
condition (96). This pouane condition is given in finite 
difference form and will be discussed in Chapter 7. 


The boundary condition for the concentration is 
aC / oY = 0 (6.19) 
due to no mass transfer occurring at the plate. 


6.2.3 IMPINGEMENT PLATE 

Again, because the impingement plate is impermeable, 
the value of stream function does not change along the 
plate. The stream function can be arbitrary set to zero to 


give 
b= 0 (6.20) 


The boundary condition for the vorticity at the 
impingment plate is the same as that of the confinement 


plate. 


> . au hike | 

edt 2h evisistev-ese ves inthe telog genabs 
igi giles cae a ele, at atsig 4 
| ee “ ea SD yeenyed od watt sey 
tesqBe) a: Ansanoel hee il ae baie — SOR 
“eo sie cal” Aeteibaos yistewed 3 


Dy an 
> i. = 


(VW: wal | 4) = T¢ 5‘ 2t 
i 
; jatata ef? ta embetegse Seagrass way: on. 
we ; eit 8 . oes G 


7 i 
~<A THIEME Mie 
alcuemiggmi 2] einke iogerneraryys ot enusaad ‘ : 


#43 eavels wteic 100 aaah, notaetat neo1 se ‘Yo sul 


13 ovat 09 240 peeuhide a. deo folloust maaaze SdF 
5" eh be ee : 


(95.91 
te de tava: ves ae ot sheng ‘oad 
- ae | fe 


7 re 
Yo jokt ow oone “ia See 
ee a eae ee ar 


ed al Pin 


_ - 
if 7 


89 


The boundary condition for the concentration is given 


by 


Ere 4 (6.21) 


6.2.4 AXIS OF SYMMETRY 
The axis of symmetry.1sS given by X=0 and the stream 
Monecion 1S avcOonstantvarong™ rt." Arbitrarily the €onstant is 


set to zero to give 


b= 0 (6.22) 


The boundary condition of vorticity can be easily 


determined by setting U=3V/3X=0U/0Y=0 and noting that 


V(X,Y)=V(=KVY) 2) The vorticity along x=0 is then given by 


2 = 0 (62230 


Since the concentration is also symmetric on both sides 


orVrhe=ax7s> *1t*leads ts 


00a /s0kt =" 0 (6.24) 


along the axis of symmetry. 
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6.2.5 OUTFLOW REGION 

This boundary is located sufficiently far from the jet 
such that the flow is nearly fully developed where the 
pniilvence,ofethe impinging. zone is «not, felt. Usually, a 
fully developed flow is assumed at this boundary (96). If, 
instead a developing flow profile for the parallel plates 
channel is assumed, the outflow boundary can then be located 
anywhere in the x-direction as long as it is far enough and 
not influenced by the impinging jet flow. With an assumption 
of a developing flow profile at the outflow boundary, this 
boundary can be located closer to the stagnation point than 
that with an assumption of a fully developed flow profile. 
inmeOther,words, the. nuniber.of.gridipoints.in x<direction can 
then be minimized. 

The developing stream function and vorticity profiles 
are derived from the developing velocity profile given by 
Sparrow et. al. (92) for their study of velocity development 
for the parallel-plate channel. By using only the most 
dominant term in their solution series, the stream function 


and the vorticity at the outflow region are given by 


1 sin [a,(2(¥/L)-1)] Ts L 


fd a,?L COS a, 20 2 


e (6.25) 
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2 Sin [a,(2(Y/L)-1)] ~Ba,?(X/L) /Re 
Q= QQ - {———_—___———--} e 2 
td a,L? cos a, 
(6.26) 


where a,=4.49341. The subscript "fd" denotes fully developed 


flow, where 


le OeGY bom tN/ ob) = 05 C6327) 


oa 
" 


fd 


" Sonetea momen t)s), 6/052 C6e2on 

The developing concentration profile is analogous to 
the developing temperature profile for parallel plates given 
by McCuen (62) and Shah and London (85). The boundary 
conditions used in this work are similar to those of the 
fundamental solution of the third kind in their work. The 


concentration at the outflow region is given by 


= —2,?(X/2L)/(Re, Sc) 
ee 


Cuenc + Ye (6.29) 
fd t 1 


i=1 


where Y,'s are nlunctiyons. Ory . E.'s and A'S are 
eigenconstants and eigenvalues, respectively and they are 
evaluated at the confinement plate for the first four terms 


in the series “as (62,85): 
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Due to the oscillating behavior of the series in Equation 
6.29 (from the values of E.Y,'s), there is no dominant term 
in this series. The fully developed concentration profile is 


given by 


C = 1 (653 1) 


The finite difference form of the boundary conditions 
in the outflow region will be discussed in Chapter 7. 

The location of the outflow region boundary is chosen 
depending on the jet Reynolds number and jet-to-plate 


Spacing. It 1S given in Table 6.1. 


TABLE 6.1 : LOCATION OF OUTFLOW REGION BOUNDARY 
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7. NUMERICAL FORMULTION 


Initially the numerical technique used in this work was 
Similar to the method introduced by Joseph, Smith and Adler 
(35) which is that of the Marker-and-Cell (MAC) method used 
by other investigators dealing with numerical studies 
eo, 28,30,65,72, 93,99): The method is to solve the unsteady 
State primitive equations. A steady state solution is 
Obtained by advancing the velocities and pressure from one 
time interval to another until the solution no longer 
changes with time. This numerical technique was later 
abandoned due to the failure to obtain a converged solution 
for higher Reynolds number unless the time interval was 
Steadily reduced. Large CPU time was required in order to 
obtain a steady state solution. 

The second numerical technique used in this work was 
the central finite-difference representation of the steady 
state vorticity transport equation. This numerical technique 
waS again abandoned due to the failure to obtain a converged 
solution for Re,? 100. Severe under-relaxation was required 
for the run of Re,=100. Converged solution for Re,=100 was 
not obtained until after 2700 iterations. 

The third and successful numerical technique used in 
this work is the hybrid differencing schemes, so called 
"upstream-weighted" and “upstream” differencing schemes, 


introduced by Raithby and Torrance (69). The detail of 
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derivations of the general finite-difference equations is 
given in Section 7.1. The derivations of the 
finite-difference equations at the boundaries are given in 
Section 7.2. The finite-difference equations for 
"upstream-weighted” and "upstream" differencing schemes is 
discussed in Section 7.3. Finally, the stability properties 
of the finite-difference equations are discussed in Section 


1.4. 


7.1 FINITE-DIFFERENCE EQUATIONS 


The governing equations are: 


a(U o) a(v >) a? a? 
stl Mal gala al i (Tes) 
aX ay aX? ay? 
07 a? 
v ‘ On ar () A) 
dX? ave 


where — = U and — = -V (e7ne3:) 


Equation ee a more general equation than Equations 6.7 
and 6.9. The variable ¢ represents vorticity, 2, in Equation 
6.7 and concentration, C, in Equation 6.9. The coefficient a 
for these two cases becomes (Re,)~' and (Re, Soom, 


respectively. 
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The region of interest is represented by a rectangular 
grid network shown in Figure 7.1 with grid lines parallel to 
the X and Y coordinates. Grid lines in X and Y directions 
are designated by i and j, respectively. As shown in Figure 
7.1, the total numbers of node (the intersection of the grid 
lines) in X and Y directions are nx and ny, respectively. 
The numbers of node cover the jet nozzle exit is nj. 

For the time being, let us restrict Our attention to 
the region surrounding the typical node (i,j) shown in 
Becure 1,2. The variables @, WY and C are defined at each 
node, while the velocity components, U and V, are defined at 
points midway between these nodes. In other words, the 
velocity components are on the boundaries of the control 
volume. The control volume of a node is the region bounded 
by the dashed lines shown in Figure 7.2, the sides of which 
lie midway between the neighbouring nodes. One of the 
advantages of using sucha finite-difference grid-is that it 
simplifies the computation of mass flux into and out of the 
control volume, since the velocity components are located at 
the control volume boundary itself. 

The finite-difference equations of the velocity 


components are obtained as follows 


CW(it1,j+1) + WCi, gti) - WCi+1,j-1) - Yi, j-1)1/4 
SU / 2a) AS A) (7.4) 
[w(itt,j+1) + H(ie1,j) - WCi-1,j+1) - PU-1,5)1/4 


= - V(i,j+1/2) 4x(i) (7.5) 
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such that U(i+1/2,j) is the mean velocity in the X-direction 
in the region bounded by grid lines paSeinaeneetas nodes 
ee a) Vi, je uel, Jo) and. (1+1, 3-1), and Vis jei72%is 
the mean velocity in the Y-direction in the region bounded 
by the grid lines passing through nodes (i+1,j+1), 
Sieve i) e(i=t /A)eandi(i44,5). With these approximate 
equations, the mass balance is satisfied exactly over the 


control volume. From Equations 7.4 and 7.5, one obtains 


CU Gi ze ene 1/285 )' AY C5) 


+ [V(i,9+1/2) - Vli,j-1/2)) xi) = 0 (7.6) 


Dividing by AX(i) AY(j), Equation 7.6 becomes the central 
finite-difference representation of the equation of 
continuity over the control volume of node (i,j). 

The convective term 3(U $)/9X of Equation 7.1 is 
represented as the difference between two fluxes, one 
through the right hand side boundary, the (i+1/2) face, and 
one through the left hand side boundary, the (i-1/2) face, 
of the control volume shown in Figure 7.2. The flux, Ug, 
through the right hand side boundary is approximated by 
U(it1/2,j)P(it1/2-a,,j), where g(i+1/2-a,,j) is equal to the 
"true" value of the variable ¢? half-way between nodes ieee 
and (i+1,j) and is located at (i+1/2-a,,j) if one linearly 
interpolates between O(i,j) and O(i+1,j) (69). Since both 
(i,j), O(it+1/2-a.,j) and @(i+1,j) lie on the same straight 


line, therefore 
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Proceeding similarly at the left hand side boundary, using 
the parameter Or the flux, UO, is approximated by 


Pitt 2,3) 0(1-1/2-0, 4) where 


(0.5 + a) O(i-1,5) + (0.5 - ay) (4,5) 
Sai 1/2 een ) (7.8) 


Therefore, the convective term in X-direction becomes 


a(U >) 
AZ “> Azo? 
aX 


UN / 2h) Oa ea eg eaZeO kK) C229) 


The convective term a(V $)/aY of Equation 7.1 is 
represented as the difference between fluxes, one through 
the top boundary, the (j-1/2) face, and one through the 
bottom boundary, the (j+1/2) face, of the control volume 
shown in Figure 7.2. Similarly, the convective term in the 


Y-direction becomes 


a(V ) 


dY 


= [V(i,j+1/2) (i, 5+1/2-B,) 


- V(i,j-1/2) O(i,5-1/2-B,)] / 4¥(j) (7.10) 
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where (0.5 + B ) ONG 903 i ees 18 as B) O(i,j+1) 
= O(i,j+1/2-B.) J cml td} 


CORSee SRO ena vet (0. on58. ob ina) 
= Oli, 5 228.) (7242) 


Note that all a's.-and B's are unknowns and are to be 
determined. | 

The finite-difference equations for the diffusion terms 
a(a?/oxX?) and a(a?/aY?) of Equation 7.1 were introduced by 


Raithby and Torrance (69) as follows 


a7 a OM a er Olay iD 
a. = Pele yo) La 
dX? AX(i) AX CY) 
Oi pee OP Ge tye 
SIN UA Ure means wears seperate kd (7135 
os AX(i-1) 


in the X-direction; and 


ne a OG Jay eld) 
a ‘ = cleo) 7 ae 
oY? AY(3) AY(3j) 
Aa) ry Oral, 
=EAte—=6.) 2 : ad : hg eye aa 
o AV Cj = iy 


in the Y-direction, where y's and 6's are weighting factors 
to be determined. The advantage of using these weighting 


factors is that they retain flexibility similar to that for 
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the convective terms. 

Considering the flux crossing the common boundary 
between two adjacent control volumes, conservation requires 
that a, and y, for the control volume of node (i,j) must 
yeaa ond Fistor Ene*control volume of the node (i+1,3). 
Similarly, pevand see sforathe*control volume ‘of ithe node 
(i,j) must equal 8) and 6 for the’control volume of the 


node (i,j+1). Therefore, Cee yee 


ai Owe yt be Pein epee beac 


S 
in the previous expressions can be replaced by the notations 
Bee) py ee et, 3), fy (1,5)5 BUBAIST/ 61,971), 
B(i,j) and §(i,j), respectively. 

By introducing all the foregoing flux terms into 


Equation 7.1, one obtains the following explicit equation at 


each interior node 
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The coefficients of Equation 7.15’ are defined by 
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As mentioned above, the sign of both elements in the 
maoouct sea tin) J eUai-1/72,5), ait, 3) Ulit1/2,9), 611,39) 
V(i,j-1/2), and B(i,j+1) V(i,j+1/2) must be the same. 
Therefore the products are all positive quantities and can 
be written as noted in Equation 7.16 with the understanding 
Enat all a'S and —’s used in this work are positive number, 
The finite-difference approximations of the 

vorticity-stream function relation are derived from the 
central difference approximations for non-uniform grid 


network, where 
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myerntroducing Equations 7.18 and 7.19 into Equation 7.2, 
one obtains the following explicit equation for the stream 


function at each interior-node 
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The coefficients of Equation 7.20 are defined by 
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The distributions of vorticity, stream-function and 
concentration inside the region of interest can then be 
obtained by solving Equations 7.15 and 7.20 using an 
iterative method. The values of the parameters a, 8, y and 6 


can be specified depending on the differencing scheme used. 
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These values will be given later in this chapter. Procedure 

of the iterative method is discussed in Chapter 8. Since the 
variables ©, YW and C are evaluated at nodes, for consistency 
in analysing the numerical results, the velocity components 

U and V are also evaluated at the nodes by using the 


following equations: 
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The velocity components obtained from Equations 7.22 and 
7.23 are used instead of those evaluated from Equations 7.4 


ana 7.5 during the iterations. 


7.2 BOUNDARY CONDITIONS 


The finite-difference equations at the boundaries will 
be discussed for the five different regions mentioned above 


separately. 
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This boundary is defined as 1<i<nj and j=1 in Figure 


7.2. For an initial parabolic velocity profile, 
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The boundary condition for the concentration for both 


parabolic and flat velocity profiles are 
Cid 0 Craze} 


7.2.2 CONFINEMENT PLATE 
This boundary is defined as njsi<nx and j=1 in Figure 
7.2. The boundary condition for the stream-function at the 


confinement plate is 
a7 15) = =0%5 (727 ) 


The boundary condition of the vorticity at the 
confinement plate is approximated by a finite difference 
expression with truncation error less than 0(AY)? which 
embodies the no-slip conditions, 3$/aY=U=0 and 9?7W/av*=aU/ayY 


=a(i,1) at the boundary. The finite-difference equation at 
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this boundary is derived using Taylor series as 
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The boundary condition for the concentration at the 
confinement plate is approximated by using a forward 
finite-difference expression with truncation error less than 
O(AY)?.: For 0C/ay=0 at the boundary, the finite-difference 


equation becomes: 
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7.2.3 IMPINGEMENT PLATE 
This boundary is defined as 1<1i<nx and j=ny in Figure 


7.2. The boundary condition of stream-function.is given by 
W(i,ny) = 0 (730) 
The boundary condition of the vorticity is Somilar to 


that at the confinement plate. The finite-difference 


equation for vorticity at this boundary is derived using 


Hh! my 
i 


f. Shee a 


a AN a 


i 
edt 26 series fyeone?. oi See: t wok rigees ae 
. My aie a4 
Go awed of pitied a. betainiwovag a svet@, 


ses. 2S f vera Wa rere hind reed iketaienagne * 
po5H ‘x theh-aicegs eae ait vetbabie. ony: ae Esa 3 


“4 Es ®t. ‘s ped iol “(he « x! D oe wy cl ce ct jot, f * a 7 
ee eer Sipps Nina i aieialnandieuan ~ me + ere a ; hig (eas - a \ 


Hea a) +» ; a 


empet ab gHet ‘oe wamriat ae Hontiet: et. rabid | 2 
td navie at: eo i sat 2 wnatbinter vated et 
ey” | 


107 


Taylor series as 
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The boundary condition of the concentration is given by 
Civ ny) = 0 | (7232) 
7.2.4 AXIS OF SYMMETRY 
This boundary is defined as i=1 and 1<j<ny in Figure 


7.2. The boundary conditions of stream-function and 


vorticity are given by 
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For the boundary condition of the concentration, a 
boundary control volume surrounding a boundary node (1,j) is 
introduced. Referring to Figure 7.2, the node (i,j) becomes 
(1,3) and 4X(1)=4x(0)=4X(1). The index i21 will refer to the 
variables inside the system, and i<1 will refer to the 


variables outside the system. Equation 7.15 15 now applied 
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to the nodes on this boundary. Variables such as v(1,j+1/2), 
Bip jel). 641,341) sand. (1.,4+1) can be evaluated at the 
boundary. Due to the fact that this boundary can be 
considered as a plane of symmetry, the variables outside the 


field can be evaluated as follows 


UO ee a C25.5,) 
leo eae kth 3,20,.%) 


ClOr gece 3) (7.34) 


mer retore thesconcentration at »this,boundary,»C(1;3) «can be 
evaluated from Equation 7.15 similarly to that of the 


interior node. 


7.2.5 OUTFLOW REGION 
This boundary is defined as i=nx and 1<j<ny in Figure 
7.2. The finite-difference equations for fully developed 


Sstream-function, vorticity and concentration are given by 
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The finite-difference equations for the stream-function 
and vorticity at this boundary can be obtained from 
Equations 6.25 and 6.26. After rearrangement, they become, 


respectively 
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(nx, 5) = OCG) + [Wlnx-1,5) - O05) 2, (7.36) 


A(nx,j) = 24,05) + [2(nx-1,5) - o¢,(5)] 2, (7.37) 


with 2, = exp {-8 a,* [4x(nx-1)] / (L Re,)} (7.38) 

The finite-difference equation for the concentration at 
this boundary can be obtained from Equation 6.29 by 
approximating the ratio of the series for X=X(nx) and 
| X=X(nx-1) as the ratio of is first term of these series, 


where 
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minh, 2), cexprlermieds Amish) weel2r L Ren Sc) } (7.40) 


7.3 FINITE-DIFFERENCING SCHEMES 


Three finite-differencing schemes are used in this 
work, namely, a central differencing scheme (C.D.S.), an 
upstream differencing scheme (U.D.S.) and an 


upstream-weighted differencing scheme (U.W.D.S.). 


7.3.1 CENTRAL DIFFERENCING SCHEME (C.D.S.) 

Setting a(i,3)=611,7)=7(1,3)=6(7,3)=0 im the 
finite-difference equations of Section 7.1 results ina 
central differencing scheme with truncation error 
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7.3.2 UPSTREAM DIFFERENCING SCHEME (U.D.S.) 

Ssettingsatr, |i=pur,,)=0.5 and -y(i.))=6(1)4) 20 on: the 
Pate difference equations in Section 7.1 results iman 
upstream differencing scheme. This differencing scheme was 
used by Torrance (94) and Raithby and Torrance (69), and is 
Similar to the "upwind differences" used by Runchal et. al. 
(74,75,76). In this case, the variables 2, WU and C in the 
System are determined to a large extent by the values of the 


corresponding variables prevailing immediately upstream. 


7.3.3 UPSTREAM-WEIGHTED DIFFERENCING SCHEME (U.W.D.S.) 

Im this case, the parameters a(i,j), Bli,a);-y(2, 3 ane 
6(i,j) are determined from the exact solution to the local 
one-dimensional transport equations between nodes. Consider 
the transport equation of the variable @ in the X-direction 
between nodes (i-1,j) and (i,j), the one-dimensional 
equation governing the variation of the variable ) with X 


between nodes (i-1,j3) and (i,j) is 


dd da? 
Ua 2) ema ee, 
ax ax? 
The exact solution of Equation 7.41 is given by 
OO) exp (01 1/2,.5)) (a eit) an 
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Locally exact solutions of the form of Equation 7.42 were 
introduced by Spalding (88) and Raithby and Torrance (69). 
Mie) vakue iof (deat. Cis025)9), where X=xX(i-1)+4(Ax(i-1)72), can 
be evaluated from Equation 7.42. This "true" value of the 
Varrable:Qers equal ctow(i-1/2-a(i,j),j) of Equation 7.8. 


Meeterore, Eremebquations 7.8 and 7.42, with e,=a(i',7). one 


W 
obtains 
1 exp (Rw/2) - 1 
a(i,j) = - = ———_—_ (7.43) 
2 exp (Rw) - 1 
wpere Rw = |U(i>1/2Re) pw xii) / a (7.44) 
Similarly, 
1 exp (Rs/2) - 1 
Tan Cy EA Gis 8 0 renee aera rene ne (7.45) 
2 exp (Rs) - 1 
wher elaRse= BhV,( in 1/72 bieaejo ws 7-2 (7.46) 


The absolute signs are introduced since a's and £'S are 
always positive numbers, and for Se abn lesey reasons are 
chosen to match the signs of local mean velocities. From 
Equation 7.13, the first derivative of > in R=direct10n at 


(i-1/2) is given by 


a fotig)) = Oi ie) 
dx|i-1/2 Axi 1) 
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Since by introducing the parameter y(i,j), the 
finite-difference equation should yield the same result as 
the locally exact solution. The right hand side of Equation 
4.47 cam be evaluated. from: Equation 7.42. After 


rearrangement, one obtains 


Rw exp (Rw/2) 


Sa)! EP) 0 ee cern cr oe (7.48) 
exp. (Rw) =. 1 ) 
Similarly, 
Rs exp (Rs/2) 
Octet 0) eign pe (7.49) 


expan Rode =” 1 


Clearly, in the limit of Rw approaching zero both a(i,j) and 
y(i,j) approach zero and the finite-difference equations 
reduce to central differencing scheme. On the other hand, 
for large Rw, a(i,j) approaches 0.5 and y(i,j) approaches 1. 
This limit implies that upstream differencing scheme is used 
for the convection term and that the diffusion term is 
dropped. Similarly, B(i,j) and 6(i,j) approach the same 


value with respect to the value of Rs. 
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7.4 STABILITY OF THE FINITE-DIFFERENCE EQUATIONS 


The finite-difference equations derived in Section 7.1 
are a set of non-linear algebraic equations to be solved 
iteratively. 

For a set of linear algebraic equations with constant 
coefficients such as 


K(i) Pah) ete le (7.50) 


BH 


The matrix theory states that such a set will converge toa 
solution in @ successive substitution method when the matrix 
i.) ) 15; diagonally dominant". (45974). .This condition can 
be-expressed as d ian vis at Oree lie pandelorgat (Lease 
all i 
one i such that the inequality holds. 
For non-linear equations, the above conditions are 

often sufficient, although they may not be the necessary 


Sond bt onsen 627) aethnerefore, inworderwiorsEquation 7.15 to 


Satisfy the conditions for convergence, the following 


ol G7 25 


must hold for all the interior nodes and with strict 
inequatity for at least one node. 
The stability of the finite-difference equations for 


the three finite differencing schemes will be discussed 
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separately. 


7.4.1 CENTRAL DIFFERENCING SCHEME (C.D.S.) 

For the central differencing scheme, «&(2,3)=B(1, 79 
=y(i,j)=6(i,j)=0, one can immediately observe from Equation 
fervor that convergence of Equation 7.15 is not always 
assured: the sum of the terms on left hand side of Equation 
we pecause Of the presence ‘Of the velocity terms in the 
numerator, is not bounded. As long as all the coefficients 
in Equation 7.16 are positive, the sum of these terms is 
unity. However, once this is not true, the sum of these 
terms may well exceed unity. In the region near the jet 
MozzZie Exit, the velocity in the Y-direction, V, is large. 
For the case of high Reynolds numbers, the coefficient ms, in 
Equation 7.16 becomes negative. Such a coefficient may lead 


to non-convergence or numerical instabilities. 


7.4.2 UPSTREAM DIFFERENCING SCHEME (U.D.S.) 

For the upstream differencing scheme, a(i,j)=B8(i,j)=0.5 
and y(i,j)=6(i,j)=0, one can immediately observe from 
Bauation #216 that all the coefficients are always positive 
and hence the sum of the terms on the left hand side of 
Equation 7.51 will always be unity. Therefore, the upstream 


differencing scheme has a better chance of convergence. 
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7.4.3 UPSTREAM-WEIGHTED DIFFERENCING SCHEME (U.W.D.S.) 

For the upstream-weighted differencing scheme, a(i,j), 
B(i,j), y(i,j) and 6(i,j) are unknowns and are depended on 
the local velocity components and the grid size by studying 
Bquations 4.43, 4.45, 4.48 and 4.49, The limit values of 
Pnese .parameterus jwere shown.in Section 7:3 where o( ij; j)jwand 
Boe avary from 0 to.0.5. and +(i,35) and 46i,j3)..vary.from 0 
to 1 with Rw and Rs increasing. 

In order to have a better chance of convergence, all 
the coefficients in Equation 7.16 should be positive such 
that the sum of the terms on the left hand side of Equation 
7.51 1S unity. Since a negative coefficient occurs only when 
the local velocities are large, let ‘us examine the 
coefficients for such cases. In the region of interest, the 
largest U occurs near the stagnation flow region and the 
Margest V occurs near. the jet nozzle exit. Therefore the 
warues Of o(i,7) and y(1,)) are close to 0.5 and 1 néar the 
Stagnation flow region, respectively. On the other hand, the 
evalues of B(1;3) and 6(i,j) are close to 0.5 and 1 near the 
jet nozzle exit, respectively. From Equation 7.16, all the 
coefficients should be positive for both cases with m, 
approaching zero for the first case and m, approaching zero 
for the second case. 

Once again, all the coefficients in Equation 7.16 are 
always positive. The convergence properties of this 
differencing scheme should be similar to those for the 


upstream differencing scheme. 
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8, COMPUTATIONAL PROCEDURE 


The grid design used for the numerical computations is 
discussed in Section 8.1, followed by a detailed 
presentation of the method of solution in Section 8.2. The 
selection of the convergence criteria is discussed in 
seeru20n 8.3. Finally, the construction of the computer 


program is given in Section 8.4. 


8.1 GRID DESIGN 


A non-uniform grid is used in the numerical 
computations. The grid size and the arrangement are known to 
be important factors in determining not only the accuracy of 
the solution, but also the convergence characteristics. 

Due to relatively large gradients of velocity and 
concentration along the impingement plate, the gridlines 
parallel and adjacent to this plate must be very closely 
Spaced. On the other hand, gridlines parallel and adjacent 
to the axis of symmetry are also closely spaced in order to 
ensure accurate calculation of the variables within the 
Stagnation flow region. 

In this work, different grid networks are designed for 
the numerical runs for the three different jet-to-plate 


spacings (L=2,4 and 12). These grid networks are results of 
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numerous trials and are given in Table 8.1 Tests for the 


accuracy of these grid networks will be given in Chapter 38. 


TABLE 8.1 : GRID NETWORKS FOR DIFFERENT L AND Re 


L Re, nx ny 
Z 100 - 400 role) 20 
4 100 - 400 25 25 
12 10072-3000 67 25 

400 69 25 


8.1.1 GRID ARRANGEMENT IN X-DIRECTION 

The finest grid spacings are adjacent to the axis of 
Symmetry. In such a region, nine nodes are used to cover a 
distance of b/2 (nj=9). The grid spacings are then increased 
Mnestens byadetactor of 2 in the x-direction. The Largest 
grid spacings appear in the outflow region which are 128 
times the finest one. The grid arrangement in X-direction 
are listed in Table 8.2. 

The distance in X-direction which is covered by grid 
network is a function of nx, the total number of nodes in 
X-direction used. In other words, the location of the 
outflow region boundary is a function of nx. For the grid 
networks with nx=55, 67 and 69 as shown in Table 8.1, the 


outflow region boundaries are located at X=74, 170 and 186, 
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TABLE 8.2 
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respectively (see Table 8.2). These locations of the outflow 


region boundaries are also mentioned above in Table 6.1. 


8.1.2 GRID ARRANGEMENT IN Y-DIRECTION 

Since finer grid spacings must be iocated adjacent to 
the impingement plate, this can be achieved with a 
coordinate stretching transformation of some kind. One of 
the most commonly used transformation is that of an 


exponential stretch. Choose 


Sige ween by ayy] 
Y= LL fesiga hy) 
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Such that the Y'-coordinate is transformed into the 
Benetched Yocoordinate..ib,. 1S an arbitrary constant, used to 
adjust the "stretch". Equation 8.1 satifies the boundary 
Senaitrons .ob iyi )=0, W'rGl)=0 and Y.ny)=L, (ny) =L. 

Pn Cn iS work:  D,S. ChOoOSen asS).0./757, Gaco and: 03075 stor 
numerical runs of L=2, 4 and 12, respectively. Numerical 
Huns With b,=1..0) for L=2 and b,=0.75 for Lad iare also made 
for the case of an initial parabolic velocity profile. The 
grid arrangements in Y-direction for three different 
jet-to-plate spacings are listed in Tables 8.3, 8.4 and 8.5 
for L=2, 4 and 12, respectively. A total of 25 nodes are 
used in Y-direction (ny=25) throughtout all the numerical 
runs. Sample runs with ny=17 for L=2 and 4, and also with 


ny=33 for L=12 are only used for the purpose of testing the 
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TABLE 8.3 : 


Dp, <€ 0.75 
3 Y 
1 0 
2 0 ah 
3 0.3038 
g 0.440 
5 0.569 
6 0.6911 
7 0:338.05 
8 0942 
2 1 O18 
10 1208 
1 1.9196 
12 A280 
He 1.856 
14 1.432 
ips) 15041 
16 12566 
17 it 627 
18 1.665 
19 1 39 
20 1/689 
Zi 1.1837 
Page 1 682 
ae) 1.924 
24 1205 
25 2.000 


GRID ARRANGEMENT 
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SPACING 
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IN Y-DIRECTION FOR L=2 


5, = 4 
5 ¥ 
1 0 
2 O85 
3 i359 
4 OSSuZ 
S 0.656 
6 0.788 
“ 0.970 
8 1.022 
9 deals 5 
10 Pee 
oa 1.308 
12 +.888 
gs: 1.462 
14 } oe 
15 GS abe! 
16 L600 
17 1.03 
18 toa 
19 ioe 
20 1.6358 
2a 1.67.6 
22 Iv ocoeial 
23 1 Oa 
24 NP LIES: 
ao 2.000 
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TABLE 8.4 : GRID ARRANGEMENT IN Y-DIRECTION FOR L=4 


b, e025 b, = 0.75 
RELATIVE RELATIVE 
: a SPACING 5 M4 SPACING 
1 0 2.606 1 0 1.7.9 
2 02258 2.506 Z 0.495 VO Ota 
3 0.506 2.404 3 roo 3m W750 
4 0.744 2.293 4 lv, Seb 0 Tees 
5 Oe oa 2.202 5 7eOO6 TOVS ite 
6 1.190 Ze tz 6 1.956 9.464 
y 1.400 2.030 7 ree 2a soo" 
8 1,601 1.949 8 Zevos (pera, 
9 1.794 1.869 9 2.064 6.500 
10 1997-9 ic. 97 GO 10 2.843 Saou e 
ie FBGA Fo: nt ad: Le 3000 Sens 
12 220 1.646 TZ B7145 4.464 
13 2.490 Te Ors les yeray ad) 3.964 
14 2.646 5 25 14 Seer 3.464 
(3 27 Te4a55 15 3.478 SOs 
16 2.941 1.394 16 a, 064 2.714 
17 3.079 1.343 le? 3.640 2.393 
18 Sere 1.209 18 Oye OF 2.107 
19 IS he, e232 19 3. /06 Toe 
PY, 3.461 et PSi2 20 CAs ils: 1.643 
Za Rua his. Le. 134 Ze 3.864 1.464 
ZZ 32690 16091 Zz 8905 kaZ2ou 
Zs 3./98 1.040 23 3.3940 io Las 
24 3.904 1.000 24 36972 1.000 
wha! 4.0060 25 4.000 
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TABLE 8.5 : GRID ARRANGEMENT IN Y-DIRECTION FOR L=12 


bye 2002075 
RELATIVE 
j Y SPACING 
1 0 2.369 
2 0.744 oa 208 
3 1.461 2.201 
4 2hu52 on 
5 2817 2.038 
6 3.457 1.965 
7 4.074 1.895 
8 4.669 1822 
g 5.241 1.755 
10 5.792 1.691 
te 6.323 1.631 
12 6.835 1.570 
13 7.328 2540 
14 7.802 1.455 
15 8.259 1.405 
16 8.700 1.350 
ig 9.124 1.299 
i8 9.532 1.255 
19 9.926 1.204 
20 10.304 1.162 
2) 10.669 Vertex 
22 dere 0.20 1.080 
23 11.360 1.038 
24 11.686 1.000 
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accuracy of the numerical solutions using a particular grid 


network. 


8.2 METHOD OF SOLUTION 


The finite-difference equations used to discretize the 
governing equations and the associated boundary conditions 
shown in Chapter 7 constitute a system of strongly 
non-linear algebraic equations. For a nx x ny grid network, 
there are (nx-2) x (ny-2) algebraic equations for each of 
Enenthree variables, WU, 2 and C making a total of 
3 x (nx-2) x (ny-2) equations to be solved. Because of the 
difficulty in solving a strongly non-linear system of 
equations, a linear algorithm is used such that the 
non-linear coefficients are updated periodically. The 
Gauss-Seidel iteration method coupled with the conventional 
SuccesSive over relaxation (SOR) method is used. This method 
has been established that an optimum relaxation factor 
exists which yields the maximum rate of convergence for 
linear problems. However, for a non-linear problem the 
relaxation factor varies not only from node to node but also 
with every iteration. An estimation of this relaxation 
factor is not worthwhile from a computational efficiency 
point of view. In such cases a constant value is used which 
is chosen based on numerical experimentation. In this work, 


the proper relaxation factors for the stream-function and 
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the concentration are 1.7 and 0.7, respectively. While for 
enervorticity, the relaxation factor is in. the range of 0.1 
to 0.4. 

A modified scheme of calculations, presented by Wilkes 
(100) is used during the iteration procedure. This modified 
scheme has been observed by MasSliyah and Nandakumar (49) to 
have a stabilizing effect on the convergence characteristics 
and it is twice as fast a the conventional SOR method. The 
depiction of this scheme is shown in Figure 8.1. The grid 
network is divided into two subgrids. The circles refer to 
Subgrid 1 and the crosses refer to subgrid 2. The iteration 
is carried out in two steps, one subgrid being covered first 
ance the other next. Inthe first half of iteration, the 
nodes of the subgrid 1 are computed using the neighbouring 
four nodes of the subgrid 2 as given by Equations 7.15 and 


7.20, together with 


1 


ac Nae ve ‘ Cit (2S) a ee 
a UGG gee ay (ieee [Si (Given See Cet 


Coin} 


where S'')(i.j) is the value of S(i,j) (represented 
variables 2, J and C) at the nth iteration belonging to 
Subarid 1; Sik katy It) is the value of S$(i,j) computed using 
the four neighbouring nodes of subgrid 2 and w is a constant 
relaxation factor. Therefore, all the nodes in subgrid 1 are 
computed first and only points belonging to subgrid 2 are 


used. 
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fe) subgrid 1 


X subgrid 2 


HEME 
FIGURE 8.1 : SUBGRIDS USED FOR THE ITERATIVE SC 
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In the next half of iteration, nodes of subgrid 2 are 
evaluated using the neighbouring four nodes which belong to 
Subgrid i. The values of these four nodes in subgrid 1 
computed in the first half iteration are used. In this 
second half iteration, the same relaxation factor, w, is 
used. This completes one full iteration. 

BaNeeronser. 47 ,eo, le lorand P (6vare solved 
Simultaneously to obtain ee velocity, vortreity and 
Stream-function distributions. The summary of the procedure 
used in.solving the stream-function and vorticity transport 
equations is listed below 
1. For a low Reynolds number case (Re,=1), arbitrary 

constant values are used for the starting values of 
stream-function and vorticity. For the cases of higher 
Reynolds number the converged solution obtained for a 
lower Reynolds number is used as the initial starting 
values. 

2. “Any Variable “which fs a*constant “1s *computed=berore 
entering the iteration loop. 

So SWrehin ebhebreeracion Sloop oO 1/72 )4) ‘and V(1,j3+1/2) are 
computed from Equations 7.4 and 7.5 using the guessed 
values of the stream-function. 

4. The parameters a, B, y and 6 for upstream-weighted 
differencing scheme (U.W.D.S.) are then evaluated. One 
may notice that these parameters for central 
differencing scheme (C.D.S.) and upstream differencing 


scheme (U.D.S.) are constants and can be evaluated 
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before entering the iteration loop. Even for the case of 
U.W.D.S., the values of these parameters do not need 
frequent updating (69). In this work, these parameters 
are updated after every 50 iterations. 

thePnonriineare coefficients m,° mz; my, myvand’'m, are 
updated using Equation 7.16. Then the vorticity at each 
interior node is iterated in the manner mentioned above 
uSing Equation 7.15. | 

Next, the stream-function is iterated in the same manner 
atoadi interior nodes using Equation 7.20. This 
iteration is carried out three times in order to get a 
smooth solution in the manner suggested by Masliyah and 
Nandakumar (49). 

The vorticity boundary conditions at both the 
impingement and confinement plates are updated. The 
vorticity and stream-function boundary conditions at the 
outflow region are also updated. 

The modified variables are treated aS improved guesses 
and steps 3 to 7 are then repeated until the process 
Satisfied specified convergence criteria. 


Oneecyole nforesteps u3ethnoughe7is referéd to®beimng one 


iteration. The convergence criterion is tested at the end of 


each iteration. 


It is found that a relatively large number of 


iterations is required for a converged solution. For 


example, for Re = 100 and L=4, 1353 iterations are needed to 


obtain the converged solution. In this case, the converged 
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solution of Re,=1 and L=4 is used as the initial guess. The 
CPU time used for this run is approximately 110 s on Amdahl 
470/V8 computer. 

Atter the converged solutions of vorticity and 
stream-function are obtained, the velocity components are 
evaluated at each interior node using Equations 7.22 and 
7.23. In addition, an important variable, the local skin 
meretion factor, is evaluated along the impingement plate 
mrom the flow field. The local skin friction factor, Ce, is 


defined as 


me = Y? 
Cya= ees COO Vv ) (83) 
where p is the density of air and t, 1s the shear stress at 
the impingement plate. The shear stress at the impingement 


plate (y=h) is given by 


tT =- (— +-) (8.4) 


Since the impingement plate is impermeable, implies that 
v=(av/ax)=0 at the impingement plate. Introducing the 
dimensionless variables of Equation 6.6 and combining 


Equations 8.3 and 8.4, yields 
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Cpeeneic be (8.5) 
Rep Y=L 


By uSing Equation Gyo, the local Skin friction factor can 
then be evaluated from the vorticity at the impingement 
plate and the jet Reynolds number. 

The transport equation of the swelling agent is solved 
tO Obtained the concentration distribultion. The procedure 
is mainly the same as mentioned above ieee in this case 
the velocity components are known, therefore the 
Scerricients m,, M,, M,, m, and m, do not vary with 
iteration. In other words, the system of algebraic equations 
to be solved is linear. 

After the converged solution of concentration is 
obtained, the bulk flow concentration is evaluated. The bulk 


flow concentration, c is defined as 


B f 
h 
Ceate cf unc ay) <7 (Wo) h) (8.6) 
° 


where U, is the average velocity in the outflow region and 
is equal to 0.5v,b/h for the jet system in this work. 
Introducing the dimensionless variables of Equation 6.6, 


yields 


; | 
Cc, = 2 f Use zdy (o-7) 
Oo 


In addition, the local Sherwood number is also evaluated 
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along the impingement plate from the concentration field. 


The definitions of the Sherwood number are 


given in 


Equations 5.1 and 5.2. Rewriting the definitions of the 


Moca. Mass transfer coefficients k and k' 


in Bauat vonse oe 


and 5.4 using ideal gas law and Dalton's law, yields 


Ne =k = Be Ls 
ath, cj) Mw 


and N ea ee 


i cog tee 
According to Fick's law, the mass transfer 


from the impingement plate is described by 


Introducing the dimensionless variables of 
combining Equations 8.8 and 8.10, and also 


S.10, yield 


dC 
Sh, = -(—) 
2 oY |Y=L 
3C 
Vie ale! yall Voom oly 
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(8.8) 
ae) 


of swelling agent 


GB .gi0) 


Equatvon .6.67and 


Equations 8.9 and 


COmalia: 


By using Equations 8.11 and 8.12, the local Sherwood numbers 


can then be evaluated from the gradient of concentration at 
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the impingement plate. 


8.3 CONVERGENCE CRITERION 


Convergence of the numerical results is assumed when 
the maximum absolute difference between two consecutive 


iterations, defined as 


éu= ee (SoA ea Bae Sanne lore CS) 
@eeless than 10° * for vorticity and stream-function, and 
10°* for concentration. The suitability of these convergence 
Criteria is confirmed by numerical experimentation. Once the 
convergence criterion is met, further iterations have no 
Significant effect on the local skin friction factor and 


local Sherwood number along the impingement plate. 


8.4 OUTLINE OF THE COMPUTER PROGRAM 


A main program with four subroutines is used in the 
Numerical study. The first subroutine, ITER1, is used to 
compute the stream-function and the vorticity distributions. 
After the converged solutions of the stream-function and the 
vorticity are obtained from ITER1, the second subroutine, 
CALCT, 1s used to compute the velocity components at all the 


nodes instead of midway points. In this subroutine, the 
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local skin friction factor along the impingement plate 6 
also evaluated. The third subroutine, ITER2, is used to 
compute the concentration distribution. After the converged 
Sorutronrof=the=concentration is obtained’ from ?PTER2-~the 
fourth subroutine, CALC2, is used to compute the bulk flow 
concentration and the local Sherwood number along the 
impingement plate. 

Only the computational procedure of subroutine ITER1.is 
shown in Figure 8.2 in the form of a flow diagram. The 
computational procedure of subroutine ITER2 is very similar 
to that of subroutine ITER1. The computational procedures of 
subroutines CALC1 and CALC2 are straightforward calculations 
only. 

The program is listed in Appendix F, together with a 


myprcal output listing. 


it i Trt fi > a i 4 


nat 
ia Yee 
_ Ss ea ise re sits 
i sani heaoe bea fi fig 69 ay * 98 be nee 


a ; et a4 ‘phat m) cr ws - ; 
L why 2 a"R Galena’: ere aa wate 8.5 . 
, ho he 3) Saw fat ge ew op Sf 4 ie . 
: | 7 6 
: ‘ : y : f is r Wa . ie _ 7 
ie oe wate tO. 
ee! j om 
fea 
i 
vey : : 
=p TTA i 
; a aed 
. | Th = - = 
r - > 2 > un ’ r) 4 oe ae 
i 
i r= . * ey 
eee hie 
i ) 
ail * 
" : : r 
hs ate a 
£ 7 : at 
¢ a . 
'd ” 
f 
=“ 
t fe rs 
bm ‘ 
é 
uN , 
ri . : 
> 1 
” 
al 
are 
a 
\ 
{ i 
: i 
a aed > Ae ee OS -— 
“ a at 7 
2 ss * a r os ar + an q = 
7 7 ie ean 
oT Po. : al 
y 
\ 


START 


READ PROBLEM PARAMETERS: Re, AND L 
READ INFORMATIONS FOR GRID DESIGN: nx, ny AND b, | 
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CHOOSE DIFFERENCING SCHEME 


PCIE iN Naty Ges a, "Sete id 


OP Wh — 


1. SET-UP GRID ARRANGEMENT 
2. COMPUTE PROGRAM CONSTANTS 
3. COMPUTE BOUNDARY CONDITIONS } 


- 
| COMPUTE U AND V USING EQUATIONS 7.4 AND 7.5 | 


Se 
{eo Rates son Che ncaa 


ITERATE ON 2 AND JU USING EQUATIONS 7.15 AND 7.20 | 
| UPDATE BOUNDARY CONDITIONS AND VELOCITY coMPONENTS | 


tb AND 2 
CONVERGED ? 


STORE RESULTS 


FIGURE 8.2 : COMPUTATIONAL FLOW DIAGRAM FOR 
SUBROUTINE ITER1 
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9. VALIDITY OF THE NUMERICAL SOLUTION 


Accuracy of numerical solution is dependent on the 
choice of grid network and differencing scheme. The 
influences of grid network and differencing scheme on the 
numerical solution are studied separately in Sections 9.1 


and 9.2, respectively. 


9.1 INFLUENCE OF GRID NETWORK 


The influences of the grid network on the numerical 
solutions for three different jet-to-plate spacings (L=2, 4 
and 12) are studied separately. 

For L=2, numerical runs for different Reynolds number 
Mcind two Giffierent grid networkseof $55 x 17 and 55 x «25 
with b,=0.75 are made. The differencing scheme used is 
UeW Dwse. he Skinatfrictiom factor evaluated along the 
impingement plate from these two grid networks are plotted 
in Figure 9.1 for Re,=100, 200, 309 and 400 with an initial 
Parabolic velocity profile. There is little ‘difference 
between the values of C,Re, obtained from the two different 
grid networks for Re,,= 100 and 200. For Re,=300 and 400, 
disagreement between the values of C,Re, obtained from these 


two grid networks is found only in a small region near the 


Stagnation flow region with a maximum difference of 
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FIGURE 9.1 3: 


INFLUENCE OF GRID NETWORK ON THE NUMERICAL 
SOLUTIONS FOR L=2 USING U.W.D.S. 
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approximately 5% only. Therefore, a conclusion can be made 
that both grid networks are suitable for the computation for 
the case of L=2. In all subsequent computations for the case 
Of L-2 the network of 55. 4° 25 with the smaller Grid’ size is 
used. 

For L=4, numerical runs for different Reynolds number 
mSang uwo grid’ networks of 55 x 17 and 55 # 25 with b,=0925 
are made. Again, the Panetencii scheme used 1S U.W.D.S.. 
The skin friction factors evaluated along the impingement 
Plate from these two networks are plotted in Figure 9.2 for 
Re,=100, 200, 300 and 400ewith.an initial. parabolic, velocity 
profile. There is little difference between the values of 
C, Re, obtained from these two grid networks for Re,=100. But 
for the Re,=200, 300 and 400, disagreement is found mainly 
in the stagnation flow region with the worst case occured 
for Re,=400. The maximum differences are 6%, 15% and 20% for 
Re,,=200, 300 and 400, respectively. Therefore, the two grid 
network are only suitable when Re,<200. But the network of 
55 xSl7 with the coarser geid size™ius no/donger suitablestor 
high Reynolds numbers. In all subsequent computations for 
Pheucasesor h<-4 the network of 554.25 with the finer “grid 
size 1s used. 

From the conclusion above, it is obvious that a grid 
network of 55 x 17 will be too coarse for the case of L=12 
because of the large jet-to-plate spacing. For L=12, 
numerical runs for Re ,=100, 200 and 300 using two different 


grid networks of 67 x 25 and 67 x 33, and for Re,=400 using 
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FIGURE 9.2 : INFLUENCE OF GRID NETWORK ON THE NUMERICAL 
SOLUTIONS FOR L=4 USING U.W.D.S. 
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two grid networks of 69 x 25 and 69 x 33 are made. For all 
these grid networks, the value of b, is chosen to be 0.075. 
The differencing scheme used is U.W.D.S.. Unfortunately, for 
Re ,,=400, no converged solution can be obtained using the 
grid network of 69 x 33. The skin friction factors evaluated 
along the impingement plate from the two grid networks of 
67 x 25 and 67 x 33 are plotted in Figure 923 for Re =100, 
200 and 300 with anyinitial) parabolic profile. There is 
again little ditference between the values of C,Re, obtained 
from these two grid networks for Re,=100. For the cases of 
Re ,=200 and 300, disagreement is found mainly in the 
Stagnation flow region. In all subsequent computations for 
the case of L=12, the metworki of 67 % 25%6r 69 x 25 with the 
coarser grid size is used simply because a converged 
solution can be obtained from this network for all Reynolds 
numbers. It is noted that in order to obtain a converged 
Solution for the case fof L=12, a little sacrifice on ~the 
accuracy of the numerical solution in the stagnation flow 
region cannot be avoided. 

FOr a partictlarinumerical run iInmithis work, only one 
grid network is used. These grid networks for different 
jet-to-plate spacing and Reynolds number are listed in Table 
Bail. 

The influence of grid arrangement on the numerical 
solution is not ae obvious as that of the grid network. For 


b=2 numerical wsolvutions Of flow and concentration fields 


evaluated from grid network of 55 x 25 with b,=0.75 and 1 
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FIGURE 9.3 : INFLUENCE OF GRID NETWORK ON THE NUMERICAL 
SOLUTIONS FOR L=12 USING U.W.D.S. 
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are nearly the same for all Reynolds numbers. The variables 
such as local skin friction factor and local Sherwood number 
evaluated along the impingement plate from the two different 
grid arrangements are nearly the same for all Reynolds 
numbers. For L=4, numerical solutions of flow and 
concentration fields evaluated from grid network of 55 x 25 
with b,=0.25 and 0.75 are also nearly the same for all 
Reynolds numbers. Only in evaluating thet locadarskiny frtictaon 
factor and the local Sherwood number along the impingement 
plate, thus disagreement in these two arrangements is found 
near the stagnation flow region. In this work, the variables 
evaluated from the grid network of 55 x 25 with b,=0.75 are 
used because these variables should be more accurate due to 
the finer grid spacings adjacent to the impingement plate. 
For L=12, only one grid arrangement with b,=0.075 is used 
for both, the.grid network of 67 »x<25..and, the<gridinetwork of 


Sux. 25; 


9.2 INFLUENCE OF DIFFERENCING SCHEME 


The influence of differencing scheme on the numerical 
solution is studied briefly for the case of L=2. To conserve 
on the computational costs, a grid network of 55 x 17 with 
b,=0.75 is used. This grid network was found in Section 9.1 
to be suitable for the computations for the case of L=2. 


Numerical runs with an initial parabolic velocity profile 
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are made using three different differencing schemes (C.D.S., 
U.D.S. and U.W.D.S.). Converged solutions are obtained for 
@uiertns Using U.DiS. and ULW.D.S., while the converged 
solutions can only be obtained for Re, <200 using C.D.S.. 

The comparisons of the converged local skin friction 
factors obtained from the numerical runs using the three 
differencing schemes are shown in Figure 9.4 for Re,= 100, 
200, 300 and 400. Iteas ein that the differencing scheme 
has little effect on the solutions except in a small region 
near the stagnation flow region. In the range of 
disagreement, some of the observations by Raithby and 
Torrance (69) on comparison of differencing schemes are 
confirmed. According to their study, solutions from the 
C.D.S. are the most accurate solutions compared to the exact 
solutions as long as converged solutions can be found. 
Although U.D.S. gives converged solutions for all cases, in 
the range where C.D.S. can be used, U.D.S. is inferior. 
Furthermore, U.W.D.S. also gives converged solutions for all 
cases and the solutions obtained from U.W.D.S. are better 
than those from U.D.S.. These observations are confirmed as 
shown in Figures 9.4(a) and 9.4(b) by assuming the solutions 
from the C.D.S. are the most accurate solutions in this 
work. The disagreements between the solutions of C.D.S. and 
U.D.S. are the largest in-both cases which indicate that the 
solutions of U.D.S. are worse than those of U.W.D.S.. 

Since the solutions of C.D.S. are the best, it is 


logical to use C.D.S. as long as it gives a converged 
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FIGURE 9.4 ; INFLUENCE OF DIFFERENCING SCHEME ON 
THE NUMERICAL SOLUTIONS FOR L=2 
USING A GRID NETWORK OF 55 x 17 
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solution. But unfortunately, the numerical solutions 
obtained by using C.D.S. for the chosen grid network of 

55 x 25 did not give a converged solution. Therefore, in all 
subsequent numerical computations, only U.W.D.S. and U.D.S. 
are used while U.W.D.S. is expected to give a better 
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10. RESULTS AND DISCUSSIONS 


Experimental and numerical results for a confined air 
jet are both discussed in this chapter. Results are given 
for two different categories: flow characteristics and mass 
transfer characteristics. Flow characteristics can only be 
Seucdied! numerical ly, ine this,work,..While mass transfer 
characteristics are studied both experimentally and 
Numerically. Experimental results in this case are used to 
verify the numerical predictions of the two-dimensional 


model. 


10.1 FLOW CHARACTERISTICS 


The flow behavior for different jet Reynolds numbers, 
jet-to-plate spacings and nozzle exit velocity profiles are 
Studied numerically. A listing of the numerical runs is 
given in Appendix G. In this section, only the numerical 
results using the upstream-weighted differencing scheme 


(U.W.D.S.) are presented. 


10.1.1 STREAMLINE CONTOURS 
The flow field is studied qualitatively by observing 
the streamline contours from the numerical solution. The 


contours of the stream-function for L=2, 4 and 12 with an 
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initial parabolic velocity profile are shown in Figures 
10.1, 10.2 and 10.3, respectively. The contours of the 
stream-function for L=4 with an initial flat velocity 
profile are shown in Figure 10.4. The jet nozzle exit is at 
the upper left hand corner with the main flow travelling 
from left to right. The upper horizontal streamline 
represents the confinement plate and the lower horizontal 
Streamline represents the impingement plate. 

For =the-case ofan initiel parabolic profile, —in 
general, a primary vortex rotating in the counter-clockwise 
direction is found near the confinement plate with its size 
increasing with jet Reynolds number. With the exception of 
the case of low Reynolds numbers for L=2, a second vortex 
rotating in clockwise direction is found near the 
impingement plate. This secondary vortex is much smaller 
than the first one in size and it is also weaker in terms of 
rotational Jintensity..eThe, variation.of the.location lof the 
vortex centre with the jet Reynolds number is shown in 
Figure 10.54 It can be observed that the centers of both the 
primary and secondary vortices move downstream with increase 
in Reynolds number. 

For the case of an initial flat profile, only the 
primary vortex is found near the confinement plate. This 
vortex behaves as the primary vortex for the case of an 
initial parabolic profile except it is smaller in size and 
also weaker in terms of rotational intensity. The variation 


of the location of the primary vortex centre with the jet 
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FIGURE 10.1 : CONTOURS OF STREAM-FUNCTION FOR L=2 WITH 
AN INITIAL PARABOLIC VELOCITY PROFILE 
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FIGURE 10.2 : CONTOURS OF STREAM-FUNCTION FOR L=4 WITH 
AN INITIAL PARABOLIC VELOCITY PROFILE 
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FIGURE 10.3 : CONTOURS OF STREAM-FUNCTION FOR L=12 WITH 
AN INITIAL PARABOLIC VELOCITY PROFILE 
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FIGURE 10.4 : CONTOURS OF STREAM-FUNCTION FOR L=4 WITH 
AN INITIAL FLAT VELOCITY PROFILE 
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FIGURE 10.5 
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NUMBER FOR THE CASE OF PARABOLIC PROFILE 
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Reynolds number is shown in Figure 10.6. Again, the primary 
vortex move downstream with increase in Reynolds number. 

A blow-up of the contours of stream-function for L=2, 4 
eed) 2 With an initial parabolic profile are shown in 
Pagures 10.7, 10.8 and 10.9, respectively, -while that for 
Bea With an=initial flat. profile is shown in Figure 10.10. 
mie -Celtre Or the jer nozzle exit 41s at the upper left hand 
corner with the main Sion SREReS le ge from left to night. The 
corresponding values of stream-function for each streamline 


numbered are listed as follows: 


ie2Ui = 0 

2. P = -0.09 

Ba, Us 20. 201 

4. W = -0.33 

5. W = -0.45 

6.50. = 50.5 

CN Sas (circulating flow) 
Site NGS al es (circulating flow) 


where streamline #1 represents the axis of symmetry and the 
impingement plate, and streamline #6 represents the 
confinement plate and the outermost free streamline of the 
Submerged jet. 

The spreading effect of the jet under the the influence 
of the confinement plate can be studied from the outermost 


free streamline (#6) in Figures 10.7, 10.8, 10. 9and 10-10% 
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FIGURE 10.6 : VARIATION OF PRIMARY VORTEX CENTRE WITH 
REYNOLDS NUMBER FOR THE CASE OF FLAT 
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FIGURE 10.7 : CONTOURS OF STREAM-FUNCTION FOR L=2 NEAR 
THE STAGNATION POINT WITH AN INITIAL 
PARABOLIC VELOCITY PROFILE 
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FIGURE 10.8 : CONTOURS OF STREAM-FUNCTION FOR L=4 NEAR 
THE STAGNATION POINT WITH AN INITIAL 
PARABOLIC VELOCITY PROFILE 
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FIGURE 10.9 


THE STAGNATION POINT WITH AN INITIAL 


PARABOLIC VELOCITY PROFILE 
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FIGURE 10.10 : CONTOURS OF STREAM-FUNCTION FOR L=4 NEAR 
THE STAGNATION POINT WITH AN INITIAL 
FLAT VELOCITY PROFILE 
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The jet contracts slightly below the nozzle exit for an 
initial parabolic profile, while for an initial fiat 
profile, the jet expands continuously. Similar spreading 
effects were observed by Van Heiningen et. al. (96) in their 
study of semi-confined air jet with different initial 
velocity profiles. The explanations of these spreading 
behaviors have already been given in Section 2.1.2. 

From the outermost eee streamlines in Figures 10.8 and 
16.10, it 1s found that the free streamline is significantly 
closer to the impingement plate for an initial parabolic 
profile. This is due to the higher momentum of the initial 
parabolic profile jet. Same observation has also been made 


by Van Heiningen et. al. (96). 


10.1.2 AXIAL VELOCITY PROFILE 

Typical axial velocity profiles at various positions in 
Y-direction for different jet Reynolds numbers for a 
confined jet are shown in Figures 10.11, LORe imme CenOats eer 
a2, 4 anduil2 Withean anltial parabolic profile, and for b=4 
With amagnitial flateprefilesin Figures 10.14. Ingallscases 
phe init@al velocrcy atwnozale exit 4(Y=-0)"are plotted. The 
axial velocities are normalized by the centerline axial 
velocity att@the Mnozzle exit, (Vi) | x=0- 

The spreading of the jet as it approachs the 
impingement plate (Y=L) is obvious for all cases. This 
spreading is more dramatic for the jet with an initial flat 


Moofile than that with an parabolic profile. 
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FIGURE 10.11 : AXIAL VELOCITY PROFILES FOR L=2 WITH AN 
INITIAL PARABOLIC VELOCITY PROFILE 
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FIGURE 10.12 : AXIAL VELOCITY PROFILES FOR L=4 WITH AN 
INITIAL PARABOLIC VELOCITY PROFILE 
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FIGURE 10.13 : AXIAL VELOCITY PROFILES FOR L=12 WITH AN 
INITIAL PARABOLIC VELOCITY PROFILE 
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FIGURE 10.14 : AXIAL VELOCITY PROFILES FOR L=4 WITH AN 
INITIAL FLAT VELOCITY PROFILE 
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The decay of the centerline axial velocity, V| x=; 
Originating from a parabolic velocity profile at nozzle exit 
for different jet Reynolds numbers is shown in Figures 
ier, “10.16 and 10.17 for L=2, 4 and 12, respectively. The 
centerline axial velocity increases slightly with axial 
distance from the nozzle exit due to the contraction of the 
jet, and then decreases slightly due to the spreading of the 
jet. Not until the jet flow is about one slot width away 
from the impingement plate, does the centerline axial 
velocity decrease rapidly to zero at the stagnation point. 
In other words, the centerline axial velocity is affected by 
the presence of the impingement plate at a distance only one 
slot width from the plate. 

The effect of Reynolds number on the decay of the 
eenterline axial velocity for different jet-to-plate 
spacings can also be studied from Figures 10.15, 10.16 and 
10.17. *For all.cases, thescenterline axial velocity decays 
less rapidly at higher Reynolds numbers. This is mostly due 
to the more penetration of the jet with higher Reynolds 
number into the surrounding fluid. 

In, the stagnation flow regivon, the centerline axial 
velocity is linearly proportional to the axial distance away 
from the stagnation point, (L-Y), as can be noted from 
Figures 10.15, 10.16 and 10.17. A similar observation has 
been made by Schlichting (78) for unconfined submerged jet. 


Introducing the dimensionless variables in Equation 6.6 into 


Equation 2.5, yields 
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FIGURE 10.15 : DECAY OF CENTERLINE AXIAL VELOCITY FOR 
L=2 WITH AN INITIAL PARABOLIC VELOCITY 
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FIGURE 10.16 : DECAY OF CENTERLINE AXIAL VELOCITY FOR 
L=4 WITH AN INITIAL PARABOLIC VELOCITY 
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FIGURE 10.17 : DECAY OF CENTERLINE AXIAL VELOCITY FOR 


L=12 WITH AN INITIAL PARABOLIC VELOCITY 
PROFILE 


(V/V) | x= a a 6 ame (Vi )ix=o] [oP TE Ona) 


The values of a, can be easily evaluated from the slopes of 
all the curves of Figures 10.15, 10.16 and 10.17 for 
different Reynolds numbers and jet-to-plate Spacings. The 
PeruessQl at are Shown in Table 10..teA plot of @, versus: L 
1S shown in Figure 10.18..The least square fitted curves for 


different Reynolds numbers are: 


ay eae eeG emt for Re,=100 
ay = 2.00mi 2? for Re,,=200 
sar WES ere = 16 tgs Wea ais -for Re,=300 
aye =e .6e wee tS for Re,=400 Cue 2) 


Equation 10.2 shows that the decay of centerline axial 
velocity in the stagnation flow region for a confined jet 
with initial parabolic profile depends not only on the 
Reynolds number but also on the jet-to-plate spacing. 

The decay of the centerline axial velocity originating 
BEOMaa i lat velocity puomile at, nozalevgexit for driterent 
Reynolds numbers is shown in Figure 10.19 for L=4. The 
centerline axial velocity decreases rapidly near the nozzle 
exit. This is due to the large spreading effects in this 
region. The decrease of centerline axial velocity becomes 
more gentle when the jet flow is further away from the 
nozzle exit. Not until the presence of the impingement plate 


is sensed, the centerline axial velocity does decrease 
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TABLE 10.1 : VALUES OF a, EVALUATED FROM EQUATION 10.1 
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FIGURE 10.18 : a, VERSUS JET-TO-PLATE SPACING FOR THE 
CASE OF PARABOLIC VELOCITY PROFILE 
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FIGURE 10.19 : DECAY OF CENTERLINE AXIAL VELOCITY FOR 
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rapidly to zero at the stagnation point. The decay of the 
centerline) axial velocity for this case is nearly 
independent of the Reynolds number. 

In the stagnation flow region, the centerline axial 
velocity is linearly proportional to the axial distance from 
pene Stagnation” pornt, (L-Y),’ as can be noted from Figure 
16219. The values of a, evaluated from. the. slopes of all 
curves of Figure 10.19 are approximately equal to 0.38 for 


all Reynolds numbers. 


10.1.3 STREAMWISE VELOCITY PROFILE 

The developments of typical streamwise velocity 
profiles with streamwise distance, X, for different Reynolds 
numbers and jet-to-plate spacings from a confined jet with 
initial parabolic velocity profile are shown in Figures 
pOe20-10.3ieeThose from a confined jét.with-initiel.flat 
Velocity profile are shown in Figures 10.32-10.35. The 
dotted line in the plots represents the fully developed 
velocity profile for the particular case. The confinement 
Dlatemis located at (L-Y)/le=1 and thesampingement plate is 
Pocated at; (Lsoy)/L=0. 

For an individual streamwise velocity profile at a 
given streamwise distance, the streamwise velocity, U, 
increases from zero at the impingement plate to a maximum, 


U within a thin layer. Such layer is refered to as the 


max’ 


viscous boundary layer in the stagnation flow region. 


Outside the viscous boundary layer, the streamwise velocity 
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FIGURE 10.20 : STREAMWISE VELOCITY PROFILES FOR L=2 
Rep=100 WITH AN INITIAL PARABOLIC 
VELOCITY PROFILE 
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FIGURE 10.21 : STREAMWISE VELOCITY PROFILES FOR L=2 
Rep=200 WITH AN INITIAL PARABOLIC 
VELOCITY PROFILE 
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FIGURE 10.22 : STREAMWISE VELOCITY PROFILES FOR L=2 
Rep=300 WITH AN INITIAL PARABOLIC 


VELOCITY PROFILE 
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FIGURE 10.23 : STREAMWISE VELOCITY PROFILES FOR L=2 
Rep=400 WITH AN INITIAL PARABOLIC 
VELOCITY PROFILE 
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FIGURE 10.25 : STREAMWISE VELOCITY PROFILES FOR L=4 
Rep=200 WITH AN INITIAL PARABOLIC 
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FIGURE 10.28 : STREAMWISE VELOCITY PROFILES FOR L=12 
Rep=100 WITH AN INITIAL PARABOLIC 
VELOCITY PROFILE 
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FIGURE 10.29 : STREAMWISE VELOCITY PROFILES FOR L=12 
Rep=200 WITH AN INITIAL PARABOLIC 
VELOCITY PROFILE 
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FIGURE 10.30 : STREAMWISE VELOCITY PROFILES FOR L=12 
Re,=300 WITH AN INITIAL PARABOLIC 
VELOCITY PROFILE 
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FIGURE 10.31 : STREAMWISE VELOCITY PROFILES FOR L=12 
Rep=400 WITH AN INITIAL PARABOLIC 
VELOCITY PROFILE 
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FIGURE 10.33 : STREAMWISE VELOCITY PROFILES FOR L=4 
Rep=200 WITH AN INITIAL FLAT 
VELOCITY PROFILE 
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FIGURE 10.34 : STREAMWISE VELOCITY PROFILES FOR L=4 
Rep=300 WITH AN INITIAL FLAT 
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decreases as the axial distance from the impingement plate 
increases. In some streamwise locations, the streamwise 
velocity becomes negative near the confinement plate. These 
negative velocities are caused by the inflow of the primary 
vortex induced by the confinement plate. 
The thickness of the viscous boundary layer in 

Stagnation flow region, o,, is defined as the distance from 
the impingement plate Aboee the streamwise velocity reaches 


S9%. OL pe and can be determined from the streamwise 


x 
velocity profile. The values of co, for=2 and 4 are given 
in Table 10.2. It is found that for a given Reynolds number 
and jet-to-plate spacing, the value of o, remains quite 
constant in the stagnation flow region. The values of so, for 
L=12 cannot be accurately determined due to the coarser grid 
arrangement used near one impingement plate for this case. 
The variation of the maximum value of the streamwise 
velocity at individual streamwise location, U,,,, with 
streamwise distance, %, in the region of XS5 is shown in 
- Figures 10.36-10.38 for the case of parabolic profile and in 
Figure 10.39 for the case of flat profile. In the stagnation 


flow region}. U iso limearly proportiomal “to the distance 


max 
from the stagnation point due to the transformation of the 
axial momentum into the streamwise momentum. A similar 
observation has been made by Schlichting (78) for unconfined 
Submerged jet. Introducing the dimensionless variables in 


Equation 6.6 into Equation 2.6 and choosing U,,,, instead of 


U, yields 
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TABLE 10.2 : THICKNESS OF VISCOUS BOUNDARY LAYER IN 
STAGNATION FLOW REGION 
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FIGURE 10.36 
PARABOLIC VELOCITY PROFILE 
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VARIATION OF Umax WITH STREAMWISE 
DISTANCE FOR L=2 WITH AN INITIAL 
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FIGURE 10.37 : VARIATION OF Umax WITH STREAMWISE 
DISTANCE FOR L=4 WITH AN INITIAL 
PARABOLIC VELOCITY PROFILE 
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FIGURE 10.38 


VARIATION OF Umax WITH STREAMWISE 
DISTANCE FOR L=12 WITH AN INITIAL 
PARABOLIC VELOCITY PROFILE 
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FIGURE 10.39 :; VARIATION OF Umax WITH STREAMWISE 
DISTANCE FOR L=4 WITH AN INITIAL 


FLAT VELOCITY PROFILE 
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The values of a, can be easily evaluated from the slopes of 
peeethe-curves of Figures 10.36-10.39 for different Reynolds 
numbers and jet-to-plate spacings. The values of a, for the 
case of parabolic profile are given in Table 10.3. Comparing 
@iecmvalUeS Of ay) an Table 10.3 to those in Table, 10.1, 
agreement within 2.5% is obtained. Thétivalues of a, for the 
ease of Flat profile aré9approximatély tequal to 0.4 for all 
Reynolds numbers. This value is within 5% agreement with 
ebadreevValuated in Sectzon 10.1.2. 

The thickness of the viscous boundary layer in the 
Stagnation flow region, 9o,, can also be evaluated from a, 


uSing Equation 2.9, 


To /bo= 2880/7 (ax Re.) 95% (10.4) 

The values of o,/b evaluated from Equation 10.4 are compared 
fo those win Table 10.2¢an Figure’ 109407 The dataypoints fall 
fairly closely on the zero error straight line with a slope 
of unity. It is noted that Equation 10.4 was obtained by 
Schlichting (78) in his study of unconfined submerged jet. 
Such a good agreement of the data points indicates that the 
presence of confinement plate has tittle eftect on the stow 
in the stagnation flow region. 

Due to the exchange of momentum with the circulating 


fluid in the primary vortex, U,,,, which increases with X in 
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TABLE 10.3 : VALUES OF a, EVALUATED FROM EQUATION 10.3 


s Bab a 
2 100 TA7e 
200 2o08 
300 Zadol 
400 2S 
4 100 eos 
200 IOS 
300 ch Hie 
400 1 Se a 
12 100 i 
200 LOR 
300 eae 


400 Wez2 


Ty 
on ty a) . 
ih pu 
mee Vs. : 


(hid? Soe cera re a2 Ne 
A eae Neo heat Gin 


158 Foe ciel : 
Rar aah, ee een) a g a ¥ 


195 


9.50 


PARABOLIC PRGFILE o 
PEAT PROP Ce A 


foul gt 
PAN 


2.38/(a, Rep)” 
0.20 0.30 0.40 


0.10 


0.00 


0.00 0.10 0.20 0.30 0.40 0.50 


FIGURE 10.40 : TEST OF GOODNESS FIT FOR EQUATION 10.4 
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the stagnation flow region must eventually reaches a maximum 
and decreases in the wall jet region. The location of the 


maximum of one 1s a function of Reynolds number and 


» 4 
jet-to-plate spacing as can be noted in Figures 10.36-10.39. 
In the wall jet region, flow separation is observed along 
the impingement plate for the case of parabolic profile with 
the exception of the cases of low Reynolds numbers for L=2. 
The locations where the flow starts to separate from the 
impingement plate are given in Table 10.4. For the case of 
flat profile, no flow separation along the impingement plate 
is observed. 


moe dccay.Ot1U aor different, Reynolds numbers min the 


x 


wall jet region are plotted as log(U ) versus log(X) in 


max 
Paguces 10.21-10.43 forja jet withean initial parabolic 
Profile, and in Figure’ 10.44 for a jet with an initial flat 
profile. Obviously, the prediction by Glauert (25) using 
Equation 2.11 in his study of the decay of maximum 
Streamwise velocity in the wall jet region originating from 
an unconfined submerged jet cannot be applied here. For a 
confined air jet, the flow which leaves the stagnation flow 
region is strongly influenced by the presence of the primary 
vortex as shown in Figures 10.1-10.4. Uma, decreases more 
gently after it reaches its maximum value. AS soon as the 
fluid flows pass the centre of the primary vortex, the flow 
expands and U,,,, decreases rapidly. Eventually, 3U a 


approaches the value of that of the fully developed profile 


between two parallel plates. The lower the Reynolds number 
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FIGURE 10.41 : DECAY OF MAXIMUM STREAMWISE VELOCITY FOR 
L=2 WITH AN INITIAL PARABOLIC VELOCITY 
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FIGURE 10.42 : DECAY OF MAXIMUM STREAMWISE VELOCITY FOR 


L=4 WITH AN INITIAL PARABOLIC VELOCITY 
PROFILE 
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: DECAY OF MAXIMUM STREAMWISE VELOCITY FOR 
L=12 WITH AN INITIAL PARABOLIC VELOCITY 
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FIGURE 10.44 : DECAY OF MAXIMUM STREAMWISE VELOCITY FOR 
L=4 WITH AN INITIAL FLAT VELOCITY 
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and the smaller the jet-to-plate spacing, the shorter the 


Streamwise distance needed to approach fully developed flow. 


10.1.4 IMPINGEMENT PLATE SKIN-FRICTION FACTOR 

gre Study Ore local’ skin-friction; factor along the 
impingement plate provides information of the local shear 
SEresco On the plate. The Local skin-friction factor is 
defined by Equation 8.3. For fully developed flow between 
myo parallel pilates, the skin-friction factor can be 
evaluated by setting Y=L in Equation 6.28 and substituting 


into Equation 8.5 to give 


Le) (1005) 


A lot. of C,Re versus X for a given jet-to-plate 


b 
spacing leads to the collapse of all curves for different 
Reynolds numbers to a single curve in the region far away 
from the stagnation flow region where the influence of the 
impinging flow is not sensed. Individual curves for 
different jet-to-plate spacings approach the value of 6/L? 
wMoichewrs equal tew1.5, 0.375 and 0.0497 for L=2) 4 and 12 
respectively. Variations of C, Re, versus X are shown in 
Figures 10.45-10.47 for a jet with an initial parabolic 
velocity profile and for a jet with an initial flat velocity 
profile in Figure 10.48. The general variation pattern of 


the local skin-friction factor is that it increases sharply 


from zero at the stagnation point to a maximum value in a 
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FOR L=2 WITH AN INITIAL PARABOLIC 
VELOCITY PROFILE 


FIGURE 10.45 : VARIATION OF LOCAL SKIN-FRICTION FACTORS 
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VARIATION OF LOCAL SKIN-FRICTION FACTORS 


FOR L=4 WITH AN INITIAL PARABOLIC 
VELOCITY PROFILE 
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FIGURE 10.47 : VARIATION OF LOCAL SKIN-FRICTION FACTORS 
FOR L=12 WITH AN INITIAL PARABOLIC 
VELOCITY PROFILE 
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FIGURE 10.48 : VARIATION OF LOCAL SKIN-FRICTION FACTORS 
FOR L=4 WITH AN INITIAL FLAT 
VELOCITY PROFILE 
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short distance and then decreases with the streamwise 
distance, X. The location of the maximum value is found at 
weve ees LOreboth L=2eand, 4, and at X=1.25 for L=12 for aiget 
with an initial parabolic velocity profile. For a jet with 
an initial flat velocity profile, the maximum value of C,Re, 
Poulocated approximately at X=1.25 for L=4, 

For impinging flow,-the skin-friction factor is 
proportional to the Re nerae number to the power of -0.5. A 
plot ‘of 0.5C,Re,°-® versus X for different Reynolds numbers 
is made so that all the curves will collapse into one 
general curve in the stagnation flow region. Such plots are 
shown in Figures 10.49-10.51 for the case of parabolic 
Prot .e and in Figure 10.52 for thevcase of flat profile. 

For the case of parabolic velocity profile, a general 
curve cannot be obtained in the range of Reynolds numbers 
studied. But there is a trend for L=2 and 4 that a general 
curve may be obtained if the Reynolds number becomes higher. 
For L=12, the value of 0.5C,Re,°-® decreases as Reynolds 


number increases indicating that the normalization of the 


local: skin-friction factor with Re,°-* may be over-corrected 


0 
b.- 
the effect of Reynolds number. This is also probably due to 
the grid effects as mentioned in Chapter 9 for the case of 
Peio eresulting an the inaccuracy of the numerical results 
in the stagnation flow region. The numerical results 
computed by Van Heiningen et. al. (96) for a semi-confined 


two-dimensional jet are also included in Figure 10.50. Their 


results are at Re, =100 and 450 for L=4. There 1s good 


cer entre afl. Ph ay A kip See Aw PA com: 


1) : it 


op bitmai aa ee wR Oth ST eTS eS 
on. nett et usin aia ror ‘ Parry oan. a 
to) & 208 Stee 203 a. tee ae bing i. frets gen Sod: 


| ee a xo ates Ce 
A ) tw sehen name oe eo, qohoolew 


nig hee at ae tem) 4h etoslant aia | 
laine Speight wets wat? vee 
2.0- to ceweg wid ee agdaaan, eblonyst edd ot. 
esadern ab lLowyet: dnaxrethes 204 & boaney . +6638 
eho otmk eaguiios Shire) zavion inet tle: sad 
sta edot@ deye .notpes. wets _potooneets ets ahs 
stiodsieg ho #aep edd 402 ey or-eeor ae | 
aiiteta 34f3, 46 “e802 att 1% tao omght ak . 9: 
tnaenay a . etitorw qiiay ie a tiedalt ag “96 eden oda 
‘Perdine . shitoiniest 30 agit oe rid ‘bantadde oe ni 
ieroetap 9 ads & Dee, Seu 108 ierat 6 ab eed aut 
,ioipiad gomooed ane ae aint api ze Sentai ad’ 
ae ee pape setae bt gira . utes ada ; 
oth 16) God oes rei eto mt ont sa aied dp teeo thar ery %: 
peer nta aa nave Pep eo ce ike “poset 01 ‘ Py atte f ae 
6 9265 ai(g pee vaddag9) RE vena sa an stasis bie ae 
eiluesy TBs iremun ‘ag to soaanadend ‘ait ot ont stuaad 
eq lug ax fests mtn ant 16 . 


Bantisea times a ief 48), (ts i: 


eset 08. Of. sca be ai S 


208 


2 

rN —+>—_+-—_+——+ 

“i Re,=100 

“ @oesesece Re,=200 
----- Rep =300 

© Pe eee Re,=400 

st Rae 

“ Le 


0S 1G, Re, 
0.6 0.8 1.0 ise 


0.4 


0.0 0.2 


0.0 0Q.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5-0 


FIGURE 10.49 : VARIATION OF LOCAL SKIN-FRICTION FACTORS 
NEAR THE STAGNATION FLOW REGION FOR 
L=2 WITH AN INITIAL PARABOLIC 
VELOCITY PROFILE 
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VARIATION OF LOCAL SKIN-FRICTION FACTORS 
NEAR THE STAGNATION FLOW REGION FOR 

L=4 WITH AN INITIAL PARABOLIC 

VELOCITY PROFILE 
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FIGURE 10.52 : VARIATION OF LOCAL SKIN-FRICTION FACTORS 
NEAR THE STAGNATION FLOW REGION FOR 
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agreement between their results and those of this work. 
FOrsrnescace-otjilatrvéelocity profile, the collapse otf 
the curves for different Reynolds numbers to nearly a single 
curve in the stagnation flow region can be noted in Figure 
Moe. rhe curve for Re, =100 deviates the most from the 
general curve probably because the Reynolds number in this 
case 1S Simply not high enough as mentioned above. For 
Re, =100, relatively sioniadeenn retardation of the submerged 
jet is found on the way from the nozzle exit to the 
impingement plate. Also included in Figure 10.52 are the 
theoretical results of Miyazaki and Silberman (58) for an 
unconfined two-dimensional jet. Their results which are 
totally independent of Reynolds number and jet-to-plate 
Spacing (L21.5), are consistently higher than those of this 
work. The disagreement is mainly due to the error introduced 
by their assumption of potential flow outside the viscous 


boundary layer. 


10.2 MASS TRANSFER CHARACTERISTICS 


A measure of local mass transfer is represented by the 
evaluation of local Sherwood number along the impingement 
plate. Local Sherwood numbers along the impingement plate 
are evaluated both experimentally and numerically. 
Experimental and numerical results of local Sherwood number 


are first discussed separately, and finally a comparison of 
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these results is made. 


10.2.1 EXPERIMENTAL RESULTS 

For a confined two-dimensional jet, experimental runs 
for various experimental variables such as the Reynolds 
number, jet-to-plate spacing and duration of the mass 
transfer experiment are made. The range of these 
experimental variables are listed below > 
foe de C=LO=plater spaciand gaigli= 2) pad Mle 
2. Reyolds number: Re, =100, 260, 306, 400 
BemeDuration Of mass transfer experiment: T=120, 180, 240, 

360, 480s 

ie 6 6TYpPGSOL Velocity iprotilé ate nozzle exit: «parabolic 
A listing of the experimental runs together with the 
experimental variables and operating conditions are given in 
Appendix E. 

THeBULLOZeM micinoe, pattern. ror run now o021-3By for 
L=2, Re, = 100 and, T=180ssaseshown inyPlate 10.)la. The “frozen 
PrInge” Sparternetor Fun mow0021-6ASfor l=2., Re = 100 and 
T=360s is shown in Plate 10.1b. These fringes are 
interpreted as contours of equal mass transfer rate. The 
local mass transfer rate can be observed in Plate 10.1 to 
decrease monotonically from the stagnation flow region. It 
will be shown later that only for this case of L=2 and 


Re,=100, a local minimum and a local maximum in local mass 


transfer rate is not observed. 
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PLATE 10.1 : CONTOURS OF EQUAL MASS TRANSFER RATE FOR 


A CONFINED TWO-DIMENSIONAL JET 
(Rep=100, E=2) 
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The "frozen £ninge”™ pattern for run nov:J023-2B- for 
L=2, Re, =306 and T=120s is shown in Plate 10.2a. The "frozen 
Prange” “pat ternwton srin no.4J023-4B for L=27 Re, =306 and 
T=240s is shown in Plate 10.2b. By observing the contours of 
equal mass transfer rate, the presence of a local minimum 
and a local maximum in the mass transfer rate can easily be 
observed in the outer region. The presence of local extrema 
in the mass transfer pes me observed for ali the 
experimental runs in this work except for the case of 
Re, =100 ana L=2 mentioned above. 

IipPLateoO,. 2,710 1S, Creinterest to note that the 
fringes exhibit spanwise fluctuations at the two ends and 
the outer region. The intensity of these fluctuations 
increases with Reynolds number and the fluctuation patterns 
are very much the same for Re ,=306 and 400. This phenomenon 
waS explained by Masliyah and Nguyen (53). The fluctuations 
are attributed to the presence of very slight roughness at 
the edge of the aluminum slot tube. 

Forwusi2eua tyorcaue feozen wiginge” pattern 1S shown: in 
Plate 10232for run ino. 23-23 gwith Re, =306 and T=120s. The 
presence of local extrema in the mass transfer rate 15 again 
observed. But in) this case the fringes are too wide for 
analysis. Invother words, the docal mass transfer rates in 
this outer region are too small to measure. Therefore, no 
quantitative study in the outer region for all the 


experimental runs of L=12 1s made. 


ie 
a es 
ie he 


{mets 


neon pas ; 


216 


PLATE 10.2 : CONTOURS OF EQUAL MASS TRANSFER RATE FOR 
A CONFINED TWO-DIMENSIONAL JET 
(Rep,=306, L=2) 
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PLATE 10.3 : CONTOURS OF EQUAL MASS TRANSFER RATE FOR 
A CONFINED TWO-DIMENSIONAL JET 
(Rep,=306, L= 12) 
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Due to the presence of local extrema in mass transfer 
rate, more complicated "frozen fringe" pattern for a 
confined two-dimensional jet is observed. In such a case, in 
order tc determine the local fringe order, the method of 
zero fringe identification used by Masliyah and Nguyen (51) 
is no longer applicable. Therefore, real time holographic 
interferometry is used. In order to determine the fringe 
order of the "frozen Fringe" pattern, a duplicate 
experimental run is made using real time holography 
interometric technique. 

As soon as the local fringe order is known, local 
Sherwood number, Shp, defined by Equation 5.1 can be easily 
determined by using Equation 5.9. In this work, only the 
local Sherwood numbers on the centerline along the 
Streamwise direction 2n6 measured. For each experimental 
run, the variation of local Sherwood number, Shp, with 
dimensionless streamwise distance, x/b or X, are determined. 
Experimental results of local Sherwood number, Sh,, and 
dimensionless streamwise distance, x/b, are given in 
Appendix E. | 

Variation of local Sherwood number with X are shown in 
Figures 10.53-10.55 for be2, 4 and. V2y in geneva wie 
variation pattern of the local Sherwood number can be 
divided into two regions. The first region is where log(Sh,) 
versus log(X) is linear. The range of such a region is 
depended on the Reynolds number as can be observed in 


Figures 10.53-10.55. Regression analysis is made to obtain a 
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FIGURE 10.53 : VARIATION OF LOCAL SHERWOOD NUMBER WITH 
DIMENSIONLESS STREAMWISE DISTANCE FOR 
L= 2 
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FIGURE 10.54 : VARIATION OF LOCAL SHERWOOD NUMBER WITH 
DIMENSIONLESS STREAMWISE DISTANCE FOR 
Le 4 
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FIGURE 10.55 ; VARIATION OF LOCAL SHERWOOD NUMBER WITH 
DIMENSIONLESS STREAMWISE DISTANCE OF 
Le 12 
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correlation for local Sherwood number with the Reynolds 
number and streamwise distance in this region. The second 
region is where the local extrema in Sherwood number 
occured. For L=12, no local Sherwood numbers are measured in 
this region because the local mass transfer rates are too 
small to measure. For L=2 and 4, although local Sherwood 
numbers are measured in this region, no regression analysis 
is attempted due to the unusual behavior of the local 
Sherwood number with the Reynolds number and streamwise 
eastance in thisaregion.« While the locationeof the local 
maximum in Sherwood number shifts further away from the 
Stagnation point when the Reynolds number increases, the 
value of the local maximum Sherwood number also increases 
with the Reynolds number. 

Two regression equations are obtained from the 
experimental data for L=2 and 4 and for L=12. The range of 
data points used is within the region where log(Sh,) versus 
HOgtX) is linearwand i§ shown in Table 10.5; The regression 
analysis is made in such a way that each data point has 
approximately the same weighting. 


The regression equation for L=2 and 4 is given by: 


= 0.66 = 05 78 (10:26) 
Shi 0.34 Rey X 
for 100sRe,<400, with an average error of 6.8%. A total of 
685 data points are used. A plot to test the applicability 


of Equation 10.6 is shown in Figure 10.56. For a perfect fit 
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TABLE 10.5 


RANGE OF 
EQUATIONS 
Eeb ie 
100 X<6 
200 X<8 
306 9p) 
400 X<10 
100 XS7 
200 Dias) 
306 > oa 
400 KSaz 
100 X< 10 
200 XSi 
306 RZ 


400 


DATA POINTS USED FOR REGRESSIONS 


NO. OF DATA POINTS 


cme ee ee 


oes 


60 


223 


Qt Gil: eres OR | 


Bie EPPO SE a ie ll ee ee A 


Me, 


DORs ARE ae L 


Tee 
, ip 
Rai ty 
fe: Va 
6 By a 
eg 
» oP 
t i 1’ 
Jt 
e- 
aa 


ne) 
0.6 


b 


Sh, /(0.34 Re 


FIGURE 10.56 
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’ TEST OF GOODNESS FIT FOR L = 2 AND 4 
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with zero experimental error, all data points should lie on 
a straight line having a slope of unity. In addition, the 
results obtained by Masliyah and Nguyen (53) presented by 
Equation 5.14 for unconfined two-dimensional jets with L=4 
and 8 are aliso plotted on Figure 10.56, as a dashed line. 
Their data are well within the range of the experimental 
Bccucacy, Of this workjhindicating that. there iss little 
effect of a confinement plate in the wall jet region. 
However, the effect of confinement plate is obvious in the 
outer region where the local extrema of the local Sherwood 
number occur. 

Phemexponenc Ot x) 0.076 ine equation 10.6) iseine fain 
agreement with the experimental findings by Masliyah and 
Nouyen, (53) of -0.73\ and the analytical solutions by 
Schwarz and Caswell (82) of -0.75. Although, the exponenti-or 
Reynolds number, 0.66, in Equation 10.6 does not agree with 
either study, it is within the range of their values of 0.55 
ana: 0.75. 


The regression equation for L=12 is given by: 
Shiesto4 hep ie Rags OVC e) 


for 100<Re,<400, with an average error of 3.8%. A total of 
195 data points are used. A plot to test the applicability 
of Equation 10.7 is shown in Figure 10.57. The exponent of 
X, -0.72, in Equation 10.7 is approximately the same as that 


of L=2 and 4. However, the exponent of Reynolds number, 


i vm + arts sata palaaciat. jane wets 


ye Vate Dee ti alee to pate ihe nea font end + xt 
: ; ae | ~_ ; 


‘ vnc: Sibel “ dywadala ae nalts) pie 


. 4 = eo 
i 


ald at see gat Panini ear nok seupt ni ate ee a 


mioadt bigode avidity ety Oh 


ad “ ie, OF Sa an 9 : 
‘ be treme icy ea) i int. taken. 

bed YEW ate8, jsaaael eros. tis 

yet aeaellh aliee oi de nd be 
Ena Bie Sik ona to ee, a ' dot 


nee boa, 26, Shee ae we pag aes i 
sud tinbs ed iS irawsnktos. 3a: ioatie 3 


he agter 

' we ee 
a. 01 op Ahi a 
an 
aT. eee ee enrianhe 1 A 


ete Togere +y Fane! aise ‘ttwo- 0 


igs se eel, eel nolsnipa Ai 129. 0 sedan 


2.0 ico sae ‘toad 4a oem bail atdity # at + aan 


het 


ty <Gt4 “ “e ig m an ‘eet = = gf 
ee te pee to ee epazeve te doi L003 gome00t 
suk tela ae teat og B0lg A ,bseu ats eiateq ated < 


3 themes abet fect tug nt nvovle at TOP sctniie 


226 


eee 


1.0 


(Ke. =100 
© Rep=200 
4 Rep,p=306 
+ Rep=400 


0.9 
+>PGBR 


0.8 


hae: 


0.7 


0.6 


0.5 


Shp /U.34 “Rep? ?) 
0.4 


0.3 


0.2 


C.0 Q.1 


| 
| 
| 
| 


0.8 0.9 


FIGURE 10.57 : TEST OF GOODNESS FIT FOR L = 12 
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0.33, is less than that for L=2 and 4. This indicates that 
for L=12, local Sherwood number is a weaker function of 


Reynolds number than that for L=2 and 4, 


10.2.2 NUMERICAL RESULTS 

Local Sherwood number, Sh, along the impingement plate 
defined by Equation 5.2 are computed numerically using 
Equation 8.12 for pete ted ete ene number, jet-toe-plate 
Spacing and type of velocity profile at the nozzle exit. The 
M@estangs of all thetnumerical runs for L=2); ¥4.sand 12 are 
given in Appendix G. In this section, only the numerical 
results using the upstream-weighted differencing scheme 
(U.W.D.S.) are presented. 

Thermal entrance solutions with a fully developed 
laminar velocity eraehate for flow between parallel plates 
were obtained by McCuen (46,62) for arbitrarily prescribed 
wall temperature or heat flux. Due to the analogy between 
heat and mass transfer, entrance solutions of mass transfer 
with fully developed laminar flow between parallel plates 
can be obtained. The boundary conditions used in this work 
are similar to those of the fundamental solution of the 
third kind. The solution of the third kind corresponding to 
the coordinate system used in this work is listed in Table 
bo. 6’. 

Therefore, a plot of 2L Sh, versus (x/ 2b) / (Ree Se) (Lora 
given jet-to-plate spacing leads to a collapse of curves for 


the different Reynolds numbers to a general curve in the 


Ves re thes At 
; 7 


Ce gaa, ind bow! 
babes yd Lge eater sahil wie. re awolgm 
te aie wetdaue atohainga auoknst ‘ae gest 

ff ine Pawo eat BY: at ine: gthoctoe 46 cays 

ett ded vied gpd! fenbenus, add vt 
staan ett ets  gabakea chek arity ab. +o xibsse 


1 at ie i aalinniaha 3 wit: ent 


ree 
Pi 


levah gag? pee atiart. hve 3 


sta ist See ie oegusad bi 


we 
44 
ha 
ay 
7 
et 
ene 
aed 
. 
A 
Ps! 
iy 
} 


mow taal: yi Eat sid: oft past peern sand: fie on/tsvegi 
rSslanaas _— ag niet vend tev! ponerse, ietensd exam, Ba 
eases Cr il Daas onan soe ees venkoss, begoteveb 1 
| sgbabnvad ad? Pontarde a 


i wr Qe he Enon: tind 


. t fe fd 


oth fo ies aig tae. Lat alton wert So saodd of satin 


a2 jo biregestie nis ‘eatitd ets Yo avisuloe: sat baa’ 6 
sida? id beeett wi thew aime gi besu wmasteye: osentbres3 
5 Ae toe gels (2o\ 6) eunree’ git 9s to talg)s. iniohedadt” - 
rot. anvubp > segellog 6s 24 alien’ eriosge Sisiqrod- sat bt 
ER seis. inrwnee. in of exe anew on or 


TABLE 10.6 : ENTRANCE SOLUTION OF THE THIRD KIND FOR 
MASS TRANSFER WITH FULLY DEVELOPED 
LAMINAR FLOW BETWEEN PARALLEL PLATES 


(X/2L)/(Re, Sc) 25. Sh: 


os Se ga Dan, nee ate 
G.0025 9i2 50 
072071 6m 259 
OF09'5 1/05 
OF025 55206 
O05 4.902 
OAG75 4.866 
ei 4.861 
Ons 4.861 
Os2o 4.861 


INFINITY 4.861 
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outflow region. The general curve approachs the fully 
developed value of 2L Shy. =4.861 in a manner similar to the 
solution curve obtained from Table 10.6. Such plots are 
shown in Figures 10.58-10.60 for a jet with an initial 
Boravolic profile and in Figure’ 10.61 for a jet with an 
moncial flat. profile. 

Lt 1S noted that in Figures 10.58-10.61, the group 
(X/2L)/(Re,Sc) fawrls’ to ia the local Sherwood number 
for the different Reynolds numbers in the region influenced 
by the impinging jet. It is only -good in the outflow region 
where the flow behavior is similar to those between two 
parallel plates. Furthermore, from Figures 10.58-10.61, the 
values of 2L Sh, approach 4.861 in a Similar manner 
regardless of the jet Reynolds number, jet-to-plate spacing 
Sea the eunitial velocity profile at the nozzle -exit™ 

The effect of Reynolds number on the local Sherwood 
number for different jet-to-plate spacings are shown in 
Figures 10.62-10.64 for the case of parabolic profile and in 
Figure 10.65 for the case of flat profile. The maximum local 
Sherwood number is found to occur at the stagnation point 
for all cases. The local Sherwood number remains quite 
constant in the stagnation flow region directly below the 
jet nozzle (x<0.5), and it then decreases with distance away 
from the stagnation point. 

For all cases with an initial parabolic profile, except 
the case of L=2 and Re,=100, the local Sherwood number is 


found to exhibit a local minimum and a local maximum in the 


| ie yen ist : 4 . 
ne Ia hie ~ agit 1 aye ‘Ab baa 


s-oimr 8 see ease: Bh enaeeaee: co-asiat € {. mn 
boone ip: hd Oh aradneh “ahtonyen | 
spetig 4 wo O48 188 pan pil eles *) i ite ai Sts eee" 

wp t aay Mae sci! * O2 “EE ei. “so svadie 

s oo eo a ste a: est . t out sivoiiverts yrs 


eI daltgeite a i ae 7B Lae doanteges. cry a 


# qa sinter iat , ehieeeet” eblenged. at ert, as" 
absent ‘etd. a9 Reinlind Ytisaiew , P ie, 
wrod hiehge | pth Ce aedioee beaainel to soothale 
+ ona 


ay - 


feogt mumeiken att aves, ‘dart Te, 9289 af? 0d aacary 
| fad ited vanpade, ‘wht ‘ss ‘Aloe 6 baie? ei vodmua 

+s bap cree hia admit “Bouwnovt2 tasof sd aoens tia 

edd woled cliaksa ratpes walt notienpata dt at 
fewG SSegRIeLS dita selene adh gant gt Soe , (2202R) otsson 
vaaiog noisenpsae” 

ti72an8 ‘eres sitecataly tai tae as adie cons tte 308 


‘io 


é 
10 Da Ma a ai mac a oe ty 5 
Rep,=100 
Soe OI0000 Rep=200 
ae at, eae 
: le < tteceeeens, eeeee e€p=4 0 
10 Pecan eecer etsy recur Af Peanes 
a 
SE a, 
SN 
24 ¥ 
10! ® 
§ 4.861 
2 
0 
10 nar 2 2 5 74 S) 2 5 
10 fOr? HOR 10°? POE 107! 
(X/2L)/(Re, Se) 
L=2 


FIGURE 10.58 ; pLrors OF (2L Shy) VERSUS (X/2L)/(Re,Sc) 
FOR L=2 WITH AN INITIAL PARABOLIC PROFILE 
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FIGURE 10.59 ; 
FOR L=4 WITH AN INITIAL PARABOLIC PROFILE 
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FIGURE 10.60 ¢ 
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FIGURE 10.61 : PLOTS OF (2L Shy) VERSUS (X/2L)/(Re, Sc) 
FOR L=4 WITH AN INITIAL FLAT PROFILE 
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L=4 WITH AN INITIAL PARABOLIC PROFILE 
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region further downstream. The locations of the extrema 
points are a function of Reynolds eno and jet-to-plate 
Spacing. Contours of the stream-function, shown in Figures 
10.1-10.3, indicate that the presence of these extrema is 
attributed to fiow recirculation in the region between the 
confinement plate and the impingement plate. In particular, 
Figure 10.66 compares the streamwise location of the 
secondary vortex centre Zhe the streamwise location of the 
minimum of Sherwood number. The data points fall fairly 
euosely On a straight line with a slope ofunity. This 
indicates that there is a strong correlation between the 
streamwise location of the secondary vortex and the local 
minimum of the Sherwood number. 

Local extrema in Sherwood number are also found for all 
Gases.with an initial flat’ profile. Only this time, the 
extrema are not as obvious as those for the case of 
parabolic profile. 

TMhemertectwmot initial jet velocity protile on local 
Sherwood number can be studied from Figures 10.63 and 10.65. 
The local Sherwood number for the case of parabolic profile 
is consistently higher than that for the case of flat 
profile at a given Reynolds number except in the region 
where the local minimum of Sherwood number occurs. 
Significant difference is found in the stagnation flow 
region where the stagnation point Sherwood number for the 
case of parabolic profile is between 1.8 and 2.5 times that 


for the case of flat profile. Such an effect Of Inrerad 
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velocity profile was also observed by numerous investigators 
(77,81,90,96) in their studies of both unconfined and 
semi-confined jets. 

Me vetfect or jet-to-plate spacing onthe local 
Sherwood number for different Reynolds numbers is shown in 
Bigure 10.67 for thécase of parabolic profile. Comparing 
the cases of L=2 and 4, there is no significant difference 
in the local Sherwood ee in the wall. jet region before 
the minimum in Sherwood number occurs. While in the 
stagnation flow region, only slight difference in Sherwood 
number is found for the case of low Reynolds numbers. This 
observation is Similar to those by various investigators 
(19,53,79) studying heat and mass transfer due to an 
unconfined impinging jet for a small jet-to-plate Spacing. 
But. torsthe. case of ae local Sherwood numbers in the 
Stagnation flow and the wall jet regions are consistently 
lower than those for L=2 andi4. This can be” explained in 
that the impingement plate is no longer located inside the 
potential core of the submerged jet, therefore the decayed 
centerline approaching velocity results in a lower mass 


transfer rate in the stagnation flow region. 


10.2.2.1 STAGNATION POINT SHERWOOD NUMBER 

Since the stagnation point Sherwood number, Shy Seeks 
very difficult to measure in the experimental set-up of this 
work, the stagnation point Sherwood numbers computed 


numerically are compared with the results in the literature. 
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FIGURE 10.67 : EFFECT OF JET-TO-PLATE SPACING ON LOCAL 
SHERWOOD NUMBER FOR THE CASE OF 
PARABOLIC VELOCITY PROFILE 
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Comparison of stagnation point Sherwood numbers 
evaluated in this work for L=2 and 4 with those in the 
literature for the case of parabolic profile is shown in 
Figure 10.68. It is noted that the difference between the 
Stagnation point Sherwood number of this work for two 
eeeterent jJét-to-plate Spacings of(L=2eand 4ievhiqher fat 
low Reynolds numbers. As Reynolds number increases, this 
difference diminishs. Also included in Figure 10.68 are the 
experimental results for an unconfined jet given by Sparrow 
and Wong (90), the theoretical results for an unconfined jet 
given by Sparrow and Lee (90) in Equation 2.24 and the 
numerical results for a semi-confined jet computed by Van 
Heiningen et. al. (96). Comparing the stagnation point 
Sherwood number for L=4 of this work with that computed by 
Van Heiningen et. al., the stagnation point Sherwood number 
ef this work is consistently Aigher. This is probably due to 
the assumption of fully developed flow at the outflow region 
of X=16.4 and 42 for Re ,= 100 and 450, respectively by Van 
Heingingen et. al.. Such streamwise locations of X=16.4 and 
42 for Re,=100 and 450, respectively.are very close to the 
locations where the local extrema in Sherwood number occur 
as can be noted in Figure 10.63. 

Comparison of stagnation point Sherwood numbers 
evaluated in this work for L=4 with those in the literature 
for the case of flat profile is shown in Figure RO 2G. 
Included in Figure 10.69 are the theoretical results for an 


unconfined jet given by Miyazaki and Silberman (58) in 
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FIGURE 10.68 ; COMPARISON OF STAGNATION POINT SHERWOOD 
NUMBER WITH LITERATURE FOR THE CASE OF 


PARABOLIC VELOCITY PROFILE (Sc=#2.74) 


ay’ =) 
\ente ; 
i; [ 3 

a dhe 


oe 


z 


= ae © 


trong’ 


10! 
6 
os 5 ° 
Y) 
4 
THIS GNORK «ae 
3} REF 58 eaten tn i cS ? 
2 
0 | 
10 2 3 4 ) 6 Tae Bt eS 
10? 103 


RE 


FIGURE 10.69 +: COMPARISON OF STAGNATION POINT SHERWOOD 
NUMBER WITH LITERATURE FOR THE CASE OF 
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Equation 2.23 and the numerical result for a semi-confined 
jet computed by Van Heiningen et. al. (96). The stagnation 
point Sherwood number predicted by Miyazaki and Silberman is 
consistently higher than that of this work. The disagreement 
1S mainly due to the error introduced by their assumption of 
potential flow outside the viscous boundary layer. On the 
other hand, the disagreement between the results computed by 


Pane Heiningen et. al. and that of this work is less obvious. 


10.2.3 COMPARISON OF EXPERIMENTAL AND NUMERICAL RESULTS 

Comparison of experimental and numerical resuits can 
only be made in the region where X>1. AS mentioned above, 
local mass transfer within the stagnation flow region (X<1) 
is difficult to measure in the experimental set-up of this 
work. Comparison of experimental and numerical results are 
SVOwneinerigures 10:.70-10.81 for L=2, 4.4and 12 with? local 
Sherwood number defined by Equation 5.1. 

For L=2 and 4, excellent agreement is obtained between 
the experimental and numerical results. Furthermore, the 
numerical results confirm the presence of local extrema in 
the Sherwood number. The numerical results from the 
upstream-weighted differencing scheme (U.W.D.S.) gave a 
better prediction as would be expected (69). 

For L=12, disagreements between experimental and 
numerical results are found in the stagnation flow and the 
wall jet regions for the cases of Re,2300. This is probably 


due to the inaccuracy of the numerical results for L=12 in 
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FIGURE 10.72 : COMPARISON OF EXPERIMENTAL AND NUMERICAL 


“RESULTS OF SHERWOOD NUMBER FOR L=2 AND 
Re,=300 
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FIGURE 10.73 : COMPARISON OF EXPERIMENTAL AND NUMERICAL 
RESULTS OF SHERWOOD NUMBER FOR L=2 AND 
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FIGURE 10.76 + COMPARISON OF EXPERIMENTAL AND NUMERICAL 
RESULTS OF SHERWOOD NUMBER FOR L=4 AND 
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FIGURE 10.77 : COMPARISON OF EXPERIMENTAL AND NUMERICAL 
RESULTS OF SHERWOOD NUMBER FOR L=4 AND 
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the stagnation flow region as mentioned in Chapter 9. On the 
other hand, this may be also due to the large jet-to-plate 
Spacing of the experimental set-up in this case. Since the 
Popect. ration ofmthemchannel! is: 0.1295a for: be42' comparing to 
those of 0.0216 and 0.043 for L=2 and 4, respectively, the 
end effect is important and the geometry of the experimental 
set-up is no longer two-dimensional. The higher the Reynolds 
number, the more severe the end effect. This could also be 
the reason why the two-dimensional model fails to predict 
the experimental results for L=12 at high Reynolds numbers 
which are measured along the centerline of a 
three-dimensional channel. To clarify the failure of the 
two-dimensional numerical model in this case, a 
three-dimensional numerical model is recommended to study 
the mass transfer eee er eee due to a confined impinging 
two-dimensional jet when the jet-to-plate spacing is large. 
Although no experimental data for L=12 are obtained in the 
region where the local extrema in Sherwood number occur, the 
presence of such extrema is confirmed by the qualitative 
Study of the contours of equal mass transfer rate shown in 
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11. CONCLUSIONS 


Local Sherwood number along the impingement plate is 
found to exhibit a local minimum and a local maximum in 
the region far away from the stagnation flow region. The 
locations of these extrema are a function of Reynolds 
number and jet-to-plate Spacing. 

Regression equations correlating local Sherwood number 
in terms of Reynolds number and streamwise distance for 
different jet-to-plate spacings are presented in the 
wall jet region up to the region where the local 
extremum of Sherwood number occurs. 

In the wall jet region, for the cases of L=2-and 4, no 
effects of soe ine sollte Spacing are found. 

Effect “or@initral nozzle “exit "velocity profile “on “the 
Stagnation point mass transfer rate is important. The 
Stagnation point mass transfer rate for an initial 
parabolic velocity profile is between 1.8 and 2.5 times 
the *valruesror am inreial ilat velocity prors les 

Mass transfer due to a confined laminar impinging 
two-dimensional jet can be successfully predicted by 
using a two-dimensional numerical model with 


upstream-weighted differencing scheme CUoW Seo eyes 
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12. RECOMMENDATIONS 


A more advanced interferometric technique can be used 
instead of the holographic interferometry employed in 
the present experimental study. A technique, so called 
"speckle interferometry", can be employed by replacing 
the hologram in the holographic interferometry by an 
electronic system. This technique deals with ordinary 
images rather than holographic reconstructions, 
therefore the intermediate production, processing and 
reconstruction of a hologram is avoided. An image of the 
object illuminated by a coherent light, is formed by 
conventional optical methods at a photo-sensitive 
surface of a vidicon tube and can be stored and handled 
electronically. In this way, the development of fringes 
can then be studied quantitatively in real time in front 
of the video monitor. 

A three-dimensional numerical model should be developed 
for the study of mass transfer due to a confined 
impinging two-dimensional jet with large jet-to-plate 
Spacing. This model can then be used to predict the 
experimental results for L=12. 

The two-dimensional numerical model, developed in this 
work, can be improved by using a high order differencing 
scheme derived from local solutions of the differential 


equation. This numerical technique, so called "single 
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cell high order differencing scheme" (SCHOS) is 
understudy by Manohar and Masliyah at the University of 


Alberta. 
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13. NOMENCLATURE 


a dimensionless parameter defined in Equation 
dou 

yee as proportional constants in Equations 2.5 and 
Zed 

A calibration constant, kg/m? 

b slot width, m | 

b, arbitrary stretch constant used to adjust 


the grid transformation in Y-direction 


B- geometric constant of optical set-up 

Cc concentration of swelling agent, kmole/m? 

ems, C5 Cg COchfa cients in Bouetions 2.15216 ee .2 
and 2.22 

Cc dimensionless concentration of swelling 


agent, Melstone PAKS odei 


C, Skinetwicticn. actor s/n.) Paes) 

dq nozzle diameter, m 

d,,d, proportional constants in Equations 2.1 and 
PB) 

D diffusion coefficient, m/s?’ 

€1,€2 proportional constants in Equation 2.2 and 
2.4 

a eigenconstants in Equation 6.29 

f coefficient in Equation 7.50 

f,,f, proportional constants in Equations 2.11 and 
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exterior flux of momentum flux) ome / se 
coefficient in Equation 7.50 

Proportional constants in Equations 2.12 and 
2.14 

dimensionless gradient of concentration at 
impingement plate 

constant value defined in Equation 5.9, m 
Jet=to=prate Spacing, mm 

local mass transfer coefficient defined in 
Equation 5.1, m/s 

local mass transfer coefficient defined in 
Equation 5.2, m/s 

constant value equals to cos, 

parameter in Equation 7.50 

dimensionless jet-to-plate spacing, h/b 
coefficients defined in Equation 7.16 
molecular weight of swelling agent, kg/kmole 
fringe order 

coefficients defined in Equation 7.21 
refractave™ index-ofVglass prism 

refractive index of swollen polymer coating 
number of nodes in X-direction which covers 
half of the slot width 

total number of nodes in X-direction 


total number of nodes in Y-direction 


mass flux, kg/m’s 
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Nussult number for two-dimensional jet 
partial vapor pressure of swelling agent in 
bulk flow, kPa 

partial vapor pressure of swelling at jet 
nozzle exit, kPa 

partial vapor pressure of swelling at 
polymer surface, kPa 

total pressure, kPa 

vapor pressure of swelling agent, kPa 
Prandlt number 

Voulmetnic)itow wrateiof air, m*/s 

radial distance measured from the jet 
centre, m 

displacement or recession of polymer 
coating, m 

molar rate of preduction of swelling agent 
per unit volume, kmole/m’s 

Reynolds number for two-dimensional jet, 

vi b/y 

Reynolds number for axisymmetric jet, vj d/» 
local Reynolds number based on the local 
velocity and distance between adjacent nodes 
in Y-direction defined in Equation 7.46 
local Reynolds number based on the local 
velocity and distance between adjacent nodes 


in X-direction defined in Equation 7.44 


parameter in Equation 8.2 
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Schmidt number, D/» 

Sherwood number, kb/D 

Sherwood number, kd/D 

Sherwood number, k'b/D 

stagnation point Sherwood number 
temperature, °C or K 

duration of mass transfer experiment, s 
Streamwise velocity in X-direction, m/s 
maximum value of streamwise velocity in 
X-direction for an individual streamwise 
velocity profile, m/s 

mean velocity in X-direction in outflow 
region, m/s 

dimensionless streamwise velocity in 
Mrdrreetivon, ~U/ Vi 

maximum value of dimensionless streamwise 
velocity InVX-alrecttonsror an ingsvvduaL 
streamwise velocity profile, ae vi 
axial velocity in Y-direction, m/s 

mean velocity of jet at nozzle exit, m/s 
dimensionless axial velocity in Y-direction, 
v/ Vj 

constant relaxation factor 

Streamwise distance measured from the jet 
centre, m 

dimensionless streamwise coordinate, x/b 


axial distance measured from the jet nozzle 
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GREEK SYMBOLS 


a; 


B, 


Be 


AX 


He 4 
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exit, m 

dimensionless axial coordinate, y/b 
eigenfunctions in Equation 6.29 
coefficients defined in Equations 7.38 and 


7.40 


weighting factor for ee eae term in 
X-direction 

eigenvalues in Equation 6.25 

weighting factor for convective term in 
Y-darection 

incident angle of light path travelling from 
glass prism to polymer coating 

refractive angle of light path travelling 
from glass prism to polymer coating 
weaghtang-factor forediffusiongterm un 
X-direction 

weaghting,.dactor forediffiusiongterm in 
Y-direction 

grid increment in X-direction measured 
between adjacent boundaries of control 
volume 

grid increment in Y-direction measured 


between adjacent boundaries of control 


volume 
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grid increment in x-direction measured 
between adjacent nodes 

grid increment in Y-direction measured 
between adjacent nodes 

maximum absolute different of parameter S 
between two consecutive iterations defined 
in Equation 8.13 

Spread of ‘jet, m 

coefficient defined in Equation 7.17 
wavelength of light, m 

elgenvalues in Equation 6.29 
Viscosityof.air, kg/m s 

kinematic viscosity of air, m?/s 
coefficient defined in Equation 7.17 
density of air, kg/m? 

density of swollen polymer coating, kg/m? 
molar density of gas mixture, kmole/m° 
viscous boundary layer thickness in wall jet 
region, m 

viscous boundary layer thickness in 
stagnation flow region, m 

shear stress at impingement plate, N/m? 
dimensionless variable in Equation 7.1 
stream-function, m’/s 

dimensionless stream-function, ¥/(b ie 
VORP CrLY; us. 


dimensionless vorticity, w(b/ vj) 
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slot width as characteristic length 

bulk flow 

nozzle diameter as characteristic length 
east, west, north and south side boundaries 
of control volume 

fully developed flow 

element number 

at jet nozzle exit 

maximum value 

in outflow region 


at coating surface or impingement plate 
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15. APPENDIX A : PHYSICAL PROPERTIES 


The physical properties used in this work in 
determination of local Sherwood number are given by Masliyah 
and Nouyen (51,52,53))...The physical..data which are 
independent of the operating conditions are given as 


follows: 


i 


density of swollen polymer, p Lo O10: ekigy tv 


Ss 


molecular weight of ethylsalicylate, Mw = 166.17 kg/kmol 


wave length of laser light, IN 6220010 20m 
Other physical properties such as vapor pressure of 
ethylsalicylate, viscosity of air, molar density of the gas 
mixture and diffusion coefficient for ethylsalicylate and 
air are functions of the operating conditions. The 
correlations for evaluating these properties at given 


Operating temperature and presSure are given as follows: 
&. The vanpou.pressure of ethylSalicylate ~PoA me a strong 
function of operating temperature, t, and can be 


evaluated from the following equation (36), 


Toqugs2? & 2s69% ae (29SNon aes) (15.1) 


where t is in Kk. 
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1.e. P°- (at? 21°C and, 93, 9kPa) = 0O800843>k pa 


The molar density of gas mixture, p°, can be 


approximated by using the molar density of air at the 


operating conditions. Applying ideal gas law, the molar 


density of gas mixture can be evaluated from the 


following equation, 
pe PA Or | Et) (a5 zy) 


where P and t are operating pressure and temperature, 
respectively. 


i C3, po? ebater? WC vanduS34.9kPanl= 0.03885akmorym: 


The viscosity of air, mu, is used in’ the evaluation: of 
the Reynolds number for an air jet. This property is 
almost independent of pressure at low pressure, but 
increases with increasing temperature (51). 


i.e. a (at 24°Ceands 934 0kPa9 =<ctuSi7xi0ertikog/mes 


The diffusion coefficient, D, is evaluated by using the 


Lennard-Jones expression for gas pairs of non-polar 


molecules. 
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16. APPENDIX B ;: CALIBRATION OF FLOWMETERS 


Iwo erotameters “are usediin this work. Rotameter A, oa 
Fischer & Porter Rotameter (Tube no. FP-0.25-09-G-6.75/61) 
is used for volumetric air flow rate less than 
1728x190 *m*/s. Rotameter B, a Brook Rotameter,;is used for 
volumetric air flow rate greater than 1 76x Os ms 7s up te 
Beox 104m 7S. 

Rotameter A was calibrated at 22.2°C and 93.45kPa. The 
calibration curve of volumetric air flow rate is shown in 
Figure 16.1 with rotameter reading at the top of the float. 
The calibration curve of volumetric air flow rate of 
Rotameter B was supplied by Masliyah and Nguyen (51) with 
calibration conditions at 22°C and 95,73kPa. This curve is 
shown in Figure 16.2. The rotameter reading is at the bottom 
Of tne efaoat. 

Theivalue of volumetric flow rate from the calibration 
eurve, OW, must be ®corrected/for the actualexperamental 


operating conditions using the following expression 
One Out tes coe CP ee (G-ty) 


where Q is actual experimental volumetric flow rate, Q' is 
volumetric flow rate from calibration curve, t 1S operating 
temperature, t' is calibration temperature, P is operating 


pressure and P' is calibration pressure. All temperatures 
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and pressures in Equation 16.1 are in K and kPa, 
respectively. 

The Reynolds number for an air jet issuing from a 
circle tube with diameter d or a slot tube with width b can 
then be evaluated from the volumetric flow rate, 0, as 


follows: 


£5007 iti oF (1-6, 2) 
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where M is the wetted perimeter of the tube. 
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17. APPENDIX C : LISTINGS OF EXPERIMENTAL RESULTS FOR 


UNCONFINED AXISYMMETRIC JET 


A listing of the experimental runs together with the 
mass transfer duration, operating pressure and operating 
temperature for unconfined axisymmetric jet are given in 
Baeble l/.). Local fringe orders of different experimental 


runs are given in Table 17.2. 


TABLE 17.1 : EXPERIMENTAL RUNS FOR UNCONFINED 
AXISYMMETRIC JET 


Reg RUN NO. TCs) I P (kPa) 

1210 CJ12-0 7 5A 90 20 93.4 
CInZe se 180 2050 O38 a2 
CJ12-6A 360 210 o0i2 

1470 CI14=a25B 90 210 94.4 
CJ14-3A 180 20.0 94.5 
CJ14-6A 360 Zin 94.5 
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TABLE 17.2 : EXPERIMENTAL RESULTS FOR UNCONFINED 
AXISYMMETRIC JET 


RUN NO n/a n RUN NO r/d n 
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SOB ed eo Les 
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18. APPENDIX D : LISTINGS OF EXPERIMENTAL RESULTS FOR 


UNCONFINED TWO-DIMENSIONAL JET 


A listing of the experimental runs together with the 
mass transfer duration, operating pressure and operating 
temperature for unconfined two-dimensional jet are given in 
Table 18.1. Local Sherwood numbers of different experimental 


runs are given in Table 18.2. 


TABLE 18.1 : EXPERIMENTAL RUNS FOR UNCONFINED 
TWO-DIMENSIONAL JET 


Rep RUN NO. TEA So) Er ieC) Pe(kPa) 
94 SJ1-2A 120 210 94.6 
SJ1-4F 240 2ie0 94.3 
SJ1-8E 480 Zi) Sao 
204 Saa- 2C 120 20.0 0236 
SJ2~-4A 240 20.0 94.2 
SJ2-8D 480 210 O36 
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TABLE 18.2 : EXPERIMENTAL RESULTS FOR UNCONFINED 
TWO-DIMENSIONAL JET 


RUN NO. x/b Shy, RUN NO. x/b Sh, 
SJ1-2A Ba13 1.789 
4.68 2.385 
B762 2.846 
2.59 3.578 
2.10 4t1%4 
1.69 A713 
SJ1-4F 22507 0.892 SJ2-4A 29.33 0.962 
$3507 1.189 19.33 14263 
9.56 1.486 14.00 1.604 
7.33 $1282 10.57 1.925 
5.85 2.080 8.67 2.245 
4.95 ONT gealoy) 2.566 
4.14 2eGTe 5.89 PL, 
3.58 elt ky 5.07 3.208 
3.05 3.268 4.36 3.528 
3.82 3.849 
3.53 517 0 
SJ1-8E 1S 0.887 SJ2-8D 23.33 12030 
14.96 1035 pce iy oti 
12.06 12183 16.66 ie epee 
SOs 1.331 14.66 1.9701 
8.64 1.478 13.00 1.618 
7.80 1.626 11.66 1.765 
6.76 Wea 10. 33 IwSi12 
5.94 1.922 9.31 2.059 
Bio 25070 8.57 2.206 
4.80 PIN) 7 a Zn 58 
ihm 2.500 
6.57 2.647 
5.97 2.795 
5.57 2.942 
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19. APPENDIX E ; LISTINGS OF EXPERIMENTAL RESULTS FOR 


CONFINED TWO-DIMENSIONAL JET 


iieeings of the experimental runs together with the 
mass transfer duration, operating pressure and operating 
temperature for confined two-dimensional jet are given in 
mepres 19.1, 19.2 and 19.3. Local Sherwood numbers for 
different experimental runs are given in Tables 19.4, 19.5 


lnaeto.6 for L=2, 4 and 12, respectively. 
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TABLE 19.1 


Pons 25 
we =ZC 
JO021-3A 
JOZI- 38 
J021-4A 
J021-4B 
J021-6A 
J021-6B 
J021-8A 


JOZU-8s 


JO22-2A 
wwe 28 
JOUZZ- SA 
Ogi SB 
J022-4A 
J022-4B 
wUZ2-6A 
J022-6B 
JO022-8A 


JUZe=—o8 


EXPERIMENTAL RUNS FOR CONFINED 
TWO-DIMENSIONAL JET (L=2) 


- 180 


21.0 
Dh 
2A 0 
205-0, 
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20 
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400 


re ee ee es 


JOZSAZA 
J023-23 
J023-3A 
JU23-4A 
J023-4B 
JG23-6A 
J023-6B 
JO023-8A 


J023-8B 


JQ24-2A 
J024-2B 
J024-3A 
J024-3B 
J024-4A 
J024-4B 
J024-6A 
J024-6B 
J024-8A 


J024-8B 


TABLE 19.1 (CONTINUED) 
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TABLE 19.2 : EXPERIMENT RUNS FOR CONFINED 


TWO-DIMENSIONAL JET (L=4) 


Rep 


—— = 


100 


RUN NO. rks) vee) P (kPa) 
J041-2A 120 7g 94.2 
J041-2B 120 2.0 93.6 
J041-3A .180 2020 93.8 
J041-3B 180 21.0 I 7 
J041-4A 240 22.0 94.1 
3041-483 240 20.5 347 
J041-4C 240 21.0 93.5 
J041-6A 360 2.6 93.4 
J041-6B 360 20.8 93.3 
J041-8A 480 20.0 93.9 
J041-8B 480 20.5 93.9 
J042-2A 120 20.0 94.2 
J042-2B 120 2020 94.2 
J042-3A 180 205 93.6 
J042-3¢ 180 201..0 93.6 
J042-4A 240 20.0 94.4 
J042-4B 240 20.0 94.4 
J042-6A 360 20)..5 94.1 
J042-6B 360 2.0 94.0 
J042-8A 480 20.0 93.6 
J042-8B 480 240 O27 
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RUN NO. 


J043-2A 
J043-2B 
J043-3A 
J043-3B 
J043-4A 
J043-4B 
J043-6A 
J043-6B 
J043-8A 


J043-8B 


J044-2A 
J044-2B 
J044-3A 
J044-3B 
J044-4A 
J044-4B 
J044-6A 
J044-6B 
J044-8A 


J044-8B 


TABLE 19.2 
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TABLE 19.3 


RUN NO. 
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EXPERIMENTAL RUNS FOR CONFINED 
TWO-DIMENSIONAL JET (L=12) 
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400 


JVZ3-2A 
323-28 
J123-4A 
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J124-2A 
J124-2B 
J124-4A 
J124-4B 
J124-6A 
J124-6B 
J124-8A 


J124-8B 


TABLE 19.3 
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EXPERIMENTAL RESULTS FOR CONFINED 
TWO-DIMENSIONAL JET (L=2) 
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TABLE 12.4 (CONTINUED) 
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20. APPENDIX F : COMPUTER PROGRAM 


The main program, PROGRAM, together with four 
Seubroucines, ITER1, CALCT, ITER2 and CALC2 are listed in 
this appendix. 

A typical output ‘for numerical run no. 2B100-f0r L=2 
and Re,=100 with an initial parabolic velocity profile are 
also listed. The differencing scheme used in this numerical 


run is the upstream-weighted differencing scheme (U.W.D.S.). 
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21. APPENDIX G : LISTINGS OF NUMERICAL RUNS 


mhevlistings of all numerical runs of which the results 
are used for studying the flow and mass transfer 
characteristics due to a confined laminar impinging 
Pyvo-dimensional jet acre given in Tables 21.1,>21.2 ana 21.3 
hon Le?, 44ang 12, respectively. Numerical runs with an 
initial parabolic velocity profile at the nozzle exit are 
studied for all three different jet-to-plate spacings. 
Numerical runs,iwith an initial flat velocity profile Jat the 


nozzle exit are studied for the case of L=4 only. 
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TABLE 21.1 : NUMERICAL RUNS FOR L=2 


Dee 0.15 
DIFFERENCING NOZZLE EXIT 
Rep RUN NO. nx ny SCHEME PROFILE 
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400 2B400 DD eae U.W.D. 5. PARABOLIC 
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TABLE 212 
O'e25 
RUN NO. nx 
4B1 55 
4F } 55 
4B100 55 
4U100 55 
4F100 55 
4B200 55 
4U200 55 
4F200 55 
4B300 55 
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4B400 35 
4U400 Ss) 
4F400 55 
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RUN NO. nx 
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